Skip to main content
  • Research article
  • Open access
  • Published:

Characterization of bovine miRNAs by sequencing and bioinformatics analysis

Abstract

Background

MicroRNAs (miRNAs) are a family of ~22 nucleotide small RNA molecules which regulate gene expression by fully or partially binding to their complementary sequences in mRNAs or promoters. A large number of miRNAs and their expression patterns have been reported in human, mouse and rat. However, miRNAs and their expression patterns in live stock species such as beef cattle are not well studied.

Results

We constructed and sequenced small-RNA libraries to yield a total of 13,541 small-RNA sequences from 11 bovine tissues including brain, subcutaneous fat, muscle, liver, kidney, spleen and thymus. In total, 228 miRNAs including 29 novel miRNA candidates were identified. Of the 199 miRNAs, 101 have been previously reported as bovine miRNAs and the other 98 are bovine orthologs of known miRNAs that have been identified in at least one other mammalian species. Of the 29 novel miRNA candidates, 17 appeared at this point in time to be bovine specific, while the remaining 12 had evidence of evolutionary conservation in other mammalian species. Five miRNAs (miR-23a, -23b, -99a, -125b and -126-5p) were very abundant across the 11 tissues, accounting for 44.3% of all small RNA sequences. The expression analysis of selected miRNAs using qRT-PCR also showed that miR-26a and -99a were highly expressed in all tissues, while miR-122 and miR-133a were predominantly expressed in liver and muscle, respectively.

Conclusion

The miRNA expression patterns among 11 tissues from beef cattle revealed that most miRNAs were ubiquitously expressed in all tissues, while only a few miRNAs were tissue specific. Only 60% miRNAs in this study were found to display strand bias, suggesting that there are some key factors for mature miRNA selection other than internal stability. Most bovine miRNAs are highly conserved in other three mammalian species, indicating that these miRNAs may have a role in different species that are potential molecular markers for evolution.

Background

MicroRNAs (miRNAs) are small non-coding RNA molecules of approximately 22 nucleotides in length, which play important regulatory roles in animals and plants [1, 2]. miRNAs have been found to down-regulate the expression of target genes by binding to the complementary sites in transcripts and causing translational repression or transcript degradation [3]. Numerous biological processes in animal development, apoptosis, fat metabolism and hematopoietic differentiation have been reported to be regulated by miRNAs [4–7]. In addition, recent studies have revealed that miRNAs can increase protein translation by binding to complementary promoter sequences, extending the important function of miRNA to protein expression [8, 9].

Many experimental techniques and computational methods have been developed to identify miRNAs [10–12]; and a large number of miRNAs have been identified in primates, rodents, birds, fish and plants [13–16]. However, the number of miRNAs from bovine species is limited with only 125 reported (miRBase 12.0, December 2008), comparing with 866 from human and 627 from mouse. Early studies suggested that most miRNAs are conserved among related species [17, 18]. However, recent studies have shown that many newly identified miRNAs are species specific [19], and the expression of miRNAs is not strictly conserved among species [20]. Although many miRNAs are known to be differentially expressed during development and across tissue types [4–6, 21], little is known about the relative abundance and specificity of expression patterns among tissues for most bovine miRNAs. In this study, we profiled bovine miRNAs and evaluated their expression patterns from 11 beef cattle tissues including muscle, kidney, liver, spleen, thymus, three fat tissues and three brain tissues. Elucidation of the expression patterns of different miRNAs among different tissues will contribute to the understanding of the roles of miRNAs in gene expression regulatory networks for particular biological functions in livestock species.

Results

Identification of miRNAs and novel miRNA candidates

From 11 small RNA libraries, 11,880 clones were sequenced and 13,541 small RNA sequences were obtained with 1,000 to 1,400 small RNA sequences covered from each library. Among the total of 12,869 sequences, 11,189 sequences were classified as previously reported bovine miRNAs, 1,578 sequences were orthologs of known miRNAs from other mammalian species, and 102 sequences were novel miRNA candidates. The rest sequences were found to be other small non-coding RNA, mRNA, tRNA or rRNA, etc. (Table 1, Additional file 1).

Table 1 An overview of sequencing results from 11 small RNA libraries

