- Research article
- Open Access
Identification and characterization of novel and conserved microRNAs in several tissues of the Chinese rare minnow (Gobiocypris rarus) based on illumina deep sequencing technology
© Hong et al. 2016
- Received: 12 November 2015
- Accepted: 28 March 2016
- Published: 12 April 2016
MicroRNAs (miRNAs), which comprise a large family of endogenous small non-coding RNA molecules, play important roles in the regulation of gene expression in various biological processes. The Chinese rare minnow (Gobiocypris rarus) is a Chinese native fish species and is used extensively as an experimental fish in China; however, relevant biological data, especially miRNA transcriptome data, have not been well documented. To discover conserved and potential novel miRNAs in Chinese rare minnows, a pool of equal amounts of RNA obtained from 6 different adult rare minnow tissues (brain, eye, gill, liver, muscle and heart) was sequenced using illumina deep sequencing technology.
In the present study, 26,930,553 raw reads, representing 2,118,439 unique high-quality reads, were obtained from the pooled small RNA library. Using bioinformatics analysis, 352 conserved and 112 novel Chinese rare minnow miRNAs were first discovered and characterized in this study. Moreover, we found extensive sequence variations (isomiRs) in rare minnow miRNAs, including internal miRNA isomiRs and terminal isomiRs at both the 5′ and 3′ ends and nucleotide variants. Six conserved and 4 novel miRNAs were selected and validated in 6 different adult rare minnow tissues using quantitative real-time PCR (qPCR). The results showed that miR-30a, miR-30b, and Novel-37 are ubiquitously expressed in a variety of tissues. miR-16a, miR-9, miR-125b, miR-34a, and Novel-69 were predominantly expressed in the brain. Novel-115 and Novel-7 were highly expressed in gills, but were relatively weakly expressed in other tissues. These results provided the expression patterns of miRNA genes in Chinese rare minnow. Finally, based on bioinformatics predictions, we mainly found that Novel-94 and Novel-1b-5p were simultaneously targeted to the 3′UTR of Dmrt1, which controls sex determination and/or sexual differentiation in a variety of metazoans at different sites. Novel-29b targeted the 3′UTR of Foxl2, which is involved in the maintenance of ovarian function and the transcriptional regulation of gonadal differentiation-related genes. Novel-62 and Novel-53 targeted the 3′UTR of ERbeta1 and ERbeta2 (which regulate the transcription of target genes), respectively.
Rare minnow is a widely used model for assessing the risk of environmental pollution in China. Identifying and characterizing rare minnow miRNA genes is necessary to discover the biological function of miRNAs and to screen for new molecule biomarkers to assess the risk of environmental pollution in the future.
- Gobiocypris rarus
- Deep sequencing
- Target prediction
MicroRNAs (miRNAs), a class of endogenous, approximately 22-nt-long, small regulatory RNAs , play important gene-regulatory roles in various biological processes [2, 3] including cell differentiation [4–6], proliferation , growth [7, 8], and aging and apoptosis [9, 10]. Since the first miRNA lin-4 was characterized in C. elegans development , thousands of miRNAs have been found in animals and plants using genetic methods and through the sequencing of small RNA libraries [12, 13].
These small molecules can reduce the production of protein by modulating the stability and/or translational potential of their mRNA targets at the post-transcriptional level . Generally, miRNAs are thought to function mainly by binding to target mRNAs through imperfect base paring with the 3′ untranslated regions (UTRs) and recruit the RNA-induced silencing complex (RISC), finally leading to down-regulation of their target genes [1, 14, 15]. Currently, a large number of miRNAs have been found in a variety of organisms from worms to humans, suggesting the evolutionary conservation of miRNA regulation mechanisms [1, 9, 11]. For instance, the let-7 miRNA gene, which was initially identified as a significant regulator that is involved in the heterochronic pathway controlling developmental timing in C. elegans, was the first miRNA known to be well-conserved from nematodes to primates [16, 17]. It is estimated that more than 60 % of mRNAs exhibit conserved miRNA-binding sites in mammals . However, teleost miRNAs were first reported in zebrafish , and various miRNAs were found to play a role during zebrafish development; some of these functions have been characterized [19, 20]. Giraldez et al.  reported that miR-430 facilitates the deadenylation and clearance of maternal mRNAs during early zebrafish embryogenesis , and numerous studies have demonstrated the potential role of miRNAs in the regulation of tumorigenesis . For example, miR-125b is an important negative regulator of p53 and p53-induced apoptosis during development and the stress response in zebrafish and humans . Furthermore, through the combined use of miRNA microarray platforms and bioinformatics analysis, Craig et al.  identified changes in miRNA expression and related biological pathways in zebrafish following 7 days of exposure to fluoxetine (FLX). Furthermore, another study evaluated miRNAs in relation to toxin/chemical exposure in fish. Brzuzan et al. used qPCR to identify changes in the expression of 9 selected miRNAs in liver samples of whitefish following exposure to microcystin-LR (MC-LR) at a dose of 100 μg·kg−1 body weight for 24 or 48 h. Among these miRNAs, miR-122 exhibited the highest constitutive level, and this was associated with a decrease in the expression of genes coding for ferritin and a Ras-like protein .
The rare minnow (Gobiocypris rarus), a small freshwater cyprinid fish, is an endemic model organism that is mostly distributed in the upstream region of the Yangtze River, Sichuan Province, China. The fish are small (30–80 mm in length), are easy to culture in the laboratory, and have a relatively short life cycle, spawning hundreds of eggs with high fertilization and hatching rates [27, 28]. To date, this Cypriniformes species has been used extensively for various biological studies in China [29, 30]. Additionally, rare minnows have been widely applied in aquatic toxicity testing as a model vertebrate species for many years in our laboratory [27, 28, 30–34]. In the field of toxicology studies, several fish models are currently in use; however, research aimed at investigating miRNA alterations is mainly conducted in vitro, and in vivo studies remain in infancy . The discovery of miRNA genes in the rare minnow genome will contribute to a better understanding of the roles played by miRNAs in regulating diverse biological processes in fish. Such studies are expected to prove useful for the future screening of novel molecule biomarkers to assess the risk of environmental pollution and in vivo toxicity.
All fish were anesthetized with 2-phenoxyethanol before being euthanized. The fish were cared for in accordance with the Regulations for the Administration of Affairs Concerning Experimental Animals for the Science and Technology Bureau of China throughout the study.
The rare minnow has been maintained in our laboratory for over 13 years. The brood stock was maintained in a flow-through system based on dechlorinated tap water and was subjected to a 16:8 h light:dark cycle at 25 ± 1 °C . The brood stock was fed with newly hatched brine shrimp (Artemia nauplii) twice daily and granulated food (TetraMin, Tetra Werke, Melle, Germany) once daily. In this study, ten healthy adult rare minnows (5 months old) were selected for study. The average body weight and body length were 1.2 ± 0.4 g and 42.3 ± 3.7 mm, respectively.
Tissue samples including the brain, eye, gill, liver, muscle and heart were collected from adult rare minnows. All samples were immersed in liquid nitrogen immediately after collection and stored at −80 °C. Total RNA was isolated from the tissue samples using a mirVana™ microRNA Isolation Kit (Ambion, USA), according to the manufacturer’s instructions. RNA integrity was examined using the Agilent 2100 Bioanalyzer system (Santa Clara, CA, USA); the samples were then stored at −80 °C until analysis.
Small RNA library construction and sequencing
In this study, one miRNA library was constructed, and an equal quantity (10 μg) of total RNA was extracted from each of the various tissue samples. Briefly, 15–35 nt molecules were purified from total RNA (mirVana™ microRNA Isolation Kit, Ambion, USA) using polyacrylamide gel electrophoresis (PAGE); then, 5′ and 3′ adaptors (TruSeq® Small RNA Sample Preparation Kit, Illumina, USA) were ligated to the 5′ and 3′ termini of the RNAs, and the samples were used as templates for cDNA synthesis. Subsequently, PCR amplification was performed using primers that were complementary to both adaptors. After purification of the amplified cDNA constructs, the products were sequenced using HiSeq2500 technology. All sequencing reads were deposited in the National Center for Biotechnology Short Read Archive (SRA) database (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/sra/) under the study accession SRP057175.
The original reads obtained from small RNA sequencing were processed by summarizing data production, evaluating sequencing quality and calculating the length distribution of the small RNA reads. The quality assessment and processing of sequenced reads were performed as recommended [36–38]. Briefly, the quality of the sequenced reads was assessed using the FastQC application (FastQC-0.11.2, http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc/). Adapter sequences from the reads were removed with FASTX-toolkit (FASTX-toolkit 0.0.14, http://hannonlab.cshl.edu/fastx_toolkit/commandline.html). After adaptor trimming, filtering low quality tags and removing reads either < 18 nt or > 30 nt, the remaining reads were mapped to the zebrafish genome with a tolerance of one mismatch in the seed sequence. Subsequently, using BLAST against the Rfam (http://rfam.sanger.ac.uk/) database, the reads that were mapped to the zebrafish genome were analyzed to annotate rRNA, tRNA, snRNA, snoRNA and other ncRNA sequences in the small sequences. Conserved rare minnow miRNAs in the sequences were then identified using BLAST searching against the miRNA database, miRBase (release 21, http://www.mirbase.org/) and no more than two mismatches can be allowed to identify homologs of known miRNAs. The sequences that did not match known miRNAs were mapped to the Ensembl (http://www.ensembl.org) zebrafish genome (Zv9) and rare minnow EST and GSS databases to identify potentially novel miRNA candidates; these were identified by folding the flanking genome sequence of unique small RNAs using MIREAP (https://sourceforge.net/projects/mireap/).
Validation of miRNAs by qPCR analysis
qPCR was used to validate six conserved and four novel selected miRNAs, and the relative expression levels were analyzed in six tissues from adult rare minnows, including the brain, liver, heart, muscle, eye and gill. Total RNA was isolated using the miRcute miRNA Isolation Kit (TIANGEN, Beijing, China) according to the manufacturer’s instructions. Then, two micrograms of total RNA from each sample were reverse-transcribed into cDNA using the miRcute miRNA First-Strand cDNA Synthesis Kit (TIANGEN, Beijing, China). After this was completed, qPCR was performed using the miRcute miRNA qPCR Detection Kit (TIANGEN, Beijing, China). The following reaction solution was prepared on ice: 10 μL 2 × miRcute miRNA Premix (containing SYBR & ROX), 0.4 μL PCR forward primer (10 μM), 0.4 μL Uni-miR qPCR Primer (10 μM, part of the Kit), 1.6 μL cDNA, and RNase-free dH2O up to 20 μL. The reaction mixtures were incubated in a 96-well plate at 94 °C for 2 min followed by 40 cycles of 94 °C for 20 s and 60 °C for 34 s. All primers used in the reverse transcription and qPCR experiments are shown in Additional file 1: Table S1. Samples (n = 3) were run simultaneously for each gene in triplicate, and a non-template control was included in each plate. The relative miRNA expression level was calculated using the 2−ΔΔCt method after the threshold cycle (Ct); the 5S rRNA (submitted to the European Nucleotide Archive, http://www.ebi.ac.uk/ena/data/view/LN878980) was used as an endogenous control. miRNA expression levels are presented as 2−ΔΔCt mean ± SE (standard error), and error bars indicate the standard error of 2−ΔΔCt mean values.
Identification of Novel miRNA targets via computational analysis
Target prediction was performed by integrating four miRNA target prediction software packages, including miRanda (http://www.microrna.org/microrna/home.do) , PITA (http://genie.weizmann.ac.il/pubs/mir07/mir07_exe.html) , RNA22 (https://cm.jefferson.edu/rna22/)  and TargetScan (http://www.targetscan.org/fish_62/) . To obtain a better miRNA target analysis result, we used the local version of these four software packages. All newly identified miRNA sequences were used to query the ten gene sequences of rare minnows using the default parameters.
Illumina sequencing of small RNAs
Read statistics of the obtained small RNAs
High-quality reads (> =18 nt)
Aligned to the zebrafish genome
Conserved miRNAs in Metazoa
Non-conserved in Metazoa
Reads for novel miRNA prediction
In summary, the total number of sequence reads generated by high-throughput sequencing using the Illumina platform was of the same order of magnitude (107) as previous results for zebrafish ; however, this number was significantly different (106) from that obtained for Nile tilapia . These differences are due to the evolutionary relationships between the species; zebrafish and rare minnows are both cyprinids, whereas Nile tilapia belongs to cichlidae.
Conserved and Novel miRNAs in rare minnow
To identify conserved miRNAs in rare minnows, perfectly matched reads from the small RNA library were searched against Metazoa mature miRNAs in miRBase (Release 21). Alignment of the sequences to the miRBase database revealed that 732,306 of all sequences yielded a positive match to known miRNAs. Finally, we identified 352 conserved miRNAs (Additional file 2: Table S2) in rare minnow.
The 10 most frequent highly expressed novel miRNA candidates
Most expressed sequence
Abundant miRNAs play essential and broad regulatory function in biological processes. In the current study, the most highly expressed miRNA in rare minnows was miR-215 (51,777 reads). Previous studies have indicated that miR-215 is involved in a variety of cancers in mammals. For example, Deng et al. found that miR-215 is significantly up-regulated in gastric cancer tissues from either gastrectomy or gastroscopy by targeting retinoblastoma tumor-suppressor gene 1 (RB1) through its 3′-UTR in gastric cancer cells . Moreover, White et al. demonstrated that miR-215 overexpression decreased cellular migration and invasion in a renal cell carcinoma (RCC) cell line model using miRNA microarray and qPCR analyses .
Multiple isomiRs in the rare minnow
Validation and analysis of miRNA expression
Identification of Novel miRNA targets via computational analysis
This work represents an initial study on the miRNA profile of the Chinese rare minnow. The differential expression of miRNAs and the prediction of their target genes provide a basis for further understanding rare minnow miRNAs and the biological processes in which they are involved. The rare minnow is a model species for various biological studies in lower vertebrates in China. The identification and functional characterization of rare minnow miRNAs reported here will provide new opportunities for functional genome research on the rare minnow and should prove useful when screening for novel molecule biomarkers for use in assessing the risk of environmental pollution.
This work was supported by the Key Program of the National Natural Science Foundation of China (21437006), the Major International Joint Research Project of the National Natural Science Foundation of China (51420105012), and the Water Pollution Control and Treatment of the National Science and Technology Major Project (2014zx-07204-008-003).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Ambros V. The functions of animal microRNAs. Nature. 2004;431:350–5.View ArticlePubMedGoogle Scholar
- Starner T, Mann S, Rhodes B, Levine J, Healey J, Kirsch D, Picard RW, Pentland A. Augmented reality through wearable computing. Presence Teleoperators Virtual Environ. 1997;6:386–98.Google Scholar
- Singh SK, Pal Bhadra M, Girschick HJ, Bhadra U. MicroRNAs--micro in size but macro in function. FEBS J. 2008;275:4929–44.View ArticlePubMedGoogle Scholar
- Miska EA. How microRNAs control cell division, differentiation and death. Curr Opin Genet Dev. 2005;15:563–8.View ArticlePubMedGoogle Scholar
- Shivdasani R. MicroRNAs: regulators of gene expression and cell differentiation. Blood. 2006;108:3646–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Brennecke J, Hipfner DR, Stark A, Russell RB, Cohen SM. Bantam encodes a developmentally regulated microRNA that controls cell proliferation and regulates the proapoptotic gene hid in Drosophila. Cell. 2003;113:25–36.View ArticlePubMedGoogle Scholar
- Hwang H-W, Mendell JT. MicroRNAs in cell proliferation, cell death, and tumorigenesis. Br J Cancer. 2006;94:776–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Yunqing L, Guessous F, Ying Z, DiPierro C, Kefas B, Johnson E, Marcinkiewicz L, Jinmai J, Yanzhi Y, Schmittgen TD, Lopes B, Schiff D, Purow B, Abounader R. MicroRNA-34a inhibits glioblastoma growth by targeting multiple oncogenes. Cancer Res. 2009;69:7569–76.View ArticleGoogle Scholar
- Somel M, Guo S, Fu N, Yan Z, Hu HY, Xu Y, Yuan Y, Ning Z, Hu Y, Menzel C, Hu H, Lachmann M, Zeng R, Chen W, Khaitovich P. MicroRNA, mRNA, and protein expression link development and aging in human and macaque brain. Genome Res. 2010;20:1207–18.Google Scholar
- Xu P, Guo M, Hay BA. MicroRNAs and the regulation of cell death. Trends Genet. 2004;20:617–24.View ArticlePubMedGoogle Scholar
- Lee RC, Feinbaum RL, Ambros V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993;75:843–54.View ArticlePubMedGoogle Scholar
- Ramachandra RK, Salem M, Gahr S, Rexroad CE, Yao J. Cloning and characterization of microRNAs from rainbow trout (Oncorhynchus mykiss): their expression during early embryonic development. BMC Dev Biol. 2008;8:41–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang R, Marshall D, Bryan GJ, Hornyik C. Identification and characterization of miRNA transcriptome in potato by high-throughput sequencing. PLoS One. 2013;8:e57233.View ArticlePubMedPubMed CentralGoogle Scholar
- Lim LP. Vertebrate MicroRNA genes. Science. 2003;299:1540.View ArticlePubMedGoogle Scholar
- Ason B, Darnell DK, Wittbrodt B, Berezikov E, Kloosterman WP, Wittbrodt J, Antin PB, Plasterk RH a. Differences in vertebrate microRNA expression. Proc Natl Acad Sci U S A. 2006;103:14385–9.Google Scholar
- Pasquinelli AE, Reinhart BJ, Slack F, Martindale MQ, Kuroda MI, Maller B, Hayward DC, Ball EE, Degnan B, Müller P, Spring J, Srinivasan a, Fishman M, Finnerty J, Corbo J, Levine M, Leahy P, Davidson E, Ruvkun G. Conservation of the sequence and temporal expression of let-7 heterochronic regulatory RNA. Nature. 2000;408(November):86–9.Google Scholar
- Roush S, Slack FJ. The let-7 family of microRNAs. Trends Cell Biol. 2008;18(September):505–16.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.View ArticlePubMedPubMed CentralGoogle Scholar
- Wienholds E, Kloosterman WP, Miska E, Alvarez-Saavedra E, Berezikov E, de Bruijn E, Horvitz HR, Kauppinen S, Plasterk RHA. MicroRNA expression in Zebrafish embryonic development. Science (80-). 2005;309:310–1.Google Scholar
- Giraldez AJ, Cinalli RM, Glasner ME, Enright AJ, Thomson JM, Baskerville S, Hammond SM, Bartel DP, Schier AF. MicroRNAs regulate brain morphogenesis in Zebrafish. Science (80-). 2005;308:833–8.Google Scholar
- Giraldez AJ, Mishima Y, Rihel J, Grocock RJ, Van Dongen S, Inoue K, Enright AJ, Schier AF. Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science. 2006;312:75–9.Google Scholar
- Lee MT, Bonneau AR, Giraldez AJ. Zygotic genome activation during the maternal-to-zygotic transition. Annu Rev Cell Dev Biol. 2014;30(August):581–613.View ArticlePubMedPubMed CentralGoogle Scholar
- Di Leva G, Garofalo M, Croce CM. MicroRNAs in cancer. Annu Rev Pathol Mech Dis. 2014;9:287–314.View ArticleGoogle Scholar
- Le MTN, Teh C, Shyh-Chang N, Xie H, Zhou B, Korzh V, Lodish HF, Lim B. MicroRNA-125b is a novel negative regulator of p53. Genes Dev. 2009;23:862–76.Google Scholar
- Craig PM, Trudeau VL, Moon TW. Profiling hepatic microRNAs in Zebrafish: fluoxetine exposure mimics a fasting response that targets AMP-Activated Protein Kinase (AMPK). PLoS One. 2014;9:e95351.View ArticlePubMedPubMed CentralGoogle Scholar
- Brzuzan P, Woêny M, Woliƒska L, Piasecka A, Łuczyƒski MK. MicroRNA expression in liver of whitefish (Coregonus lavaretus) exposed to microcystin-LR. Environ Biotechnol. 2010;6:53–60.Google Scholar
- Zha J, Wang Z, Wang N, Ingersoll C. Histological alternation and vitellogenin induction in adult rare minnow (Gobiocypris rarus) after exposure to ethynylestradiol and nonylphenol. Chemosphere. 2007;66:488–95.View ArticlePubMedGoogle Scholar
- Zhang X, Zha J, Li W, Yang L, Wang Z. Effects of 2,4-dichlorophenol on the expression of vitellogenin and estrogen receptor genes and physiology impairments in chinese rare minnow (Gobiocypris rarus). Environ Toxicol. 2008;23:694–701.View ArticlePubMedGoogle Scholar
- Jiang W, Yang Y, Zhao D, Liu X, Duan J, Xie S, Zhao H. Effects of sexual steroids on the expression of foxl2 in Gobiocypris rarus. Comp Biochem Physiol - B Biochem Mol Biol. 2011;160:187–93.Google Scholar
- Liang X, Wang M, Chen X, Zha J, Chen H, Zhu L, et al. Endocrine disrupting effects of benzotriazole in rare minnow (Gobiocypris rarus) in a sex-dependent manner. Chemosphere. 2014;112:154–62.Google Scholar
- Zhang X, Zha J, Wang Z. Influences of 4-nonylphenol on doublesex- and mab-3-related transcription factor 1 gene expression and vitellogenin mRNA induction of adult rare minnow (Gobiocypris rarus). Environ Toxicol Chem. 2008;27:196–205.View ArticlePubMedGoogle Scholar
- Yang L, Zha J, Li W, Li Z, Wang Z. Atrazine affects kidney and adrenal hormones (AHs) related genes expressions of rare minnow (Gobiocypris rarus). Aquat Toxicol. 2010;97:204–11.View ArticlePubMedGoogle Scholar
- Li W, Zha J, Spear P a, Li Z, Yang L, Wang Z. Changes of thyroid hormone levels and related gene expression in Chinese rare minnow (Gobiocypris rarus) during 3-amino-1,2,4-triazole exposure and recovery. Aquat Toxicol. 2009;92:50–7.View ArticlePubMedGoogle Scholar
- Liang X, Martyniuk CJ, Cheng G, Zha J, Wang Z. Pyruvate carboxylase as a sensitive protein biomarker for exogenous steroid chemicals. Environ Pollut. 2014;189:184–93.View ArticlePubMedGoogle Scholar
- Freeman JL, Weber GJ, Sepúlveda MS. Fishing for microRNAs in toxicology. In: microRNAs in toxicology and medicine. Chichester: Wiley; 2013. p. 49–75.View ArticleGoogle Scholar
- Chatterjee A, Rodger EJ, Stockwell PA, Weeks RJ, Morison IM. Technical considerations for reduced representation bisulfite sequencing with multiplexed libraries. J Biomed Biotechnol. 2012;2012:1–8.View ArticleGoogle Scholar
- Stockwell PA, Chatterjee A, Rodger EJ, Morison IM. DMAP: differential methylation analysis package for RRBS and WGBS data. Bioinformatics. 2014;30:1814–22.View ArticlePubMedGoogle Scholar
- Chatterjee A, Stockwell PA, Rodger EJ, Morison IM. Comparison of alignment software for genome-wide bisulphite sequence data. Nucleic Acids Res. 2012;40:e79.View ArticlePubMedPubMed CentralGoogle Scholar
- Betel D, Wilson M, Gabow A, Marks DS, Sander C. The microRNA.org resource: targets and expression. Nucleic Acids Res. 2007;36(Database):D149–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E. The role of site accessibility in microRNA target recognition. Nat Genet. 2007;39:1278–84.View ArticlePubMedGoogle Scholar
- Loher P, Rigoutsos I. Interactive exploration of RNA22 microRNA target predictions. Bioinformatics. 2012;28:3322–3.View ArticlePubMedGoogle Scholar
- Lewis BP, Burge CB, Bartel DP. Conserved Seed Pairing, Often Flanked by Adenosines, Indicates that Thousands of Human Genes are MicroRNA Targets. Cell. 2005;120:15–20.View ArticlePubMedGoogle Scholar
- Xu Z, Chen J, Li X, Ge J, Pan J, Xu X. Identification and characterization of MicroRNAs in Channel Catfish (Ictalurus punctatus) by using solexa sequencing technology. PLoS One. 2013;8:1–10.Google Scholar
- Wei C, Salichos L, Wittgrove CM, Rokas A, Patton JG. Transcriptome-wide analysis of small RNA expression in early zebrafish development. RNA. 2012;18:915–29.View ArticlePubMedPubMed CentralGoogle Scholar
- Xiao J, Zhong H, Zhou Y, Yu F, Gao Y, Luo Y, Tang Z, Guo Z, Guo E, Gan X, Zhang M, Zhang Y. Identification and characterization of microRNAs in ovary and testis of Nile tilapia (Oreochromis niloticus) by using solexa sequencing technology. PLoS One. 2014;9:e86821.Google Scholar
- Deng Y, Huang Z, Xu Y, Jin J, Zhuo W, Zhang C, Zhang X, Shen M, Yan X, Wang L, Wang X, Kang Y, Si J, Zhou T. MiR-215 modulates gastric cancer cell proliferation by targeting RB1. Cancer Lett. 2014;342:27–35.Google Scholar
- White NMA, Khella HWZ, Grigull J, Adzovic S, Youssef YM, Honey RJ, Stewart R, Pace KT, Bjarnason G a, Jewett M a S, Evans a J, Gabril M, Yousef GM. miRNA profiling in metastatic renal cell carcinoma reveals a tumour-suppressor effect for miR-215. Br J Cancer. 2011;105:1741–9.Google Scholar
- Ameres SL, Zamore PD. Diversifying microRNA sequence and function. Nat Rev Mol Cell Biol. 2013;14:475–88.View ArticlePubMedGoogle Scholar
- Morin RD, O’Connor MD, Griffith M, Kuchenbauer F, Delaney A, Prabhu A-L, Zhao Y, McDonald H, Zeng T, Hirst M, Eaves CJ, Marra M a. Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells. Genome Res. 2008;18:610–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Li SC, Liao YL, Chan WC, Ho MR, Tsai KW, Hu LY, Lai CH, Hsu CN, Lin WC. Interrogation of rabbit miRNAs and their isomiRs. Genomics. 2011;98:453–9.Google Scholar
- Janes RW, Peapus DH, Wallace B a. The crystal structure of human endothelin. Nat Struct Biol. 1994;1(April):311–9.View ArticlePubMedGoogle Scholar
- Ren G, Chen X, Yu B. Small RNAs meet their targets: when methylation defends miRNAs from uridylation. RNA Biol. 2015;11(January 2015):1099–104.Google Scholar
- Zhao Y, Yu Y, Zhai J, Ramachandran V, Dinh TT, Meyers BC, Mo B, Chen X. The arabidopsis nucleotidyl transferase HESO1 uridylates unmethylated small RNAs to trigger their degradation. Curr Biol. 2012;22:689–94.Google Scholar
- Ren G, Chen X, Yu B. Uridylation of miRNAs by HEN1 SUPPRESSOR1 in arabidopsis. Curr Biol. 2012;22:695–700.View ArticlePubMedPubMed CentralGoogle Scholar
- Xia J, Zhang W. A meta-analysis revealed insights into the sources, conservation and impact of microRNA 5′-isoforms in four model species. Nucleic Acids Res. 2014;42:1427–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Neilsen CT, Goodall GJ, Bracken CP. IsomiRs - the overlooked repertoire in the dynamic microRNAome. Trends Genet. 2012;28:544–9.View ArticlePubMedGoogle Scholar
- Siomi MC, Siomi H. Characterization of endogenous human Argonautes and their miRNA partners in RNA silencing. Nucleic Acids Symp Ser (Oxf). 2008;105:59–60.View ArticleGoogle Scholar
- Tan GC, Chan E, Molnar A, Sarkar R, Alexieva D, Isa IM, Robinson S, Zhang S, Ellis P, Langford CF, Guillot P V, Chandrashekran A, Fisk NM, Castellano L, Meister G, Winston RM, Cui W, Baulcombe D, Dibb NJ. 5′ isomiR variation is of functional and evolutionary importance. Nucleic Acids Res. 2014;42:1–12.Google Scholar
- Git A, Dvinge H, Salmon-Divon M, Osborne M, Kutter C, Hadfield J, Bertone P, Caldas C. Systematic comparison of microarray profiling, real-time PCR, and next-generation sequencing technologies for measuring differential microRNA expression. RNA. 2010;16:991–1006.Google Scholar
- Lynn FC. Meta-regulation: microRNA regulation of glucose and lipid metabolism. Trends Endocrinol Metab. 2009;20(September):452–9.View ArticlePubMedGoogle Scholar
- Esau C, Davis S, Murray SF, Yu XX, Pandey SK, Pear M, Watts L, Booten SL, Graham M, McKay R, Subramaniam A, Propp S, Lollo B a., Freier S, Bennett CF, Bhanot S, Monia BP. miR-122 regulation of lipid metabolism revealed by in vivo antisense targeting. Cell Metab. 2006;3(February):87–98.Google Scholar
- Mennigen JA. Micromanaging metabolism—a role for miRNAs in teleost energy metabolism. Comp Biochem Physiol Part B Biochem Mol Biol. 2015;1–11.Google Scholar
- Smith CA, Roeszler KN, Ohnesorg T, Cummins DM, Farlie PG, Doran TJ, Sinclair AH. The avian Z-linked gene DMRT1 is required for male sex determination in the chicken. Nature. 2009;461:267–71.Google Scholar
- Leucht C, Stigloher C, Wizenmann A, Klafke R, Folchert A, Bally-Cuif L. MicroRNA-9 directs late organizer activity of the midbrain-hindbrain boundary. Nat Neurosci. 2008;11:641–8.View ArticlePubMedGoogle Scholar
- Schulz RW, Bogerd J, Male R, Ball J, Fenske M, Olsen LC, Tyler CR. Estrogen-induced alterations in amh and dmrt1 expression signal for disruption in male sexual development in the Zebrafish estrogen-induced alterations in amh and dmrt1 expression signal for disruption in male sexual development in the Zebrafish. Environ Sci Technol. 2007;41:6305–10.Google Scholar
- Lindeman RE, Gearhart MD, Minkina A, Krentz AD, Bardwell VJ, Zarkower D. Sexual cell-fate reprogramming in the ovary by DMRT1. Curr Biol. 2015;25:764–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang DS, Kobayashi T, Zhou LY, Nagahama Y. Molecular cloning and gene expression of Foxl2 in the Nile tilapia, Oreochromis niloticus. Biochem Biophys Res Commun. 2004;320:83–9.View ArticlePubMedGoogle Scholar
- Nakamoto M, Matsuda M, Wang DS, Nagahama Y, Shibata N. Molecular cloning and analysis of gonadal expression of Foxl2 in the medaka, Oryzias latipes. Biochem Biophys Res Commun. 2006;344:353–61.View ArticlePubMedGoogle Scholar
- Baron D, Cocquet J, Xia X, Fellous M, Guiguen Y, Veitia R a. An evolutionary and functional analysis of FoxL2 in rainbow trout gonad differentiation. J Mol Endocrinol. 2004;33:705–15.View ArticlePubMedGoogle Scholar
- Yamaguchi T, Yamaguchi S, Hirai T, Kitano T. Follicle-stimulating hormone signaling and Foxl2 are involved in transcriptional regulation of aromatase gene during gonadal sex differentiation in Japanese flounder, Paralichthys olivaceus. Biochem Biophys Res Commun. 2007;359:935–40.View ArticlePubMedGoogle Scholar
- Heldring N, Pike A, Andersson S, Matthews J, Cheng G, Hartman J, Tujague M, Strom A, Treuter E, Warner M, Gustafsson J-A. Estrogen receptors: how do they signal and what are their targets. Physiol Rev. 2007;87:905–31.Google Scholar
- Weisz A, Rosales R. Identification of an estrogen response element upstream of the human c-fos gene that binds the estrogen receptor and the AP-1 transcription factor. Nucleic Acids Res. 1990;18:5097–106.View ArticlePubMedPubMed CentralGoogle Scholar
- Cicatiello L, Addeo R, Sasso A, Altucci L, Petrizzi VB, Borgo R, Cancemi M, Caporali S, Caristi S, Scafoglio C, Teti D, Bresciani F, Perillo B, Weisz A. Estrogens and progesterone promote persistent CCND1 gene activation during G1 by inducing transcriptional derepression via c-Jun/c-Fos/Estrogen receptor (Progesterone Receptor) complex assembly to a distal regulatory element and recruitment of Cyclin D1 t. Mol Cell Biol. 2004;24:7260–74.Google Scholar
- Crooks G, Hon G, Chandonia J, Brenner S. WebLogo: a sequence logo Generator Gavin. Genome Res. 2004;14:1188–90.View ArticlePubMedPubMed CentralGoogle Scholar