Parallel DNA pyrosequencing unveils new zebrafish microRNAs
© Soares et al; licensee BioMed Central Ltd. 2009
Received: 04 November 2008
Accepted: 27 April 2009
Published: 27 April 2009
MicroRNAs (miRNAs) are a new class of small RNAs of approximately 22 nucleotides in length that control eukaryotic gene expression by fine tuning mRNA translation. They regulate a wide variety of biological processes, namely developmental timing, cell differentiation, cell proliferation, immune response and infection. For this reason, their identification is essential to understand eukaryotic biology. Their small size, low abundance and high instability complicated early identification, however cloning/Sanger sequencing and new generation genome sequencing approaches overcame most technical hurdles and are being used for rapid miRNA identification in many eukaryotes.
We have applied 454 DNA pyrosequencing technology to miRNA discovery in zebrafish (Danio rerio). For this, a series of cDNA libraries were prepared from miRNAs isolated at different embryonic time points and from fully developed organs. Each cDNA library was tagged with specific sequences and was sequenced using the Roche FLX genome sequencer. This approach retrieved 90% of the 192 miRNAs previously identified by cloning/Sanger sequencing and bioinformatics. Twenty five novel miRNAs were predicted, 107 miRNA star sequences and also 41 candidate miRNA targets were identified. A miRNA expression profile built on the basis of pyrosequencing read numbers showed high expression of most miRNAs throughout zebrafish development and identified tissue specific miRNAs.
This study increases the number of zebrafish miRNAs from 192 to 217 and demonstrates that a single DNA mini-chip pyrosequencing run is effective in miRNA identification in zebrafish. This methodology also produced sufficient information to elucidate miRNA expression patterns during development and in differentiated organs. Moreover, some zebrafish miRNA star sequences were more abundant than their corresponding miRNAs, suggesting a functional role for the former in gene expression control in this vertebrate model organism.
MicroRNAs (miRNAs) are small RNAs that regulate eukaryotic gene expression at the post-transcriptional level . They are transcribed as long precursor RNA molecules (pri-miRNAs) and are successively processed by two key RNAses, namely Drosha and Dicer, into their mature forms of ~22 nucleotides . These small RNAs regulate gene expression by binding to target sites in the 3' untranslated region of mRNAs (3'UTR) [3, 4]. Recognition of the 3'UTR by miRNAs is mediated through complementary hybridization between nucleotides 2–8, numbered from the 5' end (seed sequences) of the small RNAs, and complementary sequences present in the 3'UTRs of mRNAs [3, 5, 6]. Perfect or nearly perfect complementarities between miRNAs and their 3'UTRs induce mRNA cleavage by the RNA-induced silencing complex (RISC), whereas imperfect base matching induces translational silencing through various molecular mechanisms , namely inhibition of translation initiation and activation of mRNA storage in P-bodies and/or stress granules . Interestingly, miRNAs also direct rapid deadenylation of target mRNAs, leading to decapping and rapid mRNA decay by the combined activities of the exosome (3' to 5' degradation) and the exoribonuclease Xrn1 (5' to 3'degradation) [7, 8].
Since the discovery of the first miRNA in 1993 in C. elegans , thousands of mature miRNAs have been uncovered in several species, suggesting that they appeared early in eukaryotic evolution and play fundamental roles in gene expression control. Indeed, miRNAs have been identified using a combination of bioinformatics, cloning and Sanger sequencing, and lately through new generation sequencing methods, namely the Roche 454 Pyrosequencer, the Solexa/Illumina Genome Analyzer and the Applied Biosystems SOLiD™ Sequencer, in a wide range of eukaryotes, namely plants [10–13], mammals [14–16], birds [17, 18], fish [19–21], amphibians , worms , flies , in the unicellular green algae Chlamydomonas reinhardtii  and in viruses . These small RNAs were originally thought to regulate developmental processes only [27–29], but recent studies show that they regulate a variety of other pivotal biological processes, namely differentiation , immune response , infection [32, 33] and cancer [34, 35]. The exact mechanism by which they regulate such a variety of molecular processes is not yet fully understood, however 2–3% of the human genes encode miRNAs and approximately 30% of human mRNAs contain miRNA binding sites in their 3'UTRs, suggesting major roles for these small RNAs in eukaryotic gene regulation .
The quantification of miRNA expression has been technically challenging and rather expensive due to their small size, low abundance, low stability and contamination with other cellular RNAs and mRNA fragments. Recently, the above mentioned parallel DNA sequencing methodologies have been successfully applied to both miRNA identification and quantification [11, 14, 16, 37]. The enormous sequencing power of these technologies has overcome most of the technical hurdles associated to miRNA identification and increased dramatically the number of miRNAs deposited in public databases . These new methodologies are also promoting large scale initiatives to identify most eukaryotic miRNAs, understand their evolution and identify target genes and gene networks regulated by them.
In zebrafish (ZF), 337 miRNA genes encode 192 different mature miRNAs . However, deep DNA sequencing has not yet been applied to this model organism and one is not sure whether those miRNAs represent the full ZF miRNA population [19, 20]. As in other eukaryotes, recent ZF studies highlighted critical miRNA roles in gene expression control since defective miRNA processing arrested development [39, 40]. Also, a specific subset of miRNAs is required for brain morphogenesis in ZF embryos, but not for cell fate determination or axis formation . In other words, miRNAs play an important role in ZF organogenesis and their expression at specific time points is relevant to organ formation and differentiation. Since identification of the complete set of miRNAs is fundamental to fully understand these and other fundamental biological processes, we have used high throughput 454 DNA pyrosequencing technologies to fully characterize the ZF miRNA population. This study increased the total number of ZF miRNAs from 192 to 217 and identified several star sequences (miRNA*, complementary to miRNA sequences). In addition, miRNAs predicted by homology were retrieved and novel miRNA genes encoding known miRNAs were identified.
454 DNA sequencing of zebrafish miRNAs
Distribution of microRNA reads in developmental and tissue samples.
Nr reads (perfect match to the genome)
Nr non-miRNA reads
Nr miRNA reads
Nr known miRNA reads
Nr putative new miRNA reads
Total number of miRNAs and miRNA genes identified in this study.
Number of known miRNAs detected
Number of new miRNAs identified in this study
Total number of miRNAs identified in this study
Number of known predicted miRNA genes detected
Number of novel genes identified in this study
Total of genes indentified in this study
miRNAs predicted by homology and experimentally validated.
Since 14,768 pyrosequencing reads did not match the ZF genome, we repeated the alignment using the mapping algorithm SHRiMP, which allows for introduction of insertions/deletions in the alignment. In this case, 87% (13,046 reads) of those reads produced matches in the ZF genome, however novel miRNAs were not identified. This analysis revealed some sequence variation in known miRNAs which may be related to sequencing errors or eventually to post-transcriptional miRNA editing . For example, 14 reads of the dre-miR-124 family had a (C->A) substitution at position 20, but the low number of reads of each sequence did not permit unequivocal differentiation between miRNA editing and sequencing errors. This should be investigated further in a new study.
miRNA expression patterns in zebrafish
The 24 hpf sample, corresponding to the developmental sample series, yielded low number of sequencing reads corresponding to 13 different miRNAs. The number of reads increased dramatically at late stages of development and 72 hpf sample produced the highest number of reads and the highest miRNA diversity (149 unique miRNAs) (see Additional File 3). At 96 hpf lower number of reads and lower miRNA diversity was observed. As ZF reached the adult stage, the number of reads and the number of different miRNAs raised again, which was consistent with previous studies [40, 51]. This suggested that miRNAs play an important role in differentiation and maintenance of tissue identity, rather than in tissue fate establishment . In adult fish, the brain sample produced the highest number of reads and the highest miRNA diversity (160 unique miRNAs). The gut/liver showed lower miRNA diversity (55 unique miRNAs) and the skin produced the lowest number of reads corresponding to 58 different miRNAs (see Additional File 3).
Novel zebrafish miRNAs are mostly non-conserved
A large set of both conserved and non-conserved miRNAs were previously identified in ZF by cloning/Sanger sequencing [19, 20], using miRNA cDNA libraries prepared mainly from brain and from few developmental stages, and also using bioinformatics . Our approach of isolating and preparing separated miRNA cDNA libraries from 24 hpf, 72 hpf, 96 hpf, 5 dpf, 45 dpf, total adult and from brain, eyes, heart, gills, muscle, fins, skin and gut/liver resulted in identification of 25 new miRNAs (see Additional File 4), using the miRDeep software package in combination with our data pipeline analysis (see Additional File 1). The miRDeep algorithm performed stringent searches based on the miRNA biogenesis model  and produced information on miRNA conservation, thermodynamic stability, ability to form a hairpin with the shape and sequence of the pre-miRNA molecule, number of sequences that matched mature miRNA sequences and number of sequences that matched star sequences. This algorithm alone was able to detect 153 known and 23 novel miRNAs. Our data pipeline (see Additional File 1) identified 20 known miRNAs missed by miRDeep in our samples, resulting in a total of 173 known miRNAs identified. The extension of this analysis to sequences predicted by Ensembl and by Thatcher et al  permitted detection of two additional novel miRNAs, raising the prediction number to 25.
Novel zebrafish miRNAs. miRNAs identified in this study, their level of conservation and corresponding miRNA families.
The expression of a conserved miRNA (miR_4), a non-conserved miRNA for which the star sequence was also detected (miR_8) and a non-conserved miRNA for which the star sequence was not detected (miR_9), was validated by qRT-PCR (Figure 5B). As above, there was strong correlation between qRT-PCR and 454 pyrosequencing data. Indeed, miR_8 displayed the highest number of pyrosequencing reads and had higher relative abundance in the qRT-PCR analysis. Also, miR_4 and miR_9 had similar levels of expression, but lower than miR_8.
Target prediction for novel miRNAs
The miRNA targets can be predicted computationally with high confidence for conserved miRNAs, but such predictions remain challenging for non-conserved miRNAs due to restrictions imposed by the search algorithms used in the target prediction databases . For the non-conserved miRNAs only the more extensively paired interactions can be predicted with reasonable confidence. In order to minimize noise (false predictions) in the prediction of targets, stringent criteria similar to those described by Sunkar and colleagues  were used. This was based on blast searches for antisense hits with less than 6 mismatches, with perfect seed match and thermodynamic stability using the RNA HYBRID software. Forty two putative gene targets of the 16 newly identified miRNAs, which were mainly involved in binding nucleotides, proteins or ions or had catalytic activity, were identified (see Additional File 5). The predominant biological functions of these predicted target genes were cellular processes related to metabolism and signal transduction and developmental processes, including embryonic patterning, vasculogenesis and neuron differentiation. Most miRNAs with predicted targets involved in developmental processes were detected in cDNA libraries prepared from miRNA samples collected during embryonic development. For example, miR_7 was detected at 5 dpf; miR_8 was detected at 72 hpf, 96 hpf and 5 dpf; miR_9 and miR_10 were both detected at 72 hpf and 96 hpf; miR_16 was detected at 72 hpf and 5 dpf. The MELK gene which is found in gills and is involved in erythrocyte development was predicted to be targeted by miRNA_14 raising the hypothesis that miR_14 could be involved in erythrocyte development . Other possible targets encoded hypothetical proteins and, for this reason, were not included in our analysis. This approach was unable to identify candidate targets for some of the novel miRNAs, a result that may be explained by the high stringency of the prediction algorithm because targets of conserved miR_4 and miR_11 were also missed. Alternatively, some of these low abundance non-conserved miRNAs may have appeared recently and do not yet have targets, or misannotation or incomplete annotation of the ZF genome may have prevented identification of such targets.
Star sequences of zebrafish microRNAs
To date, 192 ZF miRNAs have been identified using classical cloning and Sanger sequencing methodologies [19, 20]. In this study, 25 novel miRNAs were added to the ZF repertoire using massively parallel DNA pyrosequencing of miRNA cDNA libraries prepared from different time points of the ZF embryo development and from different tissues. This methodology retrieved 173 of the 192 known ZF miRNAs whose expression in different tissues and developmental stages were validated using Northern blot analysis and/or in situ hybridizations. This high degree of data overlapping between cloning/Sanger sequencing and DNA pyrosequencing, plus the existence of target genes for the novel miRNAs, validated our approach of miRNA libraries fractionation and provided strong support for the authenticity of the newly identified miRNAs. The pyrosequencing approach also produced important information about the relative abundance of the ZF miRNAs. During early development (24 hpf), the number of different miRNAs was low (13), but expression level was high. The number of miRNA reads was higher at 72 hpf and in young adult fish. Among differentiated organs, brain and eyes showed the highest number of miRNA reads. This confirmed previous data showing differences in temporal miRNA expression and raised the hypothesis that many miRNAs play a role in late development and are required for organ morphogenesis .
Zebrafish microRNAs expression profile
The pyrosequencing data allowed us to build a miRNA expression profile for developmental differentiation and for adult fish, based on the normalized number of reads. Most miRNAs were expressed in more than one tissue (Figure 4A), others were tissue specific or showed stronger expression in specific tissues (Figure 4B), while others were development specific. Dre-miR-135c and dre-miR-25 were highly enriched at 24 hpf, but their relative expression decreased during embryo development. The data confirmed previous studies showing that dre-miR-135 expression is higher in development than in adult fish . The miR-430 family was also present during development and was absent in adult ZF [19, 20]. The expression profile also highlighted results of Giraldez and colleagues  showing that miR-430 is essential for regulation of morphogenesis during development.
Some miRNAs were expressed ubiquitously. For example, dre-miR-124 was abundant during both development and in adult fish, and its expression increased slightly during late stages of development and in the central nervous system (both brain and eyes). This miRNA alone accounted for ~48% of the pyrosequencing reads, a result that may be explained, at least in part, by the high copy number of its gene (6 copies in various chromosomes). At 24 hpf, when a significant part of the brain development had already occurred, dre-miR-124 represented 42% of the miRNA pool and its relative abundance reached 80% at 5 dpf. In the adult tissues, it represented 80% of brain and 54% of eye miRNAs. This is in agreement with previous studies in ZF and other organisms showing that miR-124 is up-regulated during development of the nervous system and is the most abundant miRNA in the adult brain [19, 20]. Also, neuronal differentiation is enhanced followed transfection of mir-124 in mouse neuronal stem cells . Taken together, the data suggest that dre-miR-124 may play an important role in ZF development, neuronal differentiation and in regulation of brain functions [35, 58]. On the other hand, dre-miR-203a and dre-miR-203b appeared early in development and maintained high levels of expression in adult fish, in particular in gills and skin. Indeed, miR-203 is a skin-specific keratinocyte-derived miRNA involved in keratinocyte differentiation .
A subset of miRNAs was expressed in differentiated tissues only. For example, dre-miR-101a, dre-miR-130b, dre-miR-130c, dre-miR-221 and dre-miR-499 were highly enriched in the heart, in agreement with previous in situ and Northern blot studies . Dre-miR-1 and dre-miR-133a were expressed in muscle and heart, where they play an important regulatory role in other organisms [60, 61]. Indeed, deletion of miR-1 altered regulation of cardiogenesis, electrical conduction and cell cycle of cardiomyocites, and miR-133 plus miR-1 regulate cardiac hypertrophy, as their over expression inhibits it. Interestingly, dre-miR-133b and dre-miR-133c were mainly detected in muscle and were not present in the heart. Finally, dre-miR-103 was specific of gut and liver while dre-miR-122 was liver specific [40, 62]. This was not surprising because mir-122 plays important roles in regulation of metabolism and its silencing in high-fat fed mice resulted in a significant reduction of hepatic steatosis, decreased cholesterol synthesis and stimulated fatty-acid oxidation .
Expression and putative functions of the novel zebrafish miRNAs
Of the 25 novel miRNAs identified in this study, 9 belong to conserved miRNA families (existing in at least one more organism) according to the conservation criteria used in this study, and the others are non-conserved (ZF specific). Most of the novel miRNAs are encoded by a single gene, but 7 are multigenic. In the latter case, miRDeep, Ensembl and RNAfold analysis showed that different genes encoding a single miRNA produce identical miRNA hairpins. Most of the novel miRNAs produced lower number of reads than the majority of the conserved miRNAs. This was not surprising since there is good correlation between miRNA conservation and expression level . Therefore, the low abundance of the novel miRNAs identified by pyrosequencing combined with the retrieval of 90% of the known ZF miRNAs (identified by cloning/Sanger sequencing) suggests that the vast majority of miRNAs present in our samples were retrieved. However, one cannot exclude the hypothesis that new miRNAs present in our dataset escaped identification due to the high stringency of the methodology used. Also, it is possible that other low abundance and highly specific ZF miRNAs may still be discovered using other deep DNA pyrosequencing strategies, namely Solexa/Illumina or SOLiD™. Finally, cDNA libraries from tissues not screened in this study may still reveal new ZF miRNAs. Recent bioinformatics analysis of the ZF genome identified additional miRNAs , however we were unable to identify reads matching these putative miRNAs using miRDeep alone or our pipeline data analysis system. This may be due to their very low expression level. Again, other massively parallel DNA pyrosequencing approaches may overcome these limitations and uncover such putative miRNAs .
Our bioinformatics approach retrieved 41 candidate target genes of 15 novel miRNAs. Since we used stringent search criteria to minimize false positives one cannot exclude the possibility that some targets were missed. Despite this, comparative analysis of the targets of the conserved miR_15, miR_16 and miR_21 with those of known miRNAs produced significant overlapping, thus validating our target prediction approach. For example, miR_16, which belongs to the miR-107 family, has GFM2 and VOX genes as putative targets. The miRBase Targets Version 5 also retrieved these genes as candidate targets for dre-miR-107. Similar results were obtained for miR_21 where RNF11 gene was highlighted as a candidate target of this novel miRNA of the mir-222 family. This result is also supported by retrieval of miR-222 in a blast search for miRNAs that target RNF11.
Most of the predicted targets are involved in cellular and developmental processes, which is in agreement with their expression during development. Indeed, the NRP2A gene, a putative target of miR_8, is involved in the differentiation of the nervous system, neural crest cell migration  and in VEGF-mediated vessel development . This correlates with the expression pattern of this miRNA at 72 hpf, 96 hpf, 5 dpf and in the adult brain. Also, miR_9, expressed at 72 hpf and 96 hpf, is predicted to target the PRDM1 and zgc:85707 genes which play important roles in embryonic axis specification, embryonic pectoral fin morphogenesis, regulation of neuron specification, regulation of signal transduction and multicellular organism development . The SEC23B and MYST3 genes which are involved in cartilage development were retrieved as putative targets of miR_10, which was also expressed during development. MYST3 (or MOX) regulates HOX expression and segmental identity, including cartilage patterning . Finally, the gills specific miR_14 is predicted to target the MELK gene, which is also strongly expressed in the gills and is involved in erythrocyte development .
Zebrafish microRNA star sequences
Star sequences of miRNAs (miRNA*) are difficult to detect by conventional methods due to their rapid turnover. However, high throughput sequencing retrieved many of them and revealed their relative abundance in different organisms [19, 20, 57, 64]. Our DNA pyrosequencing approach identified 107 miRNA* sequences: 42 were identified previously by cloning and Sanger sequencing [19, 20], 60 were identified in this pyrosequencing study, but belong to known miRNAs, and other 5 miRNA* belong to the novel miRNAs identified in this study. Most star sequences retrieved fewer reads than the corresponding mature miRNAs which is consistent with the miRNA biogenesis model and strand selection by RISC. However, six miRNA* were more abundant than their corresponding mature miRNAs, namely dre-miR-129*, dre-miR-140*, dre-miR-142a*, dre-miR-202*, dre-miR-210* and dre-miR-214*. Similar results were observed before for dre-mir-129*, dre-mir-142a*, dre-mir-142b* and dre-mir-214* . Dre-miR-30e, dre-miR-199, dre-miR-219 and dre-miR-462 showed similar strand-bias of both mature and star strands. Since this was also observed in the chicken embryo for mir-30e and mir-219  and in Drosophila melanogaster where several miRNA* are present at physiologically relevant levels and associate with Argonaute proteins, it is likely that both strands are loaded into RISC and may guide target repression. Finally, observed alterations in the ratio of expression of mature/star molecules suggests that some star molecules are functional and their activity may vary according to cellular context [57, 64]. Obviously, the biological function of these star sequences can only be unravelled through experimental testing, but their high number of reads suggests their inclusion in future ZF miRNA chips and expression profiling studies.
This study increased the total number of ZF miRNAs from 192 to 217 and showed that miRNA cDNA libraries prepared from different developmental stages and from adult tissues is an effective methodology for miRNA discovery using low cost DNA pyrosequencing mini-chips. The methodology permitted quantitative and qualitative analysis of miRNA expression throughout the ZF life cycle, as the miRNA profile was largely in agreement with qRT-PCR and Sanger sequencing data. Most of the 25 novel miRNAs were non-conserved low abundance molecules and their targets indicated that they might be involved in developmental processes. Novel miRNA star sequences were also identified and some of them were more abundant than their corresponding mature miRNAs, suggesting that they may also be loaded into RISC and may be functional. Future deep sequencing studies may still identify additional miRNAs in ZF, however such miRNAs may be expressed at very low levels or in specific physiological or pathological conditions.
MicroRNA library construction and sequencing
Small RNA libraries were prepared from different ZF developmental stages, namely 24 hours post-fertilization (hpf), 72 hpf, 96 hpf, 5 days post-fertilization (dpf), 45 dpf, young adult and from adult brain, eyes, gills, muscle, heart, skin, fins, and gut/liver (Figure 1). Briefly, 100 μg of total RNA from each sample was isolated using TRIzol® and small RNAs were enriched by differential precipitation using polyethylene glycol. Total RNAs were fractioned using 12% denaturing PAGE and small RNAs of 15–30 nt were gel-isolated using Gel Filtration cartridges from Edge Biosystems. For cDNA synthesis, the small RNA molecules previously isolated were first ligated to a 3' adapter (AMP-5'p-5'p/CTGTAGGCACCATCAATdi-deoxyC- 3') in absence of ATP and gel excised in the range of 35 and 50 nt. A second ligation was performed with the 5' adapter ("Nelson's linker" 5'ATCGTrArGrGrCrArCrCrUrGrArArA 3'), for 1 hour at 37°C, followed by phenol extraction. First strand cDNA synthesis was then performed using a specific 3'-primer and Superscript™ III reverse transcriptase (Invitrogen). RNAseH treated cDNA was PCR-amplified with adapter specific primers. PCR products were then run on 10% denaturing PAGE containing 7 M urea and the corresponding band (100 nt) was eluted from the gel with Probe Elution Buffer from Ambion, at 37°C overnight. Parallel DNA pyrosequencing was performed using the Genome Sequencer FLX (Roche), following established protocols for DNA library sequencing .
Computational analysis of sequencing reads
Base calling and quality trimming of sequence reads was carried out using the Genome Sequencer FLX software. Raw images were processed to remove background noise and the data was normalized. TAGs and adapter sequences of ZF developmental and adult tissues samples were then identified and trimmed and those reads with correct TAGs and adapters (> 15 nt) were retrieved for downstream analysis using the miRDeep software http://www.mdc-berlin.de/en/research/research_teams/systems_biology_of_gene_regulatory_elements/projects/miRDeep/index.html, with a cut-off value of 1. Initially, miRDeep aligned the sequences against the zebrafish genome using megaBlast with seed length set at 12, the traditional blast output, and minimum local identity set at 100. The blast output was then parsed for miRDeep uploading and aligned sequences with a maximum of 2 mismatches in the 3' end were retrieved. Reads that matched more than 10 different genome loci were discarded and only those with one or more alignments were kept and, using the remaining alignments as guidelines, the potential precursors were excised from the genome. The secondary structure of putative precursors was predicted using RNAfold and signatures were created by retaining reads that aligned perfectly with those putative precursors to generate the signature format. Finally, miRDeep predicted miRNAs by discarding non-plausible Dicer products and scoring plausible ones. To assess seed conservation, plausible Dicer processing sequences were blasted against a local version of mature miRNAs from miRBase 12.0 that lacked zebrafish miRNA sequences. Borderline miRNA candidates were also resolved by determining their relative stability using Randfold. To distinguish between novel and known miRNAs, selected pre-miRNAs were blasted against Danio rerio stem loop sequences (miRBase) and those that did not produce any or produced perfect alignments were scored as novel miRNAs. Pairs of signatures and structures were used to estimate the number of false positives by randomly permutating them, using miRDeep.
To estimate the false negative rate, known mature miRNAs deposited in miRBase 12.0 were used to carry out a megaBlast search using our sequencing data set. Perfect alignments were considered as true positives. This list of miRNAs was then compared with that of the miRNAs predicted by miRDeep and the sequences present in the blast list but absent in the miRDeep list were considered false negatives. To overcome the inherent lack of sensitivity of miRDeep, novel transcripts encoding miRNAs predicted by bioinformatics were retrieved from Ensembl 5.2 using BioMart and from literature predictions . These sequences were then used to perform a megaBlast search against our data with seed length set at 12. The transcripts with perfect matches and alignment length larger than 18 nucleotides were kept for further processing. These transcripts were then compared with the mature miRNAs present in miRBase 12.0 and those that produced imperfect alignments or did not produce alignments were considered new miRNAs.
Reads without matches in the ZF genome (megaBlast) were re-aligned using SHRiMP, http://compbio.cs.toronto.edu/shrimp/ which handles short reads. Alignments were carried out using the space seed 011111111000; where 1 is the number of seed matches per window. Suboptimal alignments were retrieved and transformed into the blastparsed format for miRDeep miRNA prediction.
Statistical analysis of miRNA population
A rarefaction analysis of the detected miRNA population was carried out to assess the representativeness of the miRNA reads. A rarefaction curve, of the total number of reads obtained vs the total number of miRNA species, was plotted. The non-parametric richness estimator, Chao1  was determined to predict the total richness of the miRNA population, as a function of the observed richness (S obs ), the number of sequences observed only once (singletons, n 1 ) and the number of sequences observed twice (doubletons, n 2 ). A file containing the total number of reads of each miRNA was generated and used as input data set for the EstimateS8.0 statistical package . Both the rarefaction curve and the Chao1 statistical estimator were computed using EstimateS8.0.
Zebrafish miRNA expression profile
Where NRmiRNAXY is the number of reads of miRNAX (X = any miRNA) in sample Y, and TNRmiRNAsY is the total number of miRNAs in sample Y. 1000 is an arbitrary number of reads. The data was transformed into log2 scale to build the heat map using the MeV 4.0 software package http://www.tm4.org/mev.html.
MicroRNA expression analysis by quantitative Real-Time PCR
where Ct(sc) is the Ct expected for a sample containing a single copy template and E is the PCR efficiency. All Ct values above 35 were set to 35 and a PCR efficiency of 0.90 and Ct(sc) = 35 were assumed .
The 3'UTR sequences of ZF mRNAs were extracted from Biomart http://www.biomart.org/ and blasted against the antisense miRNA sequence for the new miRNAs or against the antisense miRNA* sequence, in the case of the star sequences. Sequences with perfect seed match between nucleotides 2 and 7, and no more than 6 mismatches in the remaining sequence, were retained for further analysis. Targets were considered positive whenever RNAhybrid confirmed them thermodynamically. Targets were discarded when RNAhybrid did not retrieve the targets obtained in the first approach.
- 3' UTR:
3' Untranslated Region
RNA induced silencing complex
hours post fertilization
days post fertilization
We are most thankful to Dr. Nicolaus Rajewsky and Marc Friedlaender from the Max Delbrück Centrum für Molekulare Medizin, Berlim, for making miRDeep available before publication and for helpful discussions of the data. ARS and PMP are supported by the Portuguese Foundation for Science and Technology (FCT). We are also thankful to the Centre for Neurosciences of the University of Coimbra for providing a PhD fellowship to ARS and to Biocant for making available the Genome Sequencer FLX. This project was funded by FCT/FEDER project SAU-MMO/55476/2004.
- Filipowicz W, Bhattacharyya SN, Sonenberg N: Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?. Nat Rev Genet. 2008, 9: 102-114. 10.1038/nrg2290.View ArticlePubMedGoogle Scholar
- Pillai RS: MicroRNA function: multiple mechanisms for a tiny RNA?. RNA. 2005, 11: 1753-1761. 10.1261/rna.2248605.PubMed CentralView ArticlePubMedGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.View ArticlePubMedGoogle Scholar
- Nilsen TW: Mechanisms of microRNA-mediated gene regulation in animal cells. Trends Genet. 2007, 23: 243-249. 10.1016/j.tig.2007.02.011.View ArticlePubMedGoogle Scholar
- Ambros V: The functions of animal microRNAs. Nature. 2004, 431: 350-355. 10.1038/nature02871.View ArticlePubMedGoogle Scholar
- Zamore PD, Haley B: Ribo-gnome: the big world of small RNAs. Science. 2005, 309: 1519-1524. 10.1126/science.1111444.View ArticlePubMedGoogle Scholar
- He L, Hannon GJ: MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004, 5: 522-531. 10.1038/nrg1379.View ArticlePubMedGoogle Scholar
- Wu L, Fan J, Belasco JG: MicroRNAs direct rapid deadenylation of mRNA. Proc Natl Acad Sci USA. 2006, 103: 4034-4039. 10.1073/pnas.0510928103.PubMed CentralView 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-854. 10.1016/0092-8674(93)90529-Y.View ArticlePubMedGoogle Scholar
- Reinhart BJ, Weinstein EG, Rhoades MW, Bartel B, Bartel DP: MicroRNAs in plants. Genes & Development. 2002, 16: 1616-1626. 10.1101/gad.1004402.View ArticleGoogle Scholar
- Sunkar R, Zhou XF, Zheng Y, Zhang WX, Zhu JK: Identification of novel and candidate miRNAs in rice by high throughput sequencing. Bmc Plant Biology. 2008, 8: 25-10.1186/1471-2229-8-25.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang JF, Zhou H, Chen YQ, Luo QJ, Qu LH: Identification of 20 microRNAs from Oryza sativa. Nucleic Acids Research. 2004, 32: 1688-1695. 10.1093/nar/gkh332.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang XJ, Reyes JL, Chua NH, Gaasterland T: Prediction and identification of Arabidopsis thaliana microRNAs and their mRNA targets. Genome Biol. 2004, 5 (9): R65-10.1186/gb-2004-5-9-r65.PubMed CentralView ArticlePubMedGoogle Scholar
- Berezikov E, Thuemmler F, van Laake LW, Kondova I, Bontrop R, Cuppen E, et al: Diversity of microRNAs in human and chimpanzee brain. Nature Genetics. 2006, 38: 1375-1377. 10.1038/ng1914.View ArticlePubMedGoogle Scholar
- Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, et al: A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 2007, 129: 1401-1414. 10.1016/j.cell.2007.04.040.PubMed CentralView ArticlePubMedGoogle Scholar
- Morin RD, O'Connor MD, Griffith M, Kuchenbauer F, Delaney A, Prabhu AL, et al: Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells. Genome Research. 2008, 18: 610-621. 10.1101/gr.7179508.PubMed CentralView ArticlePubMedGoogle Scholar
- Burnside J, Ouyang M, Anderson A, Bernberg E, Lu C, Meyers BC, et al: Deep sequencing of chicken microRNAs. Bmc Genomics. 2008, 9: 185-10.1186/1471-2164-9-185.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu HT, Wang XB, Du ZL, Li N: Identification of microRNAs from different tissues of chicken embryo and adult chicken. Febs Letters. 2006, 580: 3610-3616. 10.1016/j.febslet.2006.05.044.View ArticlePubMedGoogle Scholar
- Chen PY, Manninga H, Slanchev K, Chien MC, Russo JJ, Ju JY, et al: The developmental miRNA profiles of zebrafish as determined by small RNA cloning. Genes & Development. 2005, 19: 1288-1293. 10.1101/gad.1310605.View ArticleGoogle Scholar
- Kloosterman WP, Steiner FA, Berezikov E, de Bruijn E, Belt van de J, Verheul M, et al: Cloning and expression of new microRNAs from zebrafish. Nucleic Acids Research. 2006, 34: 2558-2569. 10.1093/nar/gkl278.PubMed CentralView 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 Developmental Biology. 2008, 8: 41-10.1186/1471-213X-8-41.PubMed CentralView ArticlePubMedGoogle Scholar
- Watanabe T, Takeda A, Mise K, Okuno T, Suzuki T, Minami N, et al: Stage-specific expression of microRNAs during Xenopus development. Febs Letters. 2005, 579: 318-324. 10.1016/j.febslet.2004.11.067.View ArticlePubMedGoogle Scholar
- Lim LP, Lau NC, Weinstein EG, Abdelhakim A, Yekta S, Rhoades MW, et al: The microRNAs of Caenorhabditis elegans. Genes & Development. 2003, 17: 991-1008. 10.1101/gad.1074403.View ArticleGoogle Scholar
- Lu J, Shen Y, Wu Q, Kumar S, He B, Shi S, et al: The birth and death of microRNA genes in Drosophila. Nat Genet. 2008, 40: 351-355. 10.1038/ng.73.View ArticlePubMedGoogle Scholar
- Zhao T, Li GL, Mi SJ, Li S, Hannon GJ, Wang XJ, et al: A complex system of small RNAs in the unicellular green alga Chlamydomonas reinhardtii. Genes & Development. 2007, 21: 1190-1203. 10.1101/gad.1543507.View ArticleGoogle Scholar
- Pfeffer S, Zavolan M, Grasser FA, Chien MC, Russo JJ, Ju JY, et al: Identification of virus-encoded microRNAs. Science. 2004, 304: 734-736. 10.1126/science.1096781.View ArticlePubMedGoogle Scholar
- Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T: Identification of novel genes coding for small expressed RNAs. Science. 2001, 294: 853-858. 10.1126/science.1064921.View ArticlePubMedGoogle Scholar
- Lau NC, Lim LP, Weinstein EG, Bartel DP: An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditis elegans. Science. 2001, 294: 858-862. 10.1126/science.1065062.View ArticlePubMedGoogle Scholar
- Lee RC, Ambros V: An extensive class of small RNAs in Caenorhabditis elegans. Science. 2001, 294: 862-864. 10.1126/science.1065329.View ArticlePubMedGoogle Scholar
- Tay YMS, Tam WL, Ang YS, Gaughwin PM, Yang H, Wang WJ, et al: MicroRNA-134 modulates the differentiation of mouse embryonic stem cells, where it causes post-transcriptional attenuation of Nanog and LRH1. Stem Cells. 2008, 26: 17-29. 10.1634/stemcells.2007-0295.View ArticlePubMedGoogle Scholar
- Tili E, Michaille JJ, Cimino A, Costinean S, Dumitru CD, Adair B, et al: Modulation of miR-155 and miR-125b levels following lipopolysaccharide/TNF-alpha stimulation and their possible roles in regulating the response to endotoxin shock. J Immunol. 2007, 179 (8): 5082-5089.View ArticlePubMedGoogle Scholar
- Gupta A, Gartner JJ, Sethupathy P, Hatzigeorgiou AG, Fraser NW: Anti-apoptotic function of a microRNA encoded by the HSV-1 latency-associated transcript. Nature. 2006, 442: 82-85.PubMedGoogle Scholar
- Jopling CL, Yi MK, Lancaster AM, Lemon SM, Sarnow P: Modulation of hepatitis C virus RNA abundance by a liver-specific microRNA. Science. 2005, 309: 1577-1581. 10.1126/science.1113329.View ArticlePubMedGoogle Scholar
- Huang QH, Gumireddy K, Schrier M, Le Sage C, Nagel R, Nair S, et al: The microRNAs miR-373 and miR-520c promote tumour invasion and metastasis. Nature Cell Biology. 2008, 10: 202-U83. 10.1038/ncb1681.View ArticlePubMedGoogle Scholar
- Silber J, Lim DA, Petritsch C, Persson AI, Maunakea AK, Yu M, et al: miR-124 and miR-137 inhibit proliferation of glioblastoma multiforme cells and induce differentiation of brain tumor stem cells. BMC Med. 2008, 6: 14-10.1186/1741-7015-6-14.PubMed CentralView ArticlePubMedGoogle Scholar
- varez-Garcia I, Miska EA: MicroRNA functions in animal development and human disease. Development. 2005, 132: 4653-4662. 10.1242/dev.02073.View ArticleGoogle Scholar
- Hafner M, Landgraf P, Ludwig J, Rice A, Ojo T, Lin C, et al: Identification of microRNAs and other small regulatory RNAs using cDNA library sequencing. Methods. 2008, 44: 3-12. 10.1016/j.ymeth.2007.09.009.PubMed CentralView ArticlePubMedGoogle Scholar
- Griffiths-Jones S: The microRNA Registry. Nucleic Acids Research. 2004, 32: D109-D111. 10.1093/nar/gkh023.PubMed CentralView ArticlePubMedGoogle Scholar
- Giraldez AJ, Mishima Y, Rihel J, Grocock RJ, van DS, Inoue K, et al: Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science. 2006, 312: 75-79. 10.1126/science.1122689.View ArticlePubMedGoogle Scholar
- Wienholds E, Kloosterman WP, Miska E, varez-Saavedra E, Berezikov E, de BE, et al: MicroRNA expression in zebrafish embryonic development. Science. 2005, 309: 310-311. 10.1126/science.1114519.View ArticlePubMedGoogle Scholar
- Giraldez AJ, Cinalli RM, Glasner ME, Enright AJ, Thomson JM, Baskerville S, et al: MicroRNAs regulate brain morphogenesis in zebrafish. Science. 2005, 308: 833-838. 10.1126/science.1109020.View ArticlePubMedGoogle Scholar
- Droege M, Hill B: The Genome Sequencer FLX trade mark System-Longer reads, more applications, straight forward bioinformatics and more complete data sets. J Biotechnol. 2008, 136 (1–2): 3-10. 10.1016/j.jbiotec.2008.03.021.View ArticlePubMedGoogle Scholar
- Friedlander MR, Chen W, Adamidi C, Maaskola J, Einspanier R, Knespel S, et al: Discovering microRNAs from deep sequencing data using miRDeep. Nature Biotechnology. 2008, 26: 407-415. 10.1038/nbt1394.View ArticlePubMedGoogle Scholar
- Thatcher EJ, Bond J, Paydar I, Patton JG: Genomic organization of zebrafish microRNAs. BMC Genomics. 2008, 9: 253-10.1186/1471-2164-9-253.PubMed CentralView ArticlePubMedGoogle Scholar
- Subramanian S, Fu Y, Sunkar R, Barbazuk WB, Zhu JK, Yu O: Novel and nodulation-regulated microRNAs in soybean roots. BMC Genomics. 2008, 9:Google Scholar
- Schloss PD, Handelsman J: The last word: Books as a statistical metaphor for microbial communities. Annual Review of Microbiology. 2007, 61: 23-34. 10.1146/annurev.micro.61.011507.151712.View ArticlePubMedGoogle Scholar
- Sogin ML, Morrison HG, Huber JA, Mark Welch D, Huse SM, Neal PR, et al: Microbial diversity in the deep sea and the underexplored "rare biosphere". Proceedings of the National Academy of Sciences of the United States of America. 2006, 103: 12115-12120. 10.1073/pnas.0605127103.PubMed CentralView ArticlePubMedGoogle Scholar
- Colwell RK, Coddington JA: Estimating Terrestrial Biodiversity Through Extrapolation. Philosophical Transactions of the Royal Society of London Series B-Biological Sciences. 1994, 345: 101-118. 10.1098/rstb.1994.0091.View ArticleGoogle Scholar
- Chao A: Estimating the Population-Size for Capture Recapture Data with Unequal Catchability. Biometrics. 1987, 43: 783-791. 10.2307/2531532.View ArticlePubMedGoogle Scholar
- Reid JG, Nagaraja AK, Lynn FC, Drabek RB, Muzny DM, Shaw CA, et al: Mouse let-7 miRNA populations exhibit RNA editing that is constrained in the 5'-seed/cleavage/anchor regions and stabilize predicted mmu-let-7a: mRNA duplexes. Genome Research. 2008, 18: 1571-1581. 10.1101/gr.078246.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Thatcher EJ, Flynt AS, Li N, Patton JR, Patton JG: miRNA expression analysis during normal zebrafish development and following inhibition of the hedgehog and notch signaling pathways. Developmental Dynamics. 2007, 236: 2172-2180. 10.1002/dvdy.21211.View ArticlePubMedGoogle Scholar
- Artzi S, Kiezun A, Shomron N: miRNAminer: a tool for homologous microRNA gene search. BMC Bioinformatics. 2008, 9: 39-10.1186/1471-2105-9-39.PubMed CentralView 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. 10.1016/j.cell.2004.12.035.View ArticlePubMedGoogle Scholar
- Cunha SR, Mohler PJ: Cardiac ankyrins: Essential components for development and maintenance of excitable membrane domains in heart. Cardiovascular Research. 2006, 71: 22-29. 10.1016/j.cardiores.2006.03.018.View ArticlePubMedGoogle Scholar
- Maziere P, Enright AJ: Prediction of microRNA targets. Drug Discovery Today. 2007, 12: 452-458. 10.1016/j.drudis.2007.04.002.View ArticlePubMedGoogle Scholar
- Saito R, Tabata Y, Muto A, Arai K, Watanabe S: Melk-like kinase plays a role in hematopoiesis in the zebra fish. Molecular and Cellular Biology. 2005, 25: 6682-6693. 10.1128/MCB.25.15.6682-6693.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Okamura K, Phillips MD, Tyler DM, Duan H, Chou YT, Lai EC: The regulatory activity of microRNA* species has substantial influence on microRNA and 3' UTR evolution. Nat Struct Mol Biol. 2008, 15: 354-363. 10.1038/nsmb.1409.PubMed CentralView ArticlePubMedGoogle Scholar
- Makeyev EV, Zhang JW, Carrasco MA, Maniatis T: The MicroRNA miR-124 promotes neuronal differentiation by triggering brain-specific alternative Pre-mRNA splicing. Molecular Cell. 2007, 27: 435-448. 10.1016/j.molcel.2007.07.015.PubMed CentralView ArticlePubMedGoogle Scholar
- Sonkoly E, Wei T, Janson PC, Saaf A, Lundeberg L, Tengvall-Linder M, et al: MicroRNAs: novel regulators involved in the pathogenesis of Psoriasis?. PLoS ONE. 2007, 2: e610-10.1371/journal.pone.0000610.PubMed CentralView ArticlePubMedGoogle Scholar
- Care A, Catalucci D, Felicetti F, Bonci D, Addario A, Gallo P, et al: MicroRNA-133 controls cardiac hypertrophy. Nat Med. 2007, 13: 613-618. 10.1038/nm1582.View ArticlePubMedGoogle Scholar
- Zhao Y, Ransom JF, Li A, Vedantham V, von Drehle M, Muth AN, et al: Dysregulation of cardiogenesis, cardiac conduction, and cell cycle in mice lacking miRNA-1-2. Cell. 2007, 129: 303-317. 10.1016/j.cell.2007.03.030.View ArticlePubMedGoogle Scholar
- Lagos-Quintana M, Rauhut R, Yalcin A, Meyer J, Lendeckel W, Tuschl T: Identification of tissue-specific microRNAs from mouse. Current Biology. 2002, 12: 735-739. 10.1016/S0960-9822(02)00809-6.View ArticlePubMedGoogle Scholar
- Esau C, Davis S, Murray SF, Yu XX, Pandey SK, Pear M, et al: miR-122 regulation of lipid metabolism revealed by in vivo antisense targeting. Cell Metabolism. 2006, 3: 87-98. 10.1016/j.cmet.2006.01.005.View ArticlePubMedGoogle Scholar
- Glazov EA, Cottee PA, Barris WC, Moore RJ, Dalrymple BP, Tizard ML: A microRNA catalog of the developing chicken embryo identified by a deep sequencing approach. Genome Research. 2008, 18: 957-964. 10.1101/gr.074740.107.PubMed CentralView ArticlePubMedGoogle Scholar
- Yu HH, Moens CB: Semaphorin signaling guides cranial neural crest cell migration in zebrafish. Developmental Biology. 2005, 280: 373-385. 10.1016/j.ydbio.2005.01.029.View ArticlePubMedGoogle Scholar
- Martyn U, Schulte-Merker S: Zebrafish neuropilins are differentially expressed and interact with vascular endothelial growth factor during embryonic vascular development. Developmental Dynamics. 2004, 231: 33-42. 10.1002/dvdy.20048.View ArticlePubMedGoogle Scholar
- Wilm TP, Solnica-Krezel L: Essential roles of a zebrafish prdm1/blimp1 homolog in embryo patterning and organogenesis. Development. 2005, 132: 393-404. 10.1242/dev.01572.View ArticlePubMedGoogle Scholar
- Miller CT, Maves L, Kimmel CB: moz regulates Hox expression and pharyngeal segmental identity in zebrafish. Development. 2004, 131: 2443-2461. 10.1242/dev.01134.View ArticlePubMedGoogle Scholar
- Kubista M, Andrade J, Bengtsson M, Forootan A, Jona J, Lind K, et al: The real-time polymerase chain reaction. Mol Aspects Med. 2006, 27: 95-125. 10.1016/j.mam.2005.12.007.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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.