In total, we identified 101 previously reported bovine miRNAs and 98 orthologs of known miRNAs from other mammalian species, which are predicted to be encoded by 191 different miRNA genes (Additional file 2). In all libraries but longissimus dorsi muscle (LDM) library, the previously reported bovine miRNA sequences represented ~75% to 95% of the total sequences of each library. Nine out of 101 known bovine miRNAs and 23 out of 98 orthologous miRNAs were found to be present only once in the total sequenced clones (Additional file 3). Orthologs of known miRNAs sequences from LDM library accounted for approximately one third of all orthologous miRNAs (Table 1). Two orthologous miRNAs (bta-miR-574-3p and -652) were detected more than two times and in at least two different tissues but were not mapped to the bovine genome (ucsc_btau4) (Additional file 3). All but two of the bovine orthologous miRNAs had orthologs in human species. These two miRNAs were found to be produced from different arms of the same pre-miRNA (bta-mir-1388) and conserved in platypus and horse (Figure 1). In addition, we identified 29 novel miRNA candidates with 22 of them found to be present only once in all sequenced clones (Table 2).

Table 2 Novel miRNA candidates identified in this study
Figure 1
figure 1

Conserved mir-1388 in three mammalian species. Mature miRNA sequences from bovine (bta) and platypus (oan) are showed in uppercase. Future name is shown as record ID for horse orthologs of mir-1388. Identical base is labeled in black or grey.

To investigate evolutionary conservation of these miRNA candidates, six vertebrate genomes including human, mouse, rat, dog, horse and chicken were taken into account. Twelve out of 29 miRNA candidates were found to be conserved in at least one of these species. Three miRNA candidates (bta-un10, bta-un13 and bta-un20) were identified in the 5' arm of bta-mir-381, -495 and -487a, respectively, while only the 3' arms of these miRNAs were reported in several species (Table 2, Additional file 4).

MicroRNA genes

A total of 207 genes encoding the identified bovine miRNAs were predicted using reported bovine miRNAs from miRBase 12.0 and the identified bovine orthologs of other mammalian miRNAs. Of these, 117 were reported in miRBase (version 12.0); 19 miRNA precursors were not reported in miRBase database previously; 71 were detected based on the orthologous miRNAs (Additional file 2). More than 30 bovine miRNAs were found to be encoded from more than one predicted hairpin precursor.

When the bovine miRNAs were compared to those from human, mouse and rat, most bovine mature miRNAs were also reported in the other three mammalian species. This agreed with that many miRNAs are conserved [17, 18] among mammalian species. Most of the conserved miRNAs in the four species were also found to have the same number of precursor(s) although they were encoded by multi-copies of genes (For example, let-7a as shown in Additional file 5), indicating that the processing of these miRNAs has been conserved among these species.

Sequence variations and end variants in mature miRNA sequences

Approximately 20% of the identified miRNA sequence were found to have mismatches with the bovine genome that were caused by post-transcriptional modification, and/or RT-PCR, and sequencing errors (Additional file 3). The prevailing sequence alterations such as 3' terminal A or U additions and A-to-G transitions could result from post-transcriptional modification as previously reported by Landgraf [13]. A-to-G transitions potentially caused by A-to-I editing were mainly identified in miR-99a and miR-376c (Figure 2), which have been reported to undergo A-to-I editing in human and mouse, respectively [22, 23].

Figure 2
figure 2

Positions of edited adenosines in two bovine miRNAs. Mature miRNA sequences are indicated in red letters. Edited adenosines are highlighted in yellow. The number of edited miRNA sequences and the number of all sequences of this miRNA are showed in brackets (edited/total).

It has been reported that a variant of end (size) different siRNAs or miRNAs could be generated during Dicer processing of double-stranded RNA (dsRNA), short hairpin RNA (shRNA) and miRNA precursors in vitro [24, 25]. Alignment of the identified sequences revealed that end variants were present in most miRNAs, showing that many miRNAs sequenced more than two times had end variants that were apparently generated from the same precursor (Additional file 3). Most end variants differed at the 3' end nucleotide(s) and a small percentage of them differed at the 5' end. The finding of the variant miRNA sequences is consistent with previous large scale miRNA profiling studies [13, 14]. More than 10 bovine miRNAs reported in miRBase were one or two nucleotides shorter at the 3' end comparing to ours, for example, bta-miR-18a and bta-miR-16b (Additional files 2 &3).

Analysis of miRNA expression in bovine tissues

Hierarchical clustering of 11 bovine tissues based on the relative cloning frequencies of miRNAs showed that miRNA expression was similar in the related tissues such as abdominal subcutaneous fat and rump subcutaneous fat, while different tissues had diverse expression profiles (Figure 3). Comparison of miRNA expression profiles among tissues revealed that very few miRNAs expression was tissue specific (e.g., miR-9, -124 in brain, miR-122 in liver, miR-1, miR-133a and -206 in muscle). Several miRNAs were found to be highly expressed in particular tissues: miR-204, -218, and -129-2-3p in brain tissues, miR-30a, -30e, -30d, -200a and -200b in kidney, miR-192 in liver, miR-451 in spleen, miR-21 in spleen and thymus, miR-193b, -378 in LDM muscle. In contrast, 15 miRNAs were identified in all tissues and several of them (e.g., miR-23b and -99a) were expressed at high levels in all tissues. Another 31 miRNAs were found in 8 or more tissues. The five most abundant miRNAs across the 11 tissues were miR-23a, -23b, -99a, -125b and -126-5p, accounting for 44.3% of all small RNA sequences (Additional file 3).

