Computational prediction and experimental validation of evolutionarily conserved microRNA target genes in bilaterian animals
© Takane et al; licensee BioMed Central Ltd. 2010
Received: 27 October 2009
Accepted: 9 February 2010
Published: 9 February 2010
In many eukaryotes, microRNAs (miRNAs) bind to complementary sites in the 3'-untranslated regions (3'-UTRs) of target messenger RNAs (mRNAs) and regulate their expression at the stage of translation. Recent studies have revealed that many miRNAs are evolutionarily conserved; however, the evolution of their target genes has yet to be systematically characterized. We sought to elucidate a set of conserved miRNA/target-gene pairs and to analyse the mechanism underlying miRNA-mediated gene regulation in the early stage of bilaterian evolution.
Initially, we extracted five evolutionarily conserved miRNAs (let-7, miR-1, miR-124, miR-125/lin-4, and miR-34) among five diverse bilaterian animals. Subsequently, we designed a procedure to predict evolutionarily conserved miRNA/target-gene pairs by introducing orthologous gene information. As a result, we extracted 31 orthologous miRNA/target-gene pairs that were conserved among at least four diverse bilaterian animals; the prediction set showed prominent enrichment of orthologous miRNA/target-gene pairs that were verified experimentally. Approximately 84% of the target genes were regulated by three miRNAs (let-7, miR-1, and miR-124) and their function was classified mainly into the following categories: development, muscle formation, cell adhesion, and gene regulation. We used a reporter gene assay to experimentally verify the downregulation of six candidate pairs (out of six tested pairs) in HeLa cells.
The application of our new method enables the identification of 31 miRNA/target-gene pairs that were expected to have been regulated from the era of the common bilaterian ancestor. The downregulation of all six candidate pairs suggests that orthologous information contributed to the elucidation of the primordial set of genes that has been regulated by miRNAs; it was also an efficient tool for the elimination of false positives from the predicted candidates. In conclusion, our study identified potentially important miRNA-target pairs that were evolutionarily conserved throughout diverse bilaterian animals and that may provide new insights into early-stage miRNA functions.
MicroRNAs (miRNAs) are a class of short (18-25 nucleotides) non-coding RNAs that regulate gene expression posttranscriptionally. Their regulatory potential relies heavily on the recognition of binding sites that are located mainly in the 3'-untranslated regions (3'-UTRs) of target messenger RNAs (mRNAs) . Currently, numerous miRNAs with diverse sequences are being characterized in a wide range of species , suggesting that this small RNA molecule has a major effect on phylogeny. The importance of miRNAs is also suggested from recent research demonstrating that miRNA-guided gene regulation is involved in diverse biological functions, such as cell differentiation, development, carcinogenesis, and tumour suppression [3–6]. For example, phylogenetically conserved miRNAs (e.g., let-7, miR-1, miR-124, and miR-125) are involved in cell differentiation and development [7–10]. In this case, let-7 regulates the expression of RAS proteins known as critical oncogene products . Moreover, miR-34, another evolutionarily conserved miRNA, is a direct downregulator of p53 and is involved in a genetic pathway that promotes cell-cycle progression .
In recent years, more than 700 miRNAs have been identified in humans , and this number is increasing. In a recent report by Friedman et al., the expression of a large number of target genes is predicted to be regulated by miRNAs ; however, relatively few of these have been verified experimentally. To overcome this problem, a series of computational methods has been developed to predict a large number of miRNA targets; e.g., TargetScan , RNAhybrid , MicroTar , PITA , miRanda , and PicTar . Nevertheless, these computational approaches often provide numerous target candidates with a large number of false positives because of the weak complementarity between miRNAs and 3'-UTRs . Recently, a phylogenetic profiling approach has been applied to overcome this limitation. For example, studies of the evolution of orthologous target sites have provided insights into the prediction of efficient miRNA targets . As for miRNAs, many miRNA families are found among various bilaterian animals, suggesting that several miRNAs and their target genes may have co-evolved; however, these features have yet to be systematically characterized.
In this study, we hypothesized that the core regulatory relationship between miRNAs and their target genes were conserved throughout the evolution of bilaterian animals. In addition, by predicting these relationships, we sought to elucidate the core function of miRNAs in the primitive gene-regulatory network of the common bilaterian ancestor. Accordingly, we focused on five miRNAs (let-7, miR-1, miR-124, miR-125/lin-4, and miR-34) that are conserved among bilaterian species (Homo sapiens (H. sapiens), Mus musculus (M. musculus), Gallus gallus (G. gallus), Drosophila melanogaster (D. melanogaster), and Caenorhabditis elegans (C. elegans)) and designed a procedure to extract conserved miRNA/target-gene pairs. We extracted evolutionarily conserved miRNA/target-gene pairs based on hybridization patterns and orthologous information. In addition, we experimentally verified several candidate pairs to support our methodology. Our results suggest a functional role of three major miRNAs (let-7, miR-1, and miR-124) that regulated genes related to development, muscle formation, and cell adhesion. These results suggest a new role for the core function of miRNAs in the primitive gene-regulatory network of the common bilaterian ancestor.
Results and Discussion
Extraction of evolutionarily conserved miRNAs among five bilaterian animals
List of miRNAs conserved among various bilaterian animals
Filtering and enrichment of the evolutionarily conserved miRNA/target-gene pairs
Summary of target gene extraction after each screening step
Ratio of predicted miRNA/target-gene pairs (%)
Ratio of experimentally verified miRNA/target-gene pairs (%)
All possible miRNA/gene pairs
STEP 1 (Free energy and seed region)
STEP 2 (Optimal parameter)
Finally, in step 3, we incorporated orthologous gene information and extracted genes that were evolutionarily conserved among more than four diverse bilaterian animals, including H. sapiens. As a result, the number of predicted miRNA/target-gene pairs was minimized substantially, from 10,356 to 31, using a significantly high Enrichment index of 25.7 (Table 2). The number of predicted miRNA/target-gene pairs was especially high for the three miRNAs let-7 (eight targets), miR-1 (seven targets), and miR-124 (eleven targets) compared with miR-125/lin-4 (three targets) and miR-34 (two targets) (Additional file 2). This suggests that let-7, miR-1, and miR-124 may have played a major role in primordial miRNA gene regulation in the common bilaterian ancestor. To verify the significance of conserved miRNA-target gene pairs, we performed same sequence analysis (from step 1 to step 3) against total 25 species-specific miRNAs (5 miRNAs each from 5 species) as a control experiment supposing that these miRNAs are also conserved in other bilaterians. For example, target prediction of hsa-miR-2277, a species-specific miRNA in human was performed in all 5 species (step 1 and step 2) and conserved targets were extracted (step 3). As a result, 11 out of 25 non-conserved miRNAs did not show any conserved miRNA-target gene pair. Furthermore, average number of the miRNA-target gene pairs of the negative control was 2.4, which is statistically lower than that of conserved miRNA-target gene pairs 6.2 based on the Welch's t-test (P < 0.05). These results support that number of genes achieved from the prediction of conserved miRNAs target genes in this study is indeed significant. In summary, we developed a new filtering method for extracting evolutionarily conserved miRNA/target-gene pairs, which was used to extract 31 reliable miRNA/target-gene pairs among the five families of miRNAs.
Evolutionarily conserved genes regulated by miRNAs
Candidates targeted by multiple miRNAs
Band 4.1-like protein 4B
Candidates targeted by single miRNAs
Lin-28 homolog B
Leucine-rich repeats and immunoglobulin-like domains protein 2 precursor
Inositol-trisphosphate 3-kinase C
Leucine-rich repeats and immunoglobulin-like domains protein 3 precursor
Eukaryotic translation initiation factor 2C 4
Insulin-like growth factor 2 mRNA-binding protein 3
Tropomyosin alpha-3 chain
Alpha-1,3-mannosyl-glycoprotein 4-beta-N-acetylglucosaminyltransferase A
La-related protein 4
Vacuolar ATP synthase subunit B, brain isoform
Laminin subunit gamma-1 precursor
STE20-related adapter protein
ETS domain-containing protein Elk-3
Long-chain-fatty-acid--CoA ligase 1
Coproporphyrinogen III oxidase, mitochondrial precursor
Integrin beta-1 precursor
Leucine-rich repeat flightless-interacting protein 2
Serine/threonine-protein kinase MARK1
Lin-28 homolog A
EPN4_HUMAN Isoform 2 of Q14677 - Homo sapiens
Experimental validation of miRNA target genes
To validate the evolutionarily conserved miRNA/target-gene candidates, we performed transfection and luciferase assays on 6 of the 31 identified evolutionarily conserved miRNA/target-gene candidates. Four of the six candidates (CNN3, LARP4, TAGLN2, and ATP6V1B2) were predicted to be regulated by miR-1, while one candidate (ELK3) was predicted to be targeted by miR-124, and the final candidate (EIF2C4) was predicted to be the target of let-7. We used the well-established downregulation of the delta-like protein 1 precursor (DLL1) gene by miR-34 as a positive control , and the LARP4/let-7 pair, which was extracted up to step 2 in our analysis, was chosen as a non-evolutionarily conserved pair. We subcloned the 3'-UTR sequence downstream from the Renilla luciferase gene (Figure 4A) and co-transfected 100 ng of the 3'-UTR reporter construct into HeLa cells using 5, 20, and 60 pmol of miRNA Mimics (hsa-let-7a, hsa-miR-1, hsa-miR-124, hsa-miR-34a, and miRIDIAN microRNA Hairpin Inhibitor Negative Control #1) (see Methods section). We observed the downregulation of six out of six candidates and of the positive control compared with the negative controls (Figure 4B, a-g). Typical results of the reporter gene assay are shown in Figure 4B for the indicated amounts of miRNAs (5, 20, and 60 pmol). The downregulation of these candidates was significant (P < 0.01), although some of these pairs represented an inhibition of only 30% under the current conditions. Among these candidates, TAGLN2 was previously suggested to be downregulated by miR-1, as assessed by microarray analysis . This feature was recently confirmed using the "pulsed stable isotope labelling with amino acids in cell culture" (pSILAC) method and a reporter gene assay . Regarding the let-7/LARP4 combination, the expression of LARP4 was not downregulated after let-7 transfection, which was supported statistically (Figure 4B, h). These experimental results suggest that our new method has the potential for efficiently extracting reliable miRNA/target-gene pairs and may be effective in the elucidation of the primordial regulatory relationships between miRNAs and their target genes during the early stage of bilaterian evolution.
Possible regulation of evolutionarily conserved miRNA targets in bilaterian animals
To provide further insight into the primary functions of evolutionarily conserved miRNAs (Table 3), we next focused on the functions of the target genes and found that evolutionarily conserved miRNA/target genes could be largely classified into four functional categories: development, differentiation, muscle movement, and gene regulation. First, we describe the function of evolutionarily conserved genes involved in development and differentiation. The laminin subunit gamma-1 precursor (LAMC1) gene, which was possibly regulated by miR-124, is one of the major components of the basement membrane. According to Smyth et al. , null mutation of LAMC1 causes embryonic lethality because of the absence of the basement membrane and failure to differentiate the endoderm. Among other candidates regulated by miR-124, the expression of the leucine-rich repeat flightless-interacting protein 2 (LRRFIP2) gene induces an extra axis in Xenopus laevis embryos . Moreover, the semaphorin-6D precursor (SEMA6D) gene, plays an important role in cardiac morphogenesis during chick embryonic dev which was predicted as a candidate of miR-34 targeting, elopment . The lin-28 gene, which regulates developmental timing in C. elegans, is reportedly controlled by lin-4, as assessed using in vivo experiments . Subsequently, the lin-28 gene was also found to be regulated by miR-125, which is an orthologous miRNA of lin-4, in H. sapiens and M. musculus . Our prediction confirmed the regulation of orthologous lin-28 genes by lin-4/miR-125 miRNA in H. sapiens, M. musculus, and C. elegans and further suggested that a similar regulatory relationship was conserved in G. gallus. Interestingly, our prediction showed that the orthologous let-7/lin-28 pair was also evolutionarily conserved among bilaterians. This feature had been validated experimentally in H. sapiens ; however, in C. elegans, lin-28 is expressed in an early stage of development, while let-7 is expressed in a later stage of development, suggesting that regulation of lin-28 by let-7 is subtle in C. elegans . It would be of interest to analyse the interaction of the let-7/lin-28 pair in other species, such as M. musculus and G. gallus, to understand whether stage-specific expression of let-7 was present in the common ancestor of bilaterian animals or if it is a trait acquired later during evolution, as many of the evolutionarily conserved miRNA/target-gene pairs were related to an essential function involved in differentiation and development.
Next, we focused on the tissue-specific miRNA/target genes. Among the 30 evolutionarily conserved target candidates, approximately one-third were expressed in a tissue-specific manner in humans, according to the BioGPS portal http://biogps.gnf.org. miR-1 is highly expressed in muscle tissues . Here, three candidate genes regulated by miR-1 (i.e., TAGLN2, CNN3, and TPM3) are also expressed in muscle tissues, according to BioGPS. TAGLN2 is a homolog of TAGLN, which encodes an actin-binding protein and is a diagnostic marker of breast and colon carcinoma in humans . The CNN3 gene also encodes an actin-binding protein that represses bone morphogenetic protein (BMP) signalling in chondrocytes, which is important for bone formation . Moreover, the TPM3 gene encodes yet another actin-binding protein that modulates muscle contraction. The other miRNA, miR-124, is expressed in the nervous system . Similarly, one of the evolutionarily conserved target candidates, the glucosamine-fructose-6-phosphate aminotransferase 2 (GFPT2) gene, is expressed in the central nervous system . In addition, the annexin A7 (ANXA7) gene, which is another candidate target of miR-124, is involved in the development of the murine brain . The striking overlap between the tissue specificity of evolutionarily conserved miRNA and that of their target genes suggests that one of the main functions of primordial miRNAs may have been the regulation of genes implicated in the temporary control of the development of muscle and of the nervous system, in a tissue-specific manner.
Finally, we found two interesting candidate genes, EIF2C4 and LARP4, which encode translation-related proteins. It is well accepted that miRNAs are regulators of gene expression, mostly at the translational level . EIF2C4 is also known as Argonaute 4 (AGO4). Although the function of AGO4 is unknown, other Argonaute protein family members are involved in the RNA-induced silencing complex (RISC), which is essential for the miRNA or siRNA pathways. A previous microarray analysis performed in HepG2 cells revealed that the EIF2C4 gene was affected by let-7 . In the present study, we demonstrated for the first time the direct downregulation of EIF2C4 by let-7, as assessed using a reporter gene assay in HeLa cells (Figure 4B, e); therefore, we speculate that negative-feedback regulation of EIF2C4 by let-7 exists in the miRNA pathway. It has been reported that regulation of the AGO1 mRNA, which is a major component of the RISC complex, in the miRNA pathway by miR-168 controls plant development in Arabidopsis thaliana . Another candidate, LARP4, encodes a member of La-motif protein family that controls translational efficiency . We also demonstrated the downregulation of the LARP4 gene via miR-1 using a reporter gene assay (Figure 4B, c), which further supports our contention that some of the evolutionarily conserved miRNAs may play an important role in the regulation of translation by controlling the expression levels of translation factors and by negatively regulating their own miRNA pathway.
We developed a procedure to extract potential evolutionarily conserved miRNA/target-gene pairs based on orthologous gene information from five bilaterian animals, and efficiently extracted 31 evolutionarily conserved miRNA/target-gene pairs from 357,430 pairs. We experimentally validated the downregulation of six candidate pairs (out of six tested pairs) in HeLa cells, which suggests that our method using orthologous information was efficient in extracting evolutionarily conserved miRNA target-gene candidates. Our findings reveal that miRNA target sites were conserved among various species, and demonstrate that especially let-7-, miR-1-, and miR-124-mediated gene regulation may have played an important role throughout evolution, in processes such as development, differentiation, and muscle movement. Moreover, our results indicate that miRNA-mediated translational regulation as well as tissue-specific expression of miRNA/target-gene pairs may have already existed in the common bilaterian ancestor. In conclusion, our study will provide new insights into the early stages of miRNA function.
miRNA and 3'-UTR sequence data
We downloaded 2,404 mature miRNA sequences (885 for H. sapiens, 689 for M. musculus, 520 for G. gallus, 153 for D. melanogaster, and 157 for C. elegans) from the miRBase, version 13.0 http://microrna.sanger.ac.uk/sequences/. We downloaded sequences corresponding to 3'-UTR using the Ensembl transcript ID annotation in FASTA format (40,498 transcripts for H. sapiens, 3,332 for M. musculus, 13,089 for G. gallus, 16,822 for D. melanogaster, and 13,560 for C. elegans) from the Ensmart database Ensembl release 53 http://www.ensembl.org/index.html. Orthologous gene information was also downloaded from the Ensmart database Ensembl release 53. We obtained 145 experimentally verified miRNA/target-gene pairs from TarBase Version 5.0.1 http://diana.cslab.ece.ntua.gr/tarbase/.
Identification of miRNAs conserved among bilaterian animals
The 2,404 miRNA sequences were aligned using ClustalX  with the following alignment parameters: gap opening, 22.50; gap extension, 0.83; and bootstrap value, 100. We checked the conservation of 2,404 miRNA sequences to extract evolutionarily conserved miRNAs. We defined the conservation threshold as an "overall sequence identity > 75% with complete matching of the seed sequence (1-7, 2-8, or 3-9 nucleotides from the miRNA's 5' end)"; furthermore, we introduced information on the phylogenetic relationship among miRNAs to extract reliably conserved miRNAs and used highly conserved miRNA families (category I) .
Extraction of evolutionarily conserved miRNA/target-gene pairs among bilaterian animals
To extract evolutionarily conserved miRNA/target-gene pairs among bilaterian animals, we devised a three-step filtering approach (Figure 1). In step 1, we predicted genes targeted by each of the five miRNAs (let-7, miR-1, miR-124, miR-125/lin-4, and miR-34) using RNAhybrid, which is fast and flexible software for miRNA target prediction, with the free-energy option and the seed-sequence option . The RNA duplex free-energy filter was defined as the appropriate value that led to the efficient extraction of experimentally verified miRNA/target-gene pairs. We also considered a complete match across the seed sequence (1-7, 2-8, or 3-9 nucleotides from the miRNA's 5' end), which was used as a filter by adding the seed option of RNAhybrid.
In step 2, we used four binding parameters of the hybridization pattern of the miRNA/mRNA duplexes. According to a recent study, a binding rule is likely to exist for the recognition of target mRNAs by miRNAs . Moreover, G-U wobble pairs within miRNA/mRNA duplexes play a key role in the interaction with target mRNAs . Subsequently, potential candidates were extracted using four binding parameters (number of mismatches of mRNA within the whole miRNA sequence, number of mismatches of miRNA within the whole miRNA sequence, number of G-U wobble pairs within the whole miRNA sequence, and number of G-U wobble pairs within the seed sequence) of the miRNA/mRNA duplexes (Additional file 1). We used the hybridization pattern of experimentally verified and predicted miRNA/mRNA pairs to calculate coverage, by changing these binding parameters one by one (28-0 for the number of mismatches of mRNA within the whole miRNA sequence; 14-0 for the number of mismatches of miRNA within the whole miRNA sequence; 10-0 for the number of G-U wobble pairs within the whole miRNA sequence; and 6-0 for the number of G-U wobble pairs within the seed sequence) (Figure 2A-D). The range of each of the four binding parameters was determined based on the coverage of experimentally verified and predicted miRNA/mRNA pairs. Parameter combinations were then plotted on a 2D graph by calculating the "ratio of experimentally verified miRNA/mRNA" and "Enrichment" using the points in the four parameter ranges (Figure 2E). The criterion "Enrichment" was defined and calculated as the value of "ratio of experimentally verified miRNA/target-gene pairs" divided by the "ratio of predicted miRNA/target-gene pairs". We obtained the most effective combination of four binding parameters for extracting miRNA/target-gene pairs based on the EC value. Parameter conbination with highest EC value was selected.
In step 3, orthologous gene information was used to extract orthologous genes targeted by the same type of miRNA. For the retrieval of evolutionarily conserved miRNA/target-gene pairs from various bilaterian animals, we set the orthologous gene information criteria as orthologous genes conserved in at least four species, each containing the miRNA target site of interest.
To construct target-site reporter plasmids, each DNA fragment (3'-UTR sequence of the DLL1 gene (668 nt; accession no. AF003522), ELK3 gene (519 nt; accession no. BC017371), EIF2C4 gene (2148 nt; accession no. AB046787), TAGLN2 gene (1391 nt; accession no. D21261), LARP4 gene (1678 nt; accession no. AY004310), CNN3 gene (1391 nt; accession no. BC025372), and ATP6V1B2 gene (1208 nt; accession no. L35249)) was amplified from HeLa genomic DNA via polymerase chain reaction using site-specific primers and was inserted into the Xho I/Not I sites of the psiCHECK-2 plasmid vector (which encodes both firefly and Renilla luciferases; Promega, Madison, WI, USA) (Figure 4A). The oligonucleotide was designed to introduce Xho I and Not I sites at the 5' and 3' termini, respectively. The resulting plasmids were termed pLuc-DLL1, pLuc-CNN3, pLuc-LARP4, pLuc-ELK3, pLuc-EIF2C4, pLuc-TAGLN2, and pLuc-ATP6V1B2, respectively. miRIDIAN™ miRNA Mimic for hsa-let-7a, hsa-miR-1, hsa-miR-124, hsa-miR-34a, and negative control (miRIDIAN microRNA Hairpin Inhibitor Negative Control #1) were purchased from Dharmacon. miRNA Mimic molecules are chemically modified double-stranded RNA oligonucleotides.
The sequences of the oligodeoxyribonucleotides used for PCR were as follows:
(1) DLL1_S: 5'-TTACTCGAGAATGGAAGTGAGATGGCAAGAC-3'
(2) DLLl_A: 5'-TTAGCGGCCGCTTGTCATTCATAAAATTTATTT-3'
(3) CNN3_S: 5'-TTACTCGAGTCCACACAGAAGGAGCTCAG-3'
(4) CNN3_A: 5'-TTAGCGGCCGCCAGGAAGAGCAAATGCATCA-3'
(5) LARP4_S: 5'-TTACTCGAGGTTCCCATTTGATGGCATGT-3'
(6) LARP4_A: 5'-TTAGCGGCCGCCATAGCACCTTGGCGATGTT-3'
(7) ELK3_S: 5'-TTACTCGAGCGTCTGGCCACAATTAAGGA-3'
(8) ELK3_A: 5'-TTAGCGGCCGCTGCTTTCATATTGCCCACTG-3'
(9) EIF2C4_S: 5'-TTACTCGAGTTCTACCAGCAGCTCGGAAT-3'
(10) EIF2C4_A: 5'-TTAGCGGCCGCTTGGATTCCAGCAAGTCCTC-3'
(11) TAGLN2_S: 5'-TTACTCGAGCCCTCCCACGAATGGTTAAT-3'
(12) TAGLN2_A: 5'-TTAGCGGCCGCATGGAAAATGAGAAGCCACG-3'
(13) ATP6V1B2_S: 5'-TTACTCGAGTCCGCGCTCTTGTGAAATAC-3'
(14) ATP6V1B2_A: 5'-TTAGCGGCCGCATAATCATGCTGACTCCCCC-3'
Transfection and luciferase reporter assay
Transient transfection and luciferase assays were performed as described previously, with slight modifications . Briefly, HeLa cells were grown in 10% FBS in DMEM and seeded in 24-well plates 24 h before transfection. Cells were transfected with the indicated amounts of reporter and miRNA Mimic (100 ng of target reporter and 5, 20, and 60 pmol of miRNA Mimic) in the presence of Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA). Firefly and Renilla luciferase activities were measured consecutively using the Dual-luciferase assay system (Promega) 24 h after transfection, according to the manufacturer's instructions.
List of abbreviations used
- EC value:
multiplied value of Enrichment and Coverage
We are grateful to Mizuki Komasa, Motomu Matsui, and Hikaru Taniguchi (Keio University, Japan) for useful discussions. This research was supported in part by a Grant-in-Aid from the 21st Century Centre of Excellence (COE) Program entitled "Understanding and Control of Life's Function via Systems Biology" (Keio University); research funds from the Yamagata Prefectural Government and Tsuruoka City in Japan; and a grant from Keio University.
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.PubMedView Article
- Pasquinelli AE, Reinhart BJ, Slack F, Martindale MQ, Kuroda MI, Maller B, Hayward DC, Ball EE, Degnan B, Muller P: Conservation of the sequence and temporal expression of let-7 heterochronic regulatory RNA. Nature. 2000, 408 (6808): 86-89. 10.1038/35040556.PubMedView Article
- Shivdasani RA: MicroRNAs: regulators of gene expression and cell differentiation. Blood. 2006, 108 (12): 3646-3653. 10.1182/blood-2006-01-030015.PubMed CentralPubMedView Article
- Wienholds E, Plasterk RH: MicroRNA function in animal development. FEBS Lett. 2005, 579 (26): 5911-5922. 10.1016/j.febslet.2005.07.070.PubMedView Article
- He L, Thomson JM, Hemann MT, Hernando-Monge E, Mu D, Goodson S, Powers S, Cordon-Cardo C, Lowe SW, Hannon GJ: A microRNA polycistron as a potential human oncogene. Nature. 2005, 435 (7043): 828-833. 10.1038/nature03552.PubMed CentralPubMedView Article
- Johnson CD, Esquela-Kerscher A, Stefani G, Byrom M, Kelnar K, Ovcharenko D, Wilson M, Wang X, Shelton J, Shingara J: The let-7 microRNA represses cell proliferation pathways in human cells. Cancer Res. 2007, 67 (16): 7713-7722. 10.1158/0008-5472.CAN-07-1083.PubMedView Article
- Reinhart BJ, Slack FJ, Basson M, Pasquinelli AE, Bettinger JC, Rougvie AE, Horvitz HR, Ruvkun G: The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature. 2000, 403 (6772): 901-906. 10.1038/35002607.PubMedView Article
- Chen JF, Mandel EM, Thomson JM, Wu Q, Callis TE, Hammond SM, Conlon FL, Wang DZ: The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation. Nat Genet. 2006, 38 (2): 228-233. 10.1038/ng1725.PubMed CentralPubMedView Article
- Yu JY, Chung KH, Deo M, Thompson RC, Turner DL: MicroRNA miR-124 regulates neurite outgrowth during neuronal differentiation. Exp Cell Res. 2008, 314 (14): 2618-2633. 10.1016/j.yexcr.2008.06.002.PubMed CentralPubMedView Article
- Caygill EE, Johnston LA: Temporal regulation of metamorphic processes in Drosophila by the let-7 and miR-125 heterochronic microRNAs. Curr Biol. 2008, 18 (13): 943-950. 10.1016/j.cub.2008.06.020.PubMed CentralPubMedView Article
- Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ: RAS is regulated by the let-7 microRNA family. Cell. 2005, 120 (5): 635-647. 10.1016/j.cell.2005.01.014.PubMedView Article
- He L, He X, Lim LP, de Stanchina E, Xuan Z, Liang Y, Xue W, Zender L, Magnus J, Ridzon D: A microRNA component of the p53 tumour suppressor network. Nature. 2007, 447 (7148): 1130-1134. 10.1038/nature05939.PubMed CentralPubMedView Article
- Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, D154-158. 36 Database
- Friedman RC, Farh KK, Burge CB, Bartel DP: Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19 (1): 92-105. 10.1101/gr.082701.108.PubMed CentralPubMedView Article
- Rehmsmeier M, Steffen P, Hochsmann M, Giegerich R: Fast and effective prediction of microRNA/target duplexes. RNA. 2004, 10 (10): 1507-1517. 10.1261/rna.5248604.PubMed CentralPubMedView Article
- Thadani R, Tammi MT: MicroTar: predicting microRNA targets from RNA duplexes. BMC Bioinformatics. 2006, 7 (Suppl 5): S20-10.1186/1471-2105-7-S5-S20.PubMed CentralPubMedView Article
- Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E: The role of site accessibility in microRNA target recognition. Nat Genet. 2007, 39 (10): 1278-1284. 10.1038/ng2135.PubMedView Article
- John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human MicroRNA targets. PLoS Biol. 2004, 2 (11): e363-10.1371/journal.pbio.0020363.PubMed CentralPubMedView Article
- Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M: Combinatorial microRNA target predictions. Nat Genet. 2005, 37 (5): 495-500. 10.1038/ng1536.PubMedView Article
- Watanabe Y, Tomita M, Kanai A: Computational methods for microRNA target prediction. Methods Enzymol. 2007, 427: 65-86. 10.1016/S0076-6879(07)27004-1.PubMedView Article
- Gaidatzis D, van Nimwegen E, Hausser J, Zavolan M: Inference of miRNA targets using evolutionary conservation and pathway analysis. BMC Bioinformatics. 2007, 8: 69-10.1186/1471-2105-8-69.PubMed CentralPubMedView Article
- Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115 (7): 787-798. 10.1016/S0092-8674(03)01018-3.PubMedView Article
- Griffiths-Jones S: The microRNA Registry. Nucleic Acids Res. 2004, D109-111. 10.1093/nar/gkh023. 32 Database
- Huang Y, Gu X: A bootstrap based analysis pipeline for efficient classification of phylogenetically related animal miRNAs. BMC Genomics. 2007, 8: 66-10.1186/1471-2164-8-66.PubMed CentralPubMedView Article
- Grimson A, Srivastava M, Fahey B, Woodcroft BJ, Chiang HR, King N, Degnan BM, Rokhsar DS, Bartel DP: Early origins and evolution of microRNAs and Piwi-interacting RNAs in animals. Nature. 2008, 455 (7217): 1193-1197. 10.1038/nature07415.PubMedView Article
- 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 CentralPubMedView Article
- Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM: Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005, 433 (7027): 769-773. 10.1038/nature03315.PubMedView Article
- Selbach M, Schwanhausser 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.PubMedView Article
- Smyth N, Vatansever HS, Murray P, Meyer M, Frie C, Paulsson M, Edgar D: Absence of basement membranes after targeting the LAMC1 gene results in embryonic lethality due to failure of endoderm differentiation. J Cell Biol. 1999, 144 (1): 151-160. 10.1083/jcb.144.1.151.PubMed CentralPubMedView Article
- Liu J, Bang AG, Kintner C, Orth AP, Chanda SK, Ding S, Schultz PG: Identification of the Wnt signaling activator leucine-rich repeat in Flightless interaction protein 2 by a genome-wide functional analysis. Proc Natl Acad Sci USA. 2005, 102 (6): 1927-1932. 10.1073/pnas.0409472102.PubMed CentralPubMedView Article
- Toyofuku T, Zhang H, Kumanogoh A, Takegahara N, Suto F, Kamei J, Aoki K, Yabuki M, Hori M, Fujisawa H: Dual roles of Sema6D in cardiac morphogenesis through region-specific association of its receptor, Plexin-A1, with off-track and vascular endothelial growth factor receptor type 2. Genes Dev. 2004, 18 (4): 435-447. 10.1101/gad.1167304.PubMed CentralPubMedView Article
- Moss EG, Lee RC, Ambros V: The cold shock domain protein LIN-28 controls developmental timing in C. elegans and is regulated by the lin-4 RNA. Cell. 1997, 88 (5): 637-646. 10.1016/S0092-8674(00)81906-6.PubMedView Article
- Wu L, Belasco JG: Micro-RNA regulation of the mammalian lin-28 gene during neuronal differentiation of embryonal carcinoma cells. Mol Cell Biol. 2005, 25 (21): 9198-9208. 10.1128/MCB.25.21.9198-9208.2005.PubMed CentralPubMedView Article
- Kiriakidou M, Nelson PT, Kouranov A, Fitziev P, Bouyioukos C, Mourelatos Z, Hatzigeorgiou A: A combined computational-experimental approach predicts human microRNA targets. Genes Dev. 2004, 18 (10): 1165-1178. 10.1101/gad.1184704.PubMed CentralPubMedView Article
- Shields JM, Rogers-Graham K, Der CJ: Loss of transgelin in breast and colon tumors and in RIE-1 cells by Ras deregulation of gene expression through Raf-independent pathways. J Biol Chem. 2002, 277 (12): 9790-9799. 10.1074/jbc.M110086200.PubMedView Article
- Haag J, Aigner T: Identification of calponin 3 as a novel Smad-binding modulator of BMP signaling expressed in cartilage. Exp Cell Res. 2007, 313 (16): 3386-3394. 10.1016/j.yexcr.2007.08.003.PubMedView Article
- Oki T, Yamazaki K, Kuromitsu J, Okada M, Tanaka I: cDNA cloning and mapping of a novel subtype of glutamine:fructose-6-phosphate amidotransferase (GFAT2) in human and mouse. Genomics. 1999, 57 (2): 227-234. 10.1006/geno.1999.5785.PubMedView Article
- Rick M, Ramos Garrido SI, Herr C, Thal DR, Noegel AA, Clemen CS: Nuclear localization of Annexin A7 during murine brain development. BMC Neurosci. 2005, 6: 25-10.1186/1471-2202-6-25.PubMed CentralPubMedView Article
- Vaucheret H, Vazquez F, Crete P, Bartel DP: The action of ARGONAUTE1 in the miRNA pathway and its regulation by the miRNA pathway are crucial for plant development. Genes Dev. 2004, 18 (10): 1187-1197. 10.1101/gad.1201404.PubMed CentralPubMedView Article
- Bousquet-Antonelli C, Deragon JM: A comprehensive analysis of the La-motif protein superfamily. RNA. 2009, 15 (5): 750-764. 10.1261/rna.1478709.PubMed CentralPubMedView Article
- Kasprzyk A, Keefe D, Smedley D, London D, Spooner W, Melsopp C, Hammond M, Rocca-Serra P, Cox T, Birney E: EnsMart: a generic system for fast and flexible access to biological data. Genome Res. 2004, 14 (1): 160-169. 10.1101/gr.1645104.PubMed CentralPubMedView Article
- Papadopoulos GL, Reczko M, Simossis VA, Sethupathy P, Hatzigeorgiou AG: The database of experimentally supported targets: a functional update of TarBase. Nucleic Acids Res. 2009, D155-158. 10.1093/nar/gkn809. 37 Database
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23 (21): 2947-2948. 10.1093/bioinformatics/btm404.PubMedView Article
- Lai EC: miRNAs: whys and wherefores of miRNA-mediated regulation. Curr Biol. 2005, 15 (12): R458-460. 10.1016/j.cub.2005.06.015.PubMedView Article
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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.