- Research article
- Open Access
- Open Peer Review
This article has Open Peer Review reports available.
miR-16 promotes the apoptosis of human cancer cells by targeting FEAT
- Hongwei Liang†2,
- Zheng Fu†2,
- Xueyuan Jiang†2, 4,
- Nan Wang2,
- Feng Wang5,
- Xueliang Wang2,
- Suyang Zhang2,
- Yanbo Wang2,
- Xin Yan3,
- Wen-xian Guan5,
- Chen-Yu Zhang2,
- Ke Zen2,
- Yujing Zhang2Email author,
- Xi Chen2Email author and
- Guangxin Zhou1Email author
© Liang et al.; licensee BioMed Central. 2015
Received: 4 August 2014
Accepted: 21 May 2015
Published: 2 June 2015
Although human cancers have heterogeneous combinations of altered oncogenes, some crucial genes are universally dysregulated in most cancers. One such gene, FEAT (faint expression in normal tissues, aberrant overexpression in tumors), is uniformly overexpressed in a variety of human cancers and plays an important role in tumorigenesis by suppressing apoptosis. However, the precise molecular mechanism through which FEAT is upregulated during tumorigenesis remains largely unknown.
In this study, we used bioinformatic analyses to search for miRNAs that potentially target FEAT. We examined the expression of FEAT protein level by western blotting and miR-16 level by qRT-PCR assay. Cancer cell lines (A549, MCF-7 and Huh-7) with miR-16 upregulation and FEAT silencing were established and the effects on apoptosis of cancer cells in vitro were assessed. Luciferase reporter assay was also performed to investigate the interaction between miR-16 and FEAT.
We identified a specific target site for miR-16 in the 3′-untranslated region (3′-UTR) of FEAT. Consistent with the bioinformatic analyses, we identified an inverse correlation between the miR-16 and FEAT protein levels in lung cancer, breast cancer, and hepatocellular cancer tissues. We then experimentally validated miR-16 as a direct regulator of FEAT using cell transfection and luciferase assays. Finally, we demonstrated that the repression of FEAT by miR-16 promoted the apoptosis of cancer cells.
Our findings provide the first clues regarding the role of miR-16 as a tumor suppressor in cancer cells through the inhibition of FEAT translation.
Although our understanding of the molecular mechanisms of carcinogenesis has greatly improved, this knowledge has not led to the identification and development of effective tools for cancer screening and prevention. In theory, one of the most feasible and promising approaches for cancer screening and prevention is targeting a common oncogene that occurs in most tumors. However, the marked heterogeneity and complexity of human tumors make it difficult to identify commonalities among cancers [1, 2]. Oncogenes that contribute to the development of human cancers are highly variable among different types of cancer and among individual tumors of the same type [1, 2]. Thus, it is still poorly understood whether there are crucial oncogenes that are commonly altered in diverse cancers. Recently, Takahashi et al. investigated a previously unrecognized protein, FEAT (faint expression in normal tissues, aberrant overexpression in tumors), and identified it as a novel, prominent promoter of tumorigenesis . FEAT protein is encoded by METTL13 gene (methyltransferase like 13), and it is aberrantly overexpressed in most human cancers but weakly expressed in normal tissues . Remarkably, transgenic mice that ectopically expressed FEAT spontaneously developed tumors, indicating that the FEAT protein potently drives tumorigenesis in vivo . Gene expression profiling has suggested that FEAT drives receptor tyrosine kinase and hedgehog signaling pathways . However, despite these recent advances in our understanding of the important roles of FEAT in cancer progression, the precise molecular mechanism through which FEAT is upregulated during tumorigenesis remains largely unknown. Further studies are needed to fully elucidate the regulation of FEAT expression in normal and neoplastic tissues.
microRNAs (miRNAs) are a class of endogenously expressed, small non-coding RNAs that play an important role in the regulation of gene expression at the post-transcriptional level [4–6]. Some of these miRNAs have attracted special attention for their involvement in the initiation, progression, and metastasis of human cancers [7, 8]. One particularly well-studied example is the ubiquitously expressed and highly conserved miR-16, one of the first miRNAs that was known to be linked to human malignancies . Evidence indicates that miR-16 can modulate the cell cycle, inhibit cell proliferation, promote cell apoptosis, and suppress tumorigenicity both in vitro and in vivo . These effects can be explained by several targets of miR-16: the anti-apoptotic gene B-cell lymphoma 2 (Bcl-2) ; numerous genes involved in the G1-S transition such as CCND1 (cyclin D1), CCND3 (cyclin D3), CCNE1 (cyclin E1), and CDK6 (cyclin-dependent kinase 6) [12–14]; and genes involved in the Wnt signaling pathway, such as WNT3A (wingless-type MMTV integration site family, member 3A) . Consistently, miR-16 is frequently deleted and/or downregulated in many types of cancer, including chronic lymphocytic leukemia [9, 15], prostate cancer , and lung cancer . Thus, miR-16 is generally thought to be a key tumor-suppressive miRNA that can target numerous oncogenes in various human cancers.
Although the dysregulation of miR-16 and FEAT plays an important role in carcinogenesis, no correlation between FEAT and miR-16 in cancers has been reported. In this study, we hypothesized that FEAT is a target of miR-16. After measuring the expression levels of miR-16 and FEAT in different types of human cancer tissues and paired noncancerous tissues, we detected an inverse correlation between miR-16 and FEAT in human cancers. Furthermore, in this study, we experimentally investigated the direct regulation of FEAT by miR-16 and the biological role of miR-16 targeting FEAT in human cancer cells.
Cells and human tissues
Clinical Features of cancer patients
Lung cancer #1
Lung cancer #2
Lung cancer #3
Squamous cell carcinoma
Breast cancer #1
Invasive Ductal Carcinoma
Breast cancer #2
Invasive Ductal Carcinoma
Breast cancer #3
Invasive Ductal Carcinoma
Hepatocellular cancer #1
Hepatocellular cancer #2
Hepatocellular cancer #3
RNA isolation and quantitative RT-PCR
Total RNA was extracted from the cultured cells and human tissues using TRIzol Reagent (Invitrogen, Carlsbad, CA) according to the manufacturer’s instructions. Assays to quantify miRNAs were performed using Taqman miRNA probes (Applied Biosystems, Foster City, CA) according to the manufacturer’s instructions. Briefly, 1 μg of total RNA was reverse-transcribed to cDNA using AMV reverse transcriptase (TaKaRa, Dalian, China) and a stem-loop RT primer (Applied Biosystems). The reaction conditions were as follows: 16 °C for 30 min, 42 °C for 30 min, and 85 °C for 5 min. Real-time PCR was performed using a TaqMan PCR kit on an Applied Biosystems 7300 Sequence Detection System (Applied Biosystems). The reactions were incubated in a 96-well optical plate at 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. All of the reactions were run in triplicate. After the reaction, the cycle threshold (CT) data were determined using fixed threshold settings, and the mean CT of the triplicate PCRs was determined. A comparative CT method was used to compare each condition with the controls. The relative levels of the miRNAs in the cells and tissues were normalized to U6. The amount of miRNA relative to the internal control U6 was calculated using the 2-△△CT equation, in which △△CT = (CT miRNA - CT U6)target—(CT miRNA—CT U6)control. To quantify the FEAT mRNA, 1 μg of total RNA was reverse-transcribed to cDNA using oligo(dT) and Thermoscript (TaKaRa) in the reaction, which was performed under the following conditions: 42 °C for 60 min and 70 °C for 10 min. Next, real-time PCR was performed using the RT product, SYBER Green Dye (Invitrogen) and specific primers for FEAT and GAPDH. The sequences of the primers were as follows: FEAT (sense): 5′—CTTCACCGAGGTCAGCAGTA-3′; FEAT (antisense): 5′—CTCCATGACTCTAGCCGACA-3′; GAPDH (sense): 5′-GATATTGTTGCCATCAATGAC-3′; and GAPDH (antisense): 5′-TTGATTTTGGAGGGATCTCG-3′. The reactions were incubated at 95 °C for 5 min, followed by 40 cycles of 95 °C for 30 s, 60 °C for 30 s, and 72 °C for 30 s. After the reactions were complete, the CT values were determined by setting a fixed threshold. The relative amount of FEAT mRNA was normalized to GAPDH.
Overexpression and knockdown of miR-16
Synthetic pre-mir-16, anti-mir-16, and scrambled negative control RNAs were purchased from Ambion (Austin, TX, USA). All cells were seeded in 6-well plates or 60-mm dishes, and the cells were transfected with Lipofectamine 2000 (Invitrogen) on the following day, when the cells were approximately 70 % confluent. In each well, equal amounts of pre-mir-16, anti-mir-16, or scrambled negative control RNA were used. The cells were harvested 24 h after transfection for quantitative RT-PCR and Western blotting.
Luciferase reporter assay
To test the direct binding of miR-16 to the target gene FEAT, a luciferase reporter assay was performed as previously described . The entire 3′-untranslated region (3′-UTR) of human FEAT was amplified using PCR with human genomic DNA as a template. The PCR products were inserted into the p-MIR-reporter plasmid (Ambion), and the insertion was confirmed to be correct via sequencing. To test the binding specificity, the sequences that interacted with the miR-16 seed sequence were mutated (from UGCUGCU to ACGACGA), and the mutant FEAT 3′-UTR was inserted into an equivalent luciferase reporter. For luciferase reporter assays, A549, MCF-7, and Huh-7 cells were cultured in 24-well plates, and each well was transfected with 1 μg of firefly luciferase reporter plasmid, 1 μg of a β-galactosidase (β-gal) expression plasmid (Ambion), and equal amounts (100 pmol) of pre-mir-16, anti-mir-16, or the scrambled negative control RNA using Lipofectamine 2000 (Invitrogen). The β-gal plasmid was used as a transfection control. Twenty-four hours post-transfection, the cells were assayed using a luciferase assay kit (Promega, Madison, WI, USA).
Plasmid construction and siRNA interference assay
An siRNA sequence targeting the human FEAT cDNA was purchased from Santa Cruz (sc-88139, Santa Cruz Biotechnology, CA, USA). A scrambled siRNA was included as a negative control. A mammalian expression plasmid encoding the human FEAT open reading frame (pReceiver-M02-FEAT) was purchased from GeneCopoeia (Germantown, MD, USA). An empty plasmid served as a negative control. The FEAT expression plasmid and FEAT siRNA were transfected into A549 cells using Lipofectamine 2000 (Invitrogen) according to the manufacturer’s instructions. Both total RNA and protein were isolated 24 h post-transfection. The FEAT mRNA and protein expression levels were assessed by quantitative RT-PCR and Western blotting.
Protein extraction and western blotting
All cells were rinsed with PBS (pH 7.4) and lysed in RIPA Lysis buffer (Beyotime, China) supplemented with a Protease and Phosphatase Inhibitor Cocktail (Thermo Scientific 78440) on ice for 30 min. The tissue samples were frozen solid with liquid nitrogen, ground into a powder and lysed in RIPA Lysis buffer containing the Protease and Phosphatase Inhibitor Cocktail on ice for 30 min. When necessary, sonication was used to facilitate lysis. Cell lysates or tissue homogenates were centrifuged for 10 min (12000 g, 4 °C). The supernatant was collected, and the protein concentration was calculated using a Pierce BCA protein assay kit (Thermo Scientific, Rockford, IL, USA). The protein levels were analyzed using Western blots with the corresponding antibodies. The protein levels were normalized by probing the same blots with a GAPDH antibody. The antibodies were purchased from the following sources: anti-FEAT (X-20) (sc-101995, Santa Cruz Biotechnology, CA, USA) and anti-GAPDH (sc-365062, Santa Cruz Biotechnology, CA, USA). Protein bands were analyzed using the ImageJ software.
The apoptosis of A549 cells was tested using an Annexin V-FITC/propidium iodide (PI) staining assay. A549 cells were cultured in 12-well plates and transfected with pre-mir-16, anti-mir-16, FEAT siRNA, or the FEAT overexpression plasmid to induce apoptosis. Pre-mir-control, anti-mir-control, control siRNA, and control plasmid served as negative controls. Cells were cultured overnight with both serum-containing complete medium and serum-depleted medium; and the attached and floating cells were then harvested. Flow cytometry analysis of apoptotic cells was performed using an Annexin V-FITC/PI staining kit (BD Biosciences, CA, USA). After washes with cold PBS, the cells were resuspended in binding buffer (100 mM HEPES, pH 7.4, 100 mM NaCl, and 25 mM CaCl2), followed by staining with Annexin V-FITC/PI at room temperature in darkness for 15 min. Apoptotic cells were then evaluated by gating PI and Annexin V-positive cells on a fluorescence-activated cell-sorting (FACS) flow cytometer (BD Biosciences, San Jose, CA). All experiments were performed in triplicate.
Based on miRTarBase, there are now more than 100 experimental validated target genes of miR-16 . the Gene Ontology (GO) classification was performed to gain insights into the biological functions of miR-16 target genes, and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed to detect the potential pathway of miRNA target genes. KEGG pathway database is a recognized and comprehensive database including all kinds of biochemistry pathways . The online based software GENECODIS was utilized in those functional annotation .
All of the Western blotting images are representative of at least three independent experiments. Quantitative RT-PCR, the luciferase reporter assay, and the cell viability and apoptosis assays were performed in triplicate, and each experiment was repeated several times. The data shown are the mean ± SE of at least three independent experiments. The differences were considered statistically significant at p < 0.05 using Student’s t-test.
Upregulation of FEAT protein in human cancer tissues
Identification of conserved miR-16 target sites within the 3′-UTR of FEAT
One important mode of post-transcriptional regulation is the repression of mRNA transcripts by miRNAs. miRNAs are therefore likely to play a biologically relevant role in regulating FEAT expression in cancer. Three computational algorithms, including TargetScan , miRanda , and PicTar , were used in combination to identify potential miRNAs that can target FEAT. Using these approaches, miR-16 was identified as a candidate regulatory miRNA of FEAT. The predicted interaction between miR-16 and the target sites in the FEAT 3′-UTR are illustrated in Fig. 1c. One potential miR-16 target site was found in the 3′-UTR of the FEAT mRNA sequence. The minimum free energy value of this hybrid is −22.0 kcal/mol, which is well within the range of genuine miRNA-target pairs. Moreover, perfect base pairing occurs between the seed region (the core sequence that encompasses the first 2–8 bases of the mature miRNA) and the cognate targets. Furthermore, the miR-16 binding sequence in the FEAT 3′-UTR is highly conserved across species.
Detection of an inverse correlation between the miR-16 and FEAT levels in cancer tissues
miRNAs are generally thought to have expression patterns that are the opposite of those of their targets [4–6]. We next investigated whether miR-16 was inversely correlated with FEAT in cancer tissues. After determining the levels of miR-16 in the same three pairs of lung cancer, breast cancer, and hepatocellular cancer tissues and the corresponding noncancerous tissues, we showed that the miR-16 levels were consistently downregulated in cancer tissues (Fig. 1d). Combining the computational prediction with the detection of inverse correlation between miR-16 and FEAT in vivo, it is quite likely that miR-16 is involved in the post-transcriptional regulation of FEAT.
Validation of FEAT as a direct target of miR-16
To determine whether the negative regulatory effects that miR-16 exerted on FEAT expression were mediated through the binding of miR-16 to the presumed sites in the 3′-UTR of the FEAT mRNA, the full-length FEAT 3′-UTR that contained the sole presumed miR-16 binding site was fused downstream of the firefly luciferase gene in a reporter plasmid. The resulting plasmid was transfected into A549, MCF-7, and Huh-7 cells along with pre-mir-16 or anti-mir-16. As expected, the overexpression of miR-16 resulted in a significant reduction of luciferase reporter activity compared with transfection with pre-scramble control, whereas the inhibition of miR-16 resulted in an increase in reporter activity compared with transfection with anti-scramble control (Fig. 2e). Furthermore, we introduced point mutations into the corresponding complementary sites in the FEAT 3′-UTR to eliminate the predicted miR-16 binding site. This mutated luciferase reporter was unaffected by both the overexpression and knockdown of miR-16 (Fig. 2e). This finding suggests that the binding site strongly contributes to the interaction between miR-16 and FEAT mRNA. In conclusion, our results demonstrate that miR-16 directly recognizes and binds to the 3′-UTR of the FEAT mRNA transcript thereby inhibiting FEAT translation.
miR-16 promotes the apoptosis of cancer cells by regulating FEAT
The significance of overexpressed proteins in cancer is recognized as a potential lead for a variety of diagnostic and therapeutic approaches for cancers. Studies that identify and characterize common oncogenic proteins will hopefully advance molecular-targeted cancer screening and prevention. However, due to the marked heterogeneity and complexity of different types of human cancers, it is rather difficult to identify such proteins that are commonly altered in diverse cancers. Notably, a previously unrecognized protein, FEAT, was recently found to be highly expressed in an unusually wide range of tumors but not expressed in most normal tissues , suggesting that FEAT is a ubiquitous protein that is involved in human cancer. FEAT was originally purified from rat livers as a protein that inhibits nuclear apoptosis in vitro . Ex vivo experiments confirmed that FEAT attenuates apoptotic cell death . Studies have demonstrated that FEAT is highly oncogenic in vivo . However, despite these recent advances in our understanding of the important roles of FEAT in cancer progression, the precise molecular mechanism through which FEAT is upregulated during tumorigenesis remains largely unknown. In this study, we showed that silencing FEAT expression using siRNA could promote cell apoptosis in cancer cells, whereas overexpressing FEAT had an opposite effect, validating its role as an essential oncogenic protein during tumorigenesis. Interestingly, we identified discordance between the FEAT protein and mRNA levels in human lung cancer, breast cancer, and hepatocellular cancer tissues. The results suggest that a post-transcriptional regulation mechanism is involved in FEAT repression. One centrally important mode of post-transcriptional regulation is the repression of mRNA transcripts by miRNAs. Therefore, we searched for miRNAs that could target FEAT and experimentally validated miR-16 as a direct regulator of FEAT. The results identified miR-16 as a novel link between the FEAT regulatory pathway and the pathogenesis of cancer. Considering that miR-16 is highly expressed in normal tissues and frequently deleted and downregulated in many types of cancer tissues, the results also explain, at least in part, why FEAT is aberrantly overexpressed in most human cancers but weakly expressed in normal tissues.
In this study, we further investigated whether the cellular phenotypes especially cell apoptosis were regulated by miR-16 targeting FEAT. We showed that miR-16 could suppress FEAT expression and, in turn, promote apoptosis in cancer cells. The results reveal a critical role for miR-16 as a tumor suppressor and pro-apoptotic molecule in carcinogenesis through the repression of FEAT translation. In fact, miR-16 has been reported to act as a tumor-suppressive miRNA in many cancer types [9–16], and multiple apoptosis-related genes are targeted by miR-16, including BCL-2, CCND1, CCND3, and CCNE1 [9, 13]. We performed KEGG pathway analysis and GO annotation analysis on the experimental validated target genes of miR-16, and the results showed that most of these target genes were indeed anti-apoptotic factors (Additional file 1 Table S1 and S2). An emerging common theme is that multiple targets regulated by a single miRNA can act in concert, rather than individually, to regulate the same biological process, such as apoptosis. The coordinated regulation of many targets by a single miRNA may allow for a prompt cellular response to the apoptosis signals. In this study, it is noted that restoring FEAT expression can successfully attenuate the pro-apoptotic effects of miR-16 on cancer cells, although miR-16 has many other targets. The results suggest that targeting FEAT is a major mechanism by which miR-16 exerts its tumor-suppressive and pro-apoptotic function. Therefore, the modulation of FEAT by miR-16 might explain, at least in part, why the downregulation of miR-16 during carcinogenesis can accelerate cancer progression.
Taken together, this study delineates a novel regulatory network employing miR-16 and FEAT to fine-tune cell apoptosis in lung, breast, and hepatocellular cancer cells. This study may provide a potential novel target for future cancer therapy.
This work was supported by grants from the National Natural Science Foundation of China (Nos. 81102039, 81472508, 81401895 and 31271378) and the Natural Science Foundation of Jiangsu Province (Nos. BK2012014 and BK20141376).
- Stratton MR, Campbell PJ, Futreal PA. The cancer genome. Nature. 2009;458(7239):719–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.View ArticlePubMedGoogle Scholar
- Takahashi A, Tokita H, Takahashi K, Takeoka T, Murayama K, Tomotsune D, et al. A novel potent tumour promoter aberrantly overexpressed in most human cancers. Sci Rep. 2011;1:15.PubMedPubMed CentralGoogle Scholar
- Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.View ArticlePubMedGoogle Scholar
- Ambros V. The functions of animal microRNAs. Nature. 2004;431(7006):350–5.View ArticlePubMedGoogle Scholar
- He L, Hannon GJ. MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004;5(7):522–31.View ArticlePubMedGoogle Scholar
- Esquela-Kerscher A, Slack FJ. Oncomirs-microRNAs with a role in cancer. Nat Rev Cancer. 2006;6(4):259–69.View ArticlePubMedGoogle Scholar
- Calin GA, Croce CM. MicroRNA signatures in human cancers. Nat Rev Cancer. 2006;6(11):857–66.View ArticlePubMedGoogle Scholar
- Calin GA, Dumitru CD, Shimizu M, Bichi R, Zupo S, Noch E, et al. Frequent deletions and down-regulation of micro—RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci U S A. 2002;99(24):15524–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Aqeilan RI, Calin GA, Croce CM. miR-15a and miR-16-1 in cancer: discovery, function and future perspectives. Cell Death Differ. 2010;17(2):215–20.View ArticlePubMedGoogle Scholar
- Cimmino A, Calin GA, Fabbri M, Iorio MV, Ferracin M, Shimizu M, et al. miR-15 and miR-16 induce apoptosis by targeting BCL2. Proceedings of the National Academy of Sciences of the United States of America. 2005;102(39):13944–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Linsley PS, Schelter J, Burchard J, Kibukawa M, Martin MM, Bartz SR, et al. Transcripts targeted by the microRNA-16 family cooperatively regulate cell cycle progression. Mol Cell Biol. 2007;27(6):2240–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu Q, Fu H, Sun F, Zhang H, Tie Y, Zhu J, et al. miR-16 family induces cell cycle arrest by regulating multiple cell cycle genes. Nucleic Acids Res. 2008;36(16):5391–404.View ArticlePubMedPubMed CentralGoogle Scholar
- Bonci D, Coppola V, Musumeci M, Addario A, Giuffrida R, Memeo L, et al. The miR-15a-miR-16-1 cluster controls prostate cancer by targeting multiple oncogenic activities. Nat Med. 2008;14(11):1271–7.View ArticlePubMedGoogle Scholar
- Calin GA, Ferracin M, Cimmino A, Di Leva G, Shimizu M, Wojcik SE, et al. A MicroRNA signature associated with prognosis and progression in chronic lymphocytic leukemia. N Engl J Med. 2005;353(17):1793–801.View ArticlePubMedGoogle Scholar
- Bandi N, Zbinden S, Gugger M, Arnold M, Kocher V, Hasan L, et al. miR-15a and miR-16 are implicated in cell cycle regulation in a Rb-dependent manner and are frequently deleted or down-regulated in non-small cell lung cancer. Cancer Res. 2009;69(13):5553–9.View ArticlePubMedGoogle Scholar
- Chen X, Guo X, Zhang H, Xiang Y, Chen J, Yin Y, et al. Role of miR-143 targeting KRAS in colorectal tumorigenesis. Oncogene. 2009;28(10):1385–92.View ArticlePubMedGoogle Scholar
- Hsu SD, Tseng YT, Shrestha S, Lin YL, Khaleel A, Chou CH, et al. miRTarBase update 2014: an information resource for experimentally validated miRNA-target interactions. Nucleic acids research. 2014;42(Database issue):D78–85.View ArticlePubMedGoogle Scholar
- Altermann E, Klaenhammer TR: PathwayVoyager: pathway mapping using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Bmc Genomics. 2005;6:60.Google Scholar
- Tabas-Madrid D, Nogales-Cadenas R, Pascual-Montano A. GeneCodis3: a non-redundant and modular enrichment analysis tool for functional genomics. Nucleic acids research. 2012;40(Web Server issue):W478–483.View ArticlePubMedPubMed CentralGoogle Scholar
- Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB. Prediction of mammalian microRNA targets. Cell. 2003;115(7):787–98.View ArticlePubMedGoogle Scholar
- John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS. Human MicroRNA targets. PLoS Biol. 2004;2(11):e363.View ArticlePubMedPubMed CentralGoogle Scholar
- Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, et al. Combinatorial microRNA target predictions. Nat Genet. 2005;37(5):495–500.View ArticlePubMedGoogle 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 credited. 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.