Figure 3
figure 3

Hierarchical clustering of tissues and miRNAs using Pearson correlation. Heatmap was constructed based on the relative cloning frequencies of miRNAs. Only those miRNAs occuring at least 20 times in all libraries are shown. Green and red indicate low frequency or high frequency miRNA cloned in the library, respectively.

To validate above miRNA expression patterns, quantitative RT-PCR was performed on tissue-specific miRNAs (miR-122, -133a), high cloning frequency miRNAs (miR-26a, -99a and -150) and low cloning frequency miRNAs (miR-103, -107, -411, -423-5p, -574-3p and -652). Two animals were tested for each tissue and similar RT-PCR results were obtained from both animals (data not show). The expression profiles of the 11 miRNAs across 11 tissues confirmed that miR-26a and -99a expressed at high levels in all tissues, while miR-122 and -133a exclusively expressed in liver and muscle, respectively. Four bovine miRNAs (miR-411, -423-5p, -574-3p and -652) were expressed at low levels in all tissue (Figure 4). The qRT-PCR data showed similar expression trend of detected expression of these miRNAs from the library except miR-103.

Figure 4
figure 4

Expression patterns of 11 miRNAs in 11 tissues. qRT-PCR results from 11 miRNAs across 11 tissues are compared. miRNA expression profiles are normalized with RNU6B. High expressed miRNA has low ΔCt.

Discussion

MiRNA expression patterns in different tissues have been profiled in several vertebrates [13, 26–30]. Although there were diverse expression profiles in different tissues, most miRNAs were ubiquitously expressed. Our comparison of miRNA expression across 11 tissues from bovine revealed a few tissue specific miRNAs: miR-9, -124 in brain, miR-122 in liver, miR-1, miR-133a and -206 in muscle, which had been previously reported in mouse and human [13, 27]. Brain and muscle tissues have much more specific or enriched miRNAs than other tissues especially fat tissues, indicating that these miRNAs may play important regulatory roles in these tissues.

MiRNAs may be gained or lost during evolution, however, many miRNAs are conserved [17, 18]. Identification of common miRNA precursors between bovine and other three mammalian species (human, mice, and rat) suggests these miRNAs may have similar roles as those in other species since these conserved miRNA have higher expression levels [14, 31]. The higher expression level of three conserved miRNAs (miR-26a, -99a and -150) in all tested bovine tissues suggest that these miRNA may be more relevant to the highly conserved biological process in mammalians. Further study to discover their regulatory functions are needed. Recent developed hhigh-throughput sequencing analysis has allowed the identification of an increasing number of species-specific miRNAs [12, 32] since these miRNA may play a role in host-specific biological process. Most of the miRNA candidates identified in this study were bovine specific, although their expression was less than 1% comparing the total sequenced clones, suggesting further studies using deep-sequencing technologies [12, 14, 32] or specialized small RNA isolation and cloning procedures [19, 33] may help to identify and understand the functions of species specific miRNAs.

MiRNAs are firstly transcribed as pri-miRNAs and gave rise to short, 70-nucleotide stem-loop structures (pre-miRNAs) by the Drosha-DGCR8 complex [34, 35]. The hairpin structures are then processed by Dicer. During the process, both strands of a miRNA precursor can form functional miRNAs. However, most of the time only one arm becomes the mature miRNA or predominant miRNA, while the other degrades or generates star miRNA. Khvorova et al [36] suggested that the arm with lower thermodynamic stability at its 5' end becomes the mature miRNA. In this study, we identified 15 pairs of bovine miRNAs and star miRNAs. The stability of the initial four base pairs of these miRNA pairs were calculated using nearest-neighbour method and 2-state hybridization algorithm [37, 38]. Sixty percent of bovine miRNAs (9/15) displayed strand bias (Additional file 6). We also calculated the stability of human miRNAs and star miRNAs from miRBase 12.0 using the same method. In total, 175 pairs of miRNAs and star miRNAs were taken into account and 15 pairs of them were excluded from evaluation. Interestingly, only 60% (96/160) of human miRNAs exhibit strand bias, too (Additional file 7). We argue that there must be additional important factors other than internal stability to determine which arm of the miRNA precursor becomes the mature miRNA or miRNA* since the following observations can not be explained: (1) Most miRNAs observed in this study had a variant of isoforms generated by Dicer and a few of the 5' end variants even processed by Drosha (e.g., four bta-miR-23a variants had an additional A nucleotide at the 5' end, Additional file 3). (2) Some conserved pre-miRNAs express different mature miRNAs or star miRNAs depending on the species. For example, the stem-loop sequence of bta-miR-126 was perfectly matched to those from human, mouse and rat; however, in cattle and mouse, both strands were observed as mature miRNAs, while in human and rat, one strand generates miRNA and the other strand generates miRNA*. (3) Some mature miRNAs from the same precursors reverse to star miRNAs or vice versa in different tissues or development stage [39].

In addition, to understand the roles of miRNA, we mapped and identified Chromosome-X related miRNA. We identified five to seven miRNA clusters containing ~20 common miRNAs in bovine, mouse and human, respectively. Cluster 4 and the sequences of corresponding orthologs of these miRNAs were not found in mouse genome. We conjectured this cluster should be on some gap of chr.X in mouse. Two miRNA clusters (5 & 6) in bovine were not found on the X-chromosome but were found on the contig Un.004.53 (Figure 5, Additional file 8). Interestingly, another common miRNA (miR-652) was not mapped to the bovine genome but was also identified in the X-chromosome of human and mouse (Additional file 8). We speculate that contig Un.004.53 and bta-miR-652 belong to chromosome-X. The identification of miRNA and miRNA clusters on chromosome-X reveals that some miRNAs are conserved even in genome location between species. The BLAST search of the miRNAs on the bovine genome identified a third precursor, which was a perfect match with bta-mir-138-2 in contig Un.004.5037. However, only two precursors of miR-138 were found in human, rat and mouse. Conserved number copy of miRNA genes between species suggests that bta-mir-138-3 is the same as bta-mir-138-2 and contig Un.004.5037 may belong to chromosome 22. Therefore, miRNAs and their clustering appearance maybe provide potential molecular markers for evolution.

Figure 5
figure 5

Organization of X-linked miRNA Clusters in huaman (A), mouse (B) and bovine (C). Five to seven X-linked miRNA clusters contained ~20 miRNAs were identified from bovine, mouse and human, respectively. miRNA clusters are from ~250 bp (cluster2) to 4.9 kb (cluster5) in size. Cluster 4 was not identified in mouse. Arrow indicated the direction of transcription.

Conclusion

Our small-RNA cloning and sequencing approach was efficiently target mature miRNAs. However, to identify and understand the functions of bovine miRNAs, further studies using deep-sequencing technologies are required. The miRNA expression patterns among 11 tissues in beef cattle showed most miRNAs are ubiquitously expressed, suggesting that these miRNAs may play a role in a broad range of biological processes in various tissues. Although mature miRNAs have been reported to display strand bias, the mechanism which arm of pre-miRNA dominates still need further study. Most identified bovine miRNAs were found in other three mammalian species and are highly conserved during evolution, suggesting that these miRNAs may have similar function in mammalian species and can be potential molecular markers for evolution.

Methods

Tissues collection and RNA extraction

Bovine tissues including longissimus dorsi muscle, kidney, liver, spleen, thymus, cerebellum, medulla, hypothalamus, abdominal subcutaneous fat and rump subcutaneous fat were collected from two 16-moth-old female Angus cross-breed cattle at the University of Alberta farm, while back subcutaneous fat tissue was collected from two 16-moth-old male Charolais cattle. Tissues were immediately frozen in liquid nitrogen and stored at -80°C until use. Total RNA and highly enriched small RNA (<200 nt) from tissues were extracted using mirVana miRNA Isolation Kit (Ambion Inc, Austin, TX, USA) according to the manufacturer's instruction.

Small RNA library preparation and sequencing

Eleven cDNA libraries for small RNAs from the bovine tissues were constructed using miRCatâ„¢ Kit (IDT DNA Technologies, Coralville, IA, USA). Briefly, ~30 ug of enriched small RNA for each library was separated on 15% denatured polyacrylamide gels. The RNAs between 18 and 24 nt in size were recovered from the gels and ligated to the 3' linker using T4 RNA ligase. The ligated RNAs (~42 nt) were eluted from 15% denatured polyacrylamide gels and ligated to the 5' linker. The 5'-and 3' -linked RNAs (~62 nt) were reverse-transcribed to cDNAs. The cDNA products were amplified and the purified PCR products were digested by BanI and finally concatemerized and cloned into pCR2.1 TOPO vector using a TOPO cloning kit (Invitrogen, Carlsbad, CA, USA). The white colonies were randomly selected from the libraries and some of them were analyzed for inserts by PCR. The plasmids were extracted by MILLIPORE vacuum system and Montage Plasmid Miniprep Kit and sequenced by 3730 Sequencing Analyzer.

Small RNA sequence analysis

A bioinformatics pipeline was developed for miRNA analysis. The small RNA sequences 18-26 in size were extracted and mapped to the latest bovine genome assembly (ucsc_btau4). All sequences were searched for miRNA sequences in miRBase. Small RNA sequence that had less than two mismatches (or >90% identity) with mammalian miRNAs in miRBase 12.0 was considered as a homologous miRNA and sorted by name without species prefix. Homologous miRNAs were further divided into previously reported bovine miRNAs and bovine orthologs of known miRNAs. Other small RNA sequences (mRNA, rRNA, tRNA, snoRNA, piRNA and misc RNA) were filtered by using BLAST against the NCBI nr database http://www.ncbi.nlm.nih.gov/.

The genomic DNA sequences flanking the bovine miRNAs from miRBase 12.0, orthologs of known miRNAs and small RNA sequences mapped to the bovine genome but not aligned to any annotated RNA classes were obtained from the bovine genome assembly. Sequences were analyzed for hairpin structure using the UNAFold algorithm with default parameters or the Mfold web server at http://mfold.bioinfo.rpi.edu/cgi-bin/rna-form1.cgi[40]. The hairpin structure containing putative miRNA was regarded as pre-miRNA if it contained neither large internal loops nor bulges, and at least 16 nt of the putative miRNA (if putative miRNA sequence length <22 nt, stretch 3' end to 22 nt) were paired with the other arm of the hairpin[18].

To investigate conservation of bovine novel miRNA candidates, highly similar DNA sequences corresponding to their precursors in the human, dog, mouse, rat, horse and chicken genome assemblies were found using BLATN NCBI GENOMES with word size = 7http://blast.ncbi.nlm.nih.gov/Blast.cgi[41]. Sequence alignments covering > 70% identity and >75% of the length of precursors of bovine miRNA candidates were used for hairpin structure analyses as described above. The orthologs of bovine miRNA candidates were defined to be those who have proper hairpin structures and less than 3 nucleotides variation sequences with those in other species.

Calculation of thermodynamic stability

Bovine miRNAs were sorted by clone counts. If the clone count of one strand was 10 times more than that of the other strand, they were regarded as miRNA and miRNA*, respectively. Human star miRNAs and their correspondingly high-expressed miRNAs were downloaded from miRBase (version 12.0). We used two-state hybridization web sever to calculate the free energy at http://dinamelt.bioinfo.rpi.edu/twostate.php. The internal stability of the initial four bases of miRNA or miRNA* (5'→3') were calculated as the delta G or delta G* (Additional files Table S5 & S6).

Quantitative RT-PCR

Quantitative RT-PCR was performed using human TaqMan miRNA probes that had the same sequence as the bovine miRNAs. Reactions were performed following manufacturer's recommendations (Applied Biosystems, Foster City, CA, USA) except for 20 ng total RNA used.

Heatmap construction

Hierarchical clustering of miRNA expression was performed with PermutMatrix [42], using Pearson distance, average linkage and normalized columns (Z score). Clustering was performed only on miRNA sequenced at least 20 times. The relative cloning frequency for each miRNA was calculated as the number of sequences for each miRNA in a library divided by the total number of miRNA sequences for that library.

Abbreviations

Asf:

abdominal subcutaneous fat

Bsf:

back subcutaneous fat

Rsf:

rump subcutaneous fat

Cer:

cerebellum

Hyp:

hypothalamus

Med:

medulla

LDM:

longissimus dorsi muscle

Kid:

kidney

Liv:

liver

Spl:

spleen

Thy:

thymus.

References

  1. Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5

    Article  CAS  PubMed  Google Scholar 

  2. Mallory AC, Vaucheret H: Functions of microRNAs and related small RNAs in plants. Nat Genet. 2006, 38 (Suppl): S31-36. 10.1038/ng1791

    Article  CAS  PubMed  Google Scholar 

  3. Pillai RS, Bhattacharyya SN, Filipowicz W: Repression of protein synthesis by miRNAs: how many mechanisms?. Trends Cell Biol. 2007, 17 (3): 118-126. 10.1016/j.tcb.2006.12.007

    Article  CAS  PubMed  Google Scholar 

  4. Song L, Tuan RS: MicroRNAs and cell differentiation in mammalian development. Birth Defects Res C Embryo Today. 2006, 78 (2): 140-149. 10.1002/bdrc.20070

    Article  CAS  PubMed  Google Scholar 

  5. Lee CT, Risom T, Strauss WM: MicroRNAs in mammalian development. Birth Defects Res C Embryo Today. 2006, 78 (2): 129-139. 10.1002/bdrc.20072

    Article  CAS  PubMed  Google Scholar 

  6. Wienholds E, Plasterk RH: MicroRNA function in animal development. FEBS Lett. 2005, 579 (26): 5911-5922. 10.1016/j.febslet.2005.07.070

    Article  CAS  PubMed  Google Scholar 

  7. Alvarez-Garcia I, Miska EA: MicroRNA functions in animal development and human disease. Development. 2005, 132 (21): 4653-4662. 10.1242/dev.02073

    Article  CAS  PubMed  Google Scholar 

  8. Place RF, Li L-C, Pookot D, Noonan EJ, Dahiya R: MicroRNA-373 induces expression of genes with complementary promoter sequences. PNAS. 2008, 105 (5): 1608-1613. 10.1073/pnas.0707594105

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Vasudevan S, Tong Y, Steitz JA: Switching from Repression to Activation: MicroRNAs Can Up-Regulate Translation. Science. 2007, 318 (5858): 1931-1934. 10.1126/science.1149460

    Article  CAS  PubMed  Google Scholar 

  10. Berezikov E, Cuppen E, Plasterk RH: Approaches to microRNA discovery. Nat Genet. 2006, 38 (Suppl): S2-7. 10.1038/ng1794

    Article  CAS  PubMed  Google Scholar 

  11. Aravin A, Tuschl T: Identification and characterization of small RNAs involved in RNA silencing. FEBS Lett. 2005, 579 (26): 5830-5840. 10.1016/j.febslet.2005.08.009

    Article  CAS  PubMed  Google Scholar 

  12. Berezikov E, van Tetering G, Verheul M, Belt van de J, van Laake L, Vos J, Verloop R, Wetering van de M, Guryev V, Takada S, et al: Many novel mammalian microRNA candidates identified by extensive cloning and RAKE analysis. Genome Res. 2006, 16 (10): 1289-1298. 10.1101/gr.5159906

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, Pfeffer S, Rice A, Kamphorst AO, Landthaler M, et al: A Mammalian microRNA Expression Atlas Based on Small RNA Library Sequencing. Cell. 2007, 129 (7): 1401-1414. 10.1016/j.cell.2007.04.040

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. 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 Res. 2008, 18 (6): 957-964. 10.1101/gr.074740.107

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Lagos-Quintana M, Rauhut R, Meyer J, Borkhardt A, Tuschl T: New microRNAs from mouse and human. RNA. 2003, 9 (2): 175-179. 10.1261/rna.2146903

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Chen PY, Manninga H, Slanchev K, Chien M, Russo JJ, Ju J, Sheridan R, John B, Marks DS, Gaidatzis D, et al: The developmental miRNA profiles of zebrafish as determined by small RNA cloning. Genes & Dev. 2005, 19 (11): 1288-1293. 10.1101/gad.1310605

    Article  CAS  Google Scholar 

  17. Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T: Identification of Novel Genes Coding for Small Expressed RNAs. Science. 2001, 294 (5543): 853-858. 10.1126/science.1064921

    Article  CAS  PubMed  Google Scholar 

  18. Ambros V, Bartel B, Bartel DP, Burge CB, Carrington JC, Chen X, Dreyfuss G, Eddy SR, Griffiths-Jones SAM, Marshall M, et al: A uniform system for microRNA annotation. RNA. 2003, 9 (3): 277-279. 10.1261/rna.2183803

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Bentwich I, Avniel A, Karov Y, Aharonov R, Gilad S, Barad O, Barzilai A, Einat P, Einav U, Meiri E, et al: Identification of hundreds of conserved and nonconserved human microRNAs. Nat Genet. 2005, 37 (7): 766-770. 10.1038/ng1590

    Article  CAS  PubMed  Google Scholar 

  20. Ason B, Darnell DK, Wittbrodt B, Berezikov E, Kloosterman WP, Wittbrodt J, Antin PB, Plasterk RHA: From the Cover: Differences in vertebrate microRNA expression. PNAS. 2006, 103 (39): 14385-14389. 10.1073/pnas.0603529103

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM: Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005, 433 (7027): 769-773. 10.1038/nature03315

    Article  CAS  PubMed  Google Scholar 

  22. Kawahara Y, Zinshteyn B, Sethupathy P, Iizasa H, Hatzigeorgiou AG, Nishikura K: Redirection of Silencing Targets by Adenosine-to-Inosine Editing of miRNAs. Science. 2007, 315 (5815): 1137-1140. 10.1126/science.1138050

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Blow MJ, Grocock RJ, van Dongen S, Enright AJ, Dicks E, Futreal PA, Wooster R, Stratton MR: RNA editing of human microRNAs. Genome Biol. 2006, 7 (4): R27- 10.1186/gb-2006-7-4-r27

    Article  PubMed Central  PubMed  Google Scholar 

  24. Zhang H, Kolb FA, Jaskiewicz L, Westhof E, Filipowicz W: Single processing center models for human Dicer and bacterial RNase III. Cell. 2004, 118 (1): 57-68. 10.1016/j.cell.2004.06.017

    Article  CAS  PubMed  Google Scholar 

  25. Vermeulen A, Behlen L, Reynolds A, Wolfson A, Marshall WS, Karpilow JON, Khvorova A: The contributions of dsRNA structure to Dicer specificity and efficiency. RNA. 2005, 11 (5): 674-682. 10.1261/rna.7272305

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Lagos-Quintana M, Rauhut R, Yalcin A, Meyer J, Lendeckel W, Tuschl T: Identification of tissue-specific microRNAs from mouse. Curr Biol. 2002, 12 (9): 735-739. 10.1016/S0960-9822(02)00809-6

    Article  CAS  PubMed  Google Scholar 

  27. Shingara J, Keiger K, Shelton J, Laosinchai-Wolf W, Powers P, Conrad R, Brown D, Labourier E: An optimized isolation and labeling platform for accurate microRNA expression profiling. RNA. 2005, 11 (9): 1461-1470. 10.1261/rna.2610405

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Xu H, Wang X, Du Z, Li N: Identification of microRNAs from different tissues of chicken embryo and adult chicken. FEBS Lett. 2006, 580 (15): 3610-3616. 10.1016/j.febslet.2006.05.044

    Article  CAS  PubMed  Google Scholar 

  29. Coutinho LL, Matukumalli LK, Sonstegard TS, Van Tassell CP, Gasbarre LC, Capuco AV, Smith TPL: Discovery and profiling of bovine microRNAs from immune-related and embryonic tissues. Physiol Genomics. 2007, 29 (1): 35-43.

    Article  CAS  PubMed  Google Scholar 

  30. Murchison EP, Kheradpour P, Sachidanandam R, Smith C, Hodges E, Xuan Z, Kellis M, Grutzner F, Stark A, Hannon GJ: Conservation of small RNA pathways in platypus. Genome Res. 2008, 18 (6): 995-1004. 10.1101/gr.073056.107

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Hicks JA, Tembhurne P, Liu HC: MicroRNA Expression in Chicken Embryos. Poult Sci. 2008, 87 (11): 2335-2343. 10.3382/ps.2008-00114

    Article  CAS  PubMed  Google Scholar 

  32. Berezikov E, Thuemmler F, van Laake LW, Kondova I, Bontrop R, Cuppen E, Plasterk RH: Diversity of microRNAs in human and chimpanzee brain. Nat Genet. 2006, 38 (12): 1375-1377. 10.1038/ng1914

    Article  CAS  PubMed  Google Scholar 

  33. Takada S, Berezikov E, Yamashita Y, Lagos-Quintana M, Kloosterman WP, Enomoto M, Hatanaka H, Fujiwara S-i, Watanabe H, Soda M, et al: Mouse microRNA profiles determined with a new and sensitive cloning method. Nucleic Acids Res. 2006, 34 (17): e115- 10.1093/nar/gkl653

    Article  PubMed Central  PubMed  Google Scholar 

  34. Gregory RI, Yan KP, Amuthan G, Chendrimada T, Doratotaj B, Cooch N, Shiekhattar R: The Microprocessor complex mediates the genesis of microRNAs. Nature. 2004, 432 (7014): 235-240. 10.1038/nature03120

    Article  CAS  PubMed  Google Scholar 

  35. Han J, Lee Y, Yeom K-H, Kim Y-K, Jin H, Kim VN: The Drosha-DGCR8 complex in primary microRNA processing. Genes & Dev. 2004, 18 (24): 3016-3027. 10.1101/gad.1262504

    Article  CAS  Google Scholar 

  36. Khvorova A, Reynolds A, Jayasena SD: Functional siRNAs and miRNAs exhibit strand bias. Cell. 2003, 115 (2): 209-216. 10.1016/S0092-8674(03)00801-8

    Article  CAS  PubMed  Google Scholar 

  37. Markham NR, Zuker M: DINAMelt web server for nucleic acid melting prediction. Nucleic Acids Res. 2005, 33 (suppl_2): W577-581. 10.1093/nar/gki591..

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Mathews DH, Sabina J, Zuker M, Turner DH: Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure. J Mol Biol. 1999, 288 (5): 911-940. 10.1006/jmbi.1999.2700

    Article  CAS  PubMed  Google Scholar 

  39. Grimson A, Srivastava M, Fahey B, Woodcroft BJ, Chiang HR, King N, Degnan BM, Rokhsar DS, Bartel DP: Early origins and evolution of microRNAs and Piwi-interacting RNAs in animals. Nature. 2008, 455 (7217): 1193-1197. 10.1038/nature07415

    Article  CAS  PubMed  Google Scholar 

  40. Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucl Acids Res. 2003, 31 (13): 3406-3415. 10.1093/nar/gkg595

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Caraux G, Pinloche S: PermutMatrix: a graphical environment to arrange gene expression profiles in optimal linear order. Bioinformatics. 2005, 21 (7): 1280-1281. 10.1093/bioinformatics/bti141

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by Alberta Livestock Industry Development Fund (2007F044R). The authors thank Mrs. Y. Meng for her technical report and Dr. U. Basu for her assistance in qRT-PCR analysis.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Le Luo Guan.

Additional information

Authors' contributions

LLG, and WJ designed the experiments; WJ performed library construction, sequence analysis, and original manuscript writing; JRG and PS developed bioinformatics pipeline for data analyses. SM involved in bovine genome mapping discussion and manuscript construction. WJ and LLG wrote the paper.

Electronic supplementary material

12867_2009_458_MOESM1_ESM.pdf

Additional file 1: Table S1- Description of 11 Small RNA libraries from bovine tissues. miRNAs include reported bovine miRNAs and orthologs of known miRNAs from other mammalian species (PDF 54 KB)

12867_2009_458_MOESM2_ESM.xls

Additional file 2: Table S2 - miRNAs and miRNA hairpins identified from bovine. New bovine miRNAs identified in this study are filled in gray; miRNAs identified in this study but with a different size than those reported previous are filled in green; miRNAs not-mapped bovine genome and their hairpins listed here from human are filled in yellow. (XLS 90 KB)

12867_2009_458_MOESM3_ESM.xls

Additional file 3: Table S3 - Distributions of miRNA and novel miRNA candidates in 11 bovine tissues. miRNA sequences mismatched with their genome were regarded as modification (Md) if its 3' terminal has A or U additions, or the same site of miRNA, more than 2 sequences found A to G mutation (A-G transition). Other random mismatches were considered as mutations (Mt). (XLS 202 KB)

12867_2009_458_MOESM4_ESM.pdf

Additional file 4: Supplemental File S1 - Hairpin structures of 29 novel bovine miRNA candidates and their corresponding vertebrate orthologs. Hairpin structure of novel bovine miRNA candidates and their corresponding vertebrate orthologs are listed together. Bovine miRNA candidates and their corresponding vertebrate orthologs are shown in uppercase. Variations between bovine miRNA candidates and vertebrate orthologs are indicated in red letters. For vertebrate orthologs, species (hsa, human; mmu, mouse; rno, rat; cfa, dog; eca, horse; gga, chicken), NCBI Accession Number, Genomic location are shown. (PDF 144 KB)

12867_2009_458_MOESM5_ESM.xls

Additional file 5: Table S4 - Comparison of bovine miRNA genes with those from other three mammalian species. 1: Genes for miRNA and related miRNA counted together are highlighted in yellow; "/" means "or". 2: All miRNA genes were from miRBase 12.0 and our study. 3: Star miRNAs were not listed here. (XLS 52 KB)

12867_2009_458_MOESM6_ESM.xls

Additional file 6: Table S5 - Thermodynamic stability of bovine miRNA and miRNA*. The internal stability of the initial four bases of miRNA or miRNA* (5'→3') were calculated as the delta G or delta G*, using two-state hybridization web sever (http://dinamelt.bioinfo.rpi.edu/twostate.php) (XLS 80 KB)

12867_2009_458_MOESM7_ESM.xls

Additional file 7: Table S6 - Thermodynamic stability of human miRNA and miRNA*. The internal stability of the initial four bases of miRNA or miRNA* (5'→3') were calculated as the delta G or delta G*, using two-state hybridization web sever http://dinamelt.bioinfo.rpi.edu/twostate.php. (XLS 568 KB)

12867_2009_458_MOESM8_ESM.pdf

Additional file 8: Table S7 - X-linked miRNAs in bovine, mouse and human. X-linked miRNAs from bovine, mouse and human, were clustered if they were less than 5,000 bp apart on the same strand. (PDF 98 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Jin, W., Grant, J.R., Stothard, P. et al. Characterization of bovine miRNAs by sequencing and bioinformatics analysis. BMC Molecular Biol 10, 90 (2009). https://doi.org/10.1186/1471-2199-10-90

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2199-10-90

Keywords