Identification and profiling of conserved and novel microRNAs from Chinese Qinchuan bovine longissimus thoracis
© Sun et al.; licensee BioMed Central Ltd. 2013
Received: 1 July 2012
Accepted: 2 January 2013
Published: 18 January 2013
MicroRNAs (miRNAs) are a family of ~22 nucleotide small RNA molecules that regulate gene expression by fully or partially binding to their complementary sequences. Recently, a large number of miRNAs and their expression patterns have been identified in various species. However, to date no miRNAs have been reported to modulate muscle development in beef cattle.
Total RNAs from the Chinese Qinchuan bovine longissimus thoracis at fetal and adult stages were used to construct small RNA libraries for Solexa SBS technology sequencing. A total of 15,454,182 clean reads were obtained from the fetal bovine library and 13,558,164 clean reads from the adult bovine library. In total, 521 miRNAs including 104 novel miRNA candidates were identified. Furthermore, the nucleotide bias, base edit and family of the known miRNAs were also analyzed. Based on stem-loop qPCR, 25 high-read miRNAs were detected, and the results showed that bta-miRNA-206, miRNA-1, miRNA-133, miRNAn12, and miRNAn17 were highly expressed in muscle-related tissue or organs, suggesting that these miRNAs may play a role in the development of bovine muscle tissues.
This study confirmed the authenticity of 417 known miRNAs, discovered 104 novel miRNAs in bos taurus, and identified five muscle-specific miRNAs. The identification of novel miRNAs significantly expanded the repertoire of bovine miRNAs and could contribute to further studies on the muscle development of cattle.
KeywordsBovine Deep sequencing technology microRNA Muscle Proliferation Differentiation
MicroRNAs (miRNAs) are a new class of single-stranded endogenous non-coding small RNA molecules (~22 nucleotides) that bind primarily to the 3′UTR of target mRNAs to repress their translation and accelerate their decay , and may regulate up to 30% of genes . The majority of miRNAs are conserved across species and play essential roles in regulating many distinct processes such as brain morphogenesis , insulin secretion , virus immune defense , metabolism, and signal transduction . In addition, several studies have also revealed the significance of miRNAs in myocyte proliferation and differentiation in recent years . For example, miR-1 and miR-133 have distinct roles in modulating skeletal muscle proliferation and differentiation in cultured myoblasts in vitro and in Xenopus laevis embryos in vivo. miR-1 promotes myogenesis by targeting histone deacetylase 4 (HDAC4), a transcriptional repressor of muscle gene expression. In contrast, miR-133 enhances myoblast proliferation by repressing serum response factor (SRF) . MiR-206, miR-1, and miR-133 are muscle specific miRNAs . The transcription of miR-206 is induced by MyoD, which promotes myogenic differentiation . Bone Morphogenetic Protein-2 (BMP-2), which is known to inhibit myogenesis, represses the expression of miR-206 by inhibiting its maturation process . Similarly, overexpression of miR-181 during muscle differentiation is important to promote myogenesis by down-regulating the homeobox protein Hox-A11, an inhibitor of muscle differentiation . MiR-486 has also been shown to induce myoblast differentiation by down-regulating Pax7 , while MiR-27b regulates Pax3 protein levels and ensures myogenic differentiation . Recently, the myoblast cell line C2C12 was used for functional analysis of miR-214 in vitro. The results showed that miR-214 may target the negative regulators of Myf5, MyoD and myogenin in the corresponding stages of skeletal muscle development in vivo to regulate embryonic myogenesis. To date, miRNAs have become one of the most abundant categories of gene regulatory molecules in mammalian species, but the role of individual miRNAs in muscle development is still unknown.
Currently 19,724 mature miRNAs have been discovered from 153 species and deposited in the publicly available miRNA database miRBase (Release 17.0, April 2011) (http://www.mirbase.org). Specifically, the number of miRNAs from bovine species is limited with only 674 reported, compared with 1,921 from human and 1,157 from mouse. Despite the recognized importance of miRNAs in regulating gene expression during development and other biological processes in beef, there has been little information about miRNA expression in cattle. This is surprising given that cattle have tremendous importance not only for food production but as a mammalian model organism for comparative genomics and biological studies . Recent related studies have been conducted to provide insight into the miRNA population present in bovine species by investigating the characteristics, expression pattern and features of their target genes. For instance, 59 distinct miRNAs were identified from bovine adipose tissue and mammary gland in 2006 alone . Bta-mir424 and bta-mir-10b are highly abundant in germinal vesicle oocytes, as well as in early stage embryos (until 16-cell stage) . MiR-196a is a bona fide negative regulator of the newborn ovary homeobox gene (NOBOX) during bovine early embryogenesis . Expression of bovine nucleoplasmin 2 (NPM2) is temporally regulated during early embryogenesis by miR-181a . Approximately 20% of the miRNAs involved in adipogenesis and lipid deposition were identified as being correlated with backfat thickness. The results suggest that miRNAs play a regulatory role in white adipose tissue development in beef .
Given the emerging roles of miRNAs in development, identifying the differentially expressed miRNAs is an important first step to investigating the function of miRNAs in the course of bovine growth and development. In this study, we have constructed two small RNA cDNA libraries from Chinese Qinchuan bovine longissimus thoracis at fetal and adult stages. By high throughput sequencing of the small RNA libraries and subsequent bioinformatic analysis, miRNAs in the longissimus thoracis were identified. In addition, the expression patterns of the high-read miRNAs from differential tissues (muscle, heart, liver, lung, kidney, brain, intestine, fat, and spleen) at multiple developmental stages of bovine muscle tissues (fetal, calf, and adult) were evaluated. Elucidation of the expression patterns of different miRNAs among different tissues will contribute to understanding the roles of miRNAs in gene expression regulatory networks for particular biological functions in livestock species.
Results and discussion
Tissue collection and high-throughput sequencing of small RNAs
Skeletal muscle is composed of myofibers, intramuscular adipocytes and connective tissue. Myofibers are the structural units of skeletal muscle . In livestock, all muscle fibers are formed during the prenatal stage. Bovine prenatal myogenesis can be briefly divided into three different generations of cells, which appear at around 60, 90, and 110 days of fetal life (post-conception) . In contrast, postnatal skeletal muscle development is mainly due to the increase in muscle fiber size , and new muscle fibers are only generated during the adult stage to replace injured muscle fibers . This pattern is significantly different between prenatal and postnatal bovine muscle development. Hence, in this study fetal and adult Chinese Qinchuan bovine longissimus thoracis were collected and two miRNA libraries were constructed for Solexa SBS technology sequencing.
Summary of small RNA sequencing date
Fetal bovine muscle tissue
Adult bovine muscle tissue
Distribution of the genome-mapped sequence reads in small RNA libraries
Fetal bovine muscle tissue
Adult bovine muscle tissue
Identification of conserved bovine miRNAs
Summary of known miRNA in each sample
Unique matched to pre-miRs
Read matched to pre-miRs
Fetal bovine muscle tissue
Adult bovine muscle tissue
Position 2–8 of a mature miRNA is called the seed region, which is highly conserved. The target of a miRNA might be different with the change of nucleotides in this region. In our analysis pipeline, miRNAs that might have base edit can be detected by aligning unannotated sRNA tags with mature miRNAs from miRBase17, allowing one mismatch at a certain position. The results showed that approximately 23.41% in the adult bovine library and 29.06% in the fetal bovine library of the identified miRNA sequences were found to have mismatches that were caused by post-transcriptional modification, and/or RT-PCR, and sequencing errors (Additional files 5 and 6). The prevailing sequence alterations such as 3' terminal A or U additions and A-to-G transitions could result from post-transcriptional modification . In addition, alignment of the identified sequences revealed that end variants were present in most miRNAs, which were apparently generated from the same precursor. Most end variants differed by one or several nucleotides at the 3' end nucleotide and a small percentage of them differed at the 5' end. For example, miR-105b, miR-1185, and miR-122 only had 3' end variants. miR-142a, miR-1271, and miR-138 differed by only one nucleotide in the 5' end, but they had several 3' end variants (Additional file 7). Variation in end (size) in miRNAs may result from the processing of double-stranded RNA (dsRNA), short hairpin RNA (shRNA) and miRNA precursors by Dicer [34, 35].
To date, 2,042, 1,281, 791, 289, 360, 103, 247, and 306 conserved miRNAs were identified in human, mouse, chicken, dog, horse, sheep, zebrafish and wild boar, respectively http://www.mirbase.org/. In our study, further analysis identified a total of 289 and 407 conserved miRNAs that belonged to 200 miRNA families in adult and fetal bovine libraries. The identified miRNA families have been shown to be conserved in a variety of species. For example, let-7, miR-1, miR-34, miR-9, and miR-25 families have been found in 64, 61, 56, 65 and 62 species, respectively, while miR-2363, miR-2384, miR-2404, miR-2424, and miR-409 families have only been detected in bos taurus (Additional file 8). This may suggest a species-specific expression profile for miRNAs. The largest miRNA family size identified was miR-2284, which consisted of 12 members, and let-7, miR-30, and miR-181/376 possessed 9, 7, and 4 members, respectively; whereas other miRNA families such as miR-1, miR-31, miR-93, and miR-206 had only one member (Additional file 1). The finding that most members of conserved miRNA families were expressed in the bovine longissimus thoracis supports the idea that regulatory or functional diversification has occurred [36, 37]. Different family members also displayed drastically different expression levels. For example, the abundance of the miR-2284 family varied from 1 read (bta-miR-2284a) to 46,548 reads (bta-miR-2284x) with deep sequencing. This was also the case for some other miRNA families, such as bta-let-7 (from 6 to 1,434,682 reads), bta-miR-30 (from 34 to 12,681 reads) and bta-miR-181 (from 376 to 20,258 reads). However, the expression levels of some miRNA families were similar, such as miR-15 in which 246 and 276 reads were detected, respectively. The existence of a dominant member in a miRNA family may suggest that the regulatory role of this family was performed by the dominant member at the developmental time when the samples were collected for RNA extraction. Abundance comparisons of different members in a miRNA family may provide valuable information on the role that miRNAs play in that specific stage of bovine development.
Identification of novel bovine miRNAs
The characteristic hairpin structure of a miRNA precursor can be used to predict novel miRNA. Prediction software Mireap http://sourceforge.net/projects/mireap/ was developed to predict novel miRNA by exploring the secondary structure, the Dicer cleavage site and the minimum free energy of the unannotated small RNA tags that could be mapped to a genome. The following criteria  were used for screening the candidates for potential miRNAs or pre-miRNAs: (1) Pre-miRNA sequences can fold into an appropriate hairpin secondary structure that contains the ~22nt mature miRNA sequence within one arm of the hairpin. (2) miRNA precursors with secondary structures had higher negative minimal free energies (MFEs) and minimal free energy indexes (MFEIs) than other different types of RNAs. (3) miRNA had an AU content of 30–70%. (4) miRNA had less than six mismatches with the opposite miRNA* sequence in the other arm. (5) No loop or break in miRNA sequences was allowed. Based on Solexa sequencing, we identified 104 novel bovine miRNAs, which corresponded to 145 genomic loci. Thirty-six novel miRNAs were in the adult bovine library and 92 were in the fetal bovine library, of which 24 overlapped in both libraries (Additional file 9). In addition, an examination of pre-miRNAs and other RNAs (tRNA, rRNA, and mRNA) revealed that miRNAs were significantly different from other RNAs . Specifically, more than 90% of miRNA precursors have an MFEI greater than 0.85, significantly higher than tRNAs (0.64), rRNAs (0.59), or mRNAs (0.65). The results suggested that the MFEI can easily be used to distinguish miRNA from other non-coding and coding RNAs. This provides a more precise criterion to predict miRNAs using computational approaches, and in our database 122 had a MFEI greater than 0.85. Remarkably, the read number for each novel miRNA was much lower than that for the majority of conserved miRNAs. For instance, bta-miRn70 with the highest read number and bta-miRn96 with the lowest were both novel miRNAs, and the total reads were only 541,928 and 5, respectively. Recently, similar results in pigs  and even in maize  have also been reported.
Validation of related bovine miRNAs
To further explore the muscle-specific miRNAs involved in the muscle development of the cattle, we performed quantitative analysis of the miRNAs in fetal (day 90 bovine embryos), calf (3-day-old) and adult (2-year-old) bovine longissimus thoracis. The results showed that expression of bta-miRNA-133 and miRn12 was increased in the muscle tissues from day 90 bovine embryos to 2-year-olds, respectively. However, the expression levels of bta-miRNA-206, miRNA-1, and miRn17 did not change between day 90 bovine embryo and 2-year muscle tissues, but significantly increased in the calf muscle tissue (Figure 3). Previous studies in vitro have shown that miR-1, miR-133, and miR-206 can target multiple muscle-development-related genes. Specifically, muscle-specific miR-206, which is directly activated by MyoD, can target sequences in the Fstl1 and Utrn gene and these sequences are sufficient to suppress gene expression in the presence of miR-206 . miR-1 promotes myogenesis by targeting HDAC4, a transcriptional repressor of muscle gene expression. In contrast, miR-133 enhances myoblast proliferation by repressing SRF . Also, miR-1 and miR-206 regulate Pax7 directly. Inhibition of these two substantially enhances satellite cell proliferation and increases Pax7 protein levels in vivo. Although the bovine-specific target genes of miRNA-206, miRNA-1, miRNA-133, miRn12, and miRn17 are not known, their consistent expression pattern and high conservation indicate that they are also likely to play roles in the development of bovine muscle tissues.
We have identified 417 known miRNAs and 104 novel miRNAs in longissimus thoracis from fetal and adult Qinchuan bovine using deep sequencing technologies. This study expands the repertoire of bovine miRNAs and could initiate further study in the muscle development of cattle. In addition, the miRNA expression patterns among nine tissues in beef cattle showed that most miRNAs are ubiquitously expressed, suggesting that these miRNAs may play a role in a broad range of biological processes in various tissues. However, we have also identified some muscle-specific miRNAs, suggesting that these miRNAs are likely to play a role in the development of bovine muscle tissues and could be potential molecular markers for genetics and breeding.
All animals in this study were maintained according to the No. 5 proclamation of the Ministry of Agriculture, P. R. China. Sample collection was approved by the Animal Care Commission of the College of Animal Science and Technology, Northwest A&F University. Bovine embryos of slaughtered cows were collected from Tumen abattoir, a local slaughterhouse of Xi'An, P.R. China. A newborn Qinchuan calf and adult Qinchuan cattle were obtained from Meixian Qinbao Co., Ltd.
Tissue collection and high-throughput sequencing
Day 90 (d90) bovine embryos (gestation period 280 days) were collected into sterile physiological saline immediately after removal from the reproductive tract of slaughtered cows at a local abattoir. Fetal age was estimated based on crown-rump length . Bovine tissue samples including the longissimus thoracis, heart, liver, lung, kidney, cortex (brain), small intestine, fat and spleen were collected from fetal, calf and adult Chinese Qinchuan bovine. These tissues were snap-frozen in liquid nitrogen and stored at −80°C until use. In this study, two miRNA libraries were constructed. Total RNAs were extracted from three fetal and three adult Chinese Qinchuan bovine longissimus thoracis were pooled, respectively. Subsequently, low molecular weight RNAs were separated by 15% polyacrylamide gel electrophoresis (PAGE), and RNA molecules in the range of 18–30nt were enriched and ligated with proprietary adapters to the 5′ and 3′ termini. A reverse transcription reaction followed by low cycle PCR was performed to obtain sufficient product for Solexa technology (Beijing Genomics Institute, China).
Small RNA sequence analysis
After clearing away the 3′ adaptor sequence, removal of redundancy and reads smaller than 18nt, the clean reads were screened against and mapped to the latest bovine genome assembly http://hgdownload.cse.ucsc.edu/goldenPath/bosTau4/bigZips/bosTau4.fa.gz using the program SOAP . To identify sequences originating from protein-coding genes, repeats, rRNA, tRNA, snRNA, and snoRNA, we used bovine mRNA http://hgdownload.cse.ucsc.edu/goldenPath/bosTau4/database/refGene.txt.gz and CDS http://hgdownload.cse.ucsc.edu/goldenPath/bosTau4/bigZips/refMrna.fa.gz, RepeatMasker http://www.repeatmasker.org and Sanger Rfam data (version 10.1). Subsequently, the remaining reads were searched against the Sanger miRBase (version 18.0) to identify the conserved miRNAs. Only those small RNAs whose mature and precursor sequences perfectly matched known bovine miRNAs in miRBase were considered to be conserved miRNAs. To discover potential novel miRNA precursor sequences, unique sequences that have more than 10 hits to the genome or match to known non-coding RNAs were removed. Then the flanking sequences (150 nt upstream and downstream) of each unique sequence were extracted for secondary structure analysis with Mfold http://www.bioinfo.rpi.edu/applications/mfold and then evaluated by Mireap http://sourceforge.net/projects/mireap/. Specifically, the miRNA candidates that passed Mireap were deemed as highly probable if their corresponding miRNA*s were also found in the small RNA libraries. After prediction, the resulting potential miRNA loci were examined carefully based on the distribution and numbers of small RNAs on the entire precursor regions. Those sequences residing in the stem region of the stem-loop structure and ranging between 20–22nt with free energy hybridization lower than −20 kcal/mol were considered .
MicroRNA expression analysis
The x and y represented normalized expression levels, and the N1 and N2 represented total count of clean reads of a given miRNA in small RNA libraries of the fetal and adult stage, respectively. Stem-loop real-time reverse transcription polymerase chain reaction (RT-PCR) with SYBR Green was used for the analysis of miRNA expression . Total RNA (1 μg from tested tissues) was converted to cDNA with a RT primer mixture (250 nM) using PrimeScript® RT reagent Kit with gDNA Eraser (TaKaRa, Dalian, China). The cDNA was then used for real-time PCR quantification of miRNA using the miRNA specific primer and the universal primer. The bovine ribosomal protein S18 (RPS18) (GenBank NO. NM_001033614.1) gene was used as an endogenous control. The primers for miRNAs and the control gene are listed in Additional file 11. Real-time quantitative PCR was performed using a Bio-Rad CFX 96™ Real Time Detection System and SYBR Green PCR Master Mix (TaKaRa, Dalian, China) in a 20 μl reaction. All reactions were carried out in triplicate. The PCR mix included 100 ng cDNA for each miRNA, 0.4 μM forward and reverse primers, and 10 μl 2 × SYBR Green PCR Master Mix. The cycle conditions were as follows: 95°C for 30 s, followed by 40 cycles of 95°C for 10 s, 60°C for 10 s, and 68°C for 20 s. The threshold cycle (Ct) was defined as the cycle number at which the fluorescence intensity passed a predetermined threshold. The quantification of each miRNA relative to RPS18 gene was calculated using the equation: N = 2-ΔΔCt.
This study was supported by the National Natural Science Foundation of China (Grant Nos. 31272408, 30972080), Agricultural Science and Technology Innovation Projects of Shaanxi Province (No. 2012NKC01-13), Program of National Beef Cattle Industrial Technology System (CARS-38), and Natural Science Foundation of Jiangsu Province (BK2011206).
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.View ArticlePubMedGoogle Scholar
- Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115: 787-798. 10.1016/S0092-8674(03)01018-3.View ArticlePubMedGoogle Scholar
- Giraldez A, Cinalli R, Glasner M, Enright A, Thomson J, Baskerville S, Hammond S, Bartel D, Schier A: MicroRNAs regulate brain morphogenesis in zebrafish. Science. 2005, 308: 833-838. 10.1126/science.1109020.View ArticlePubMedGoogle Scholar
- Poy MN, Eliasson L, Krutzfeldt J, Kuwajima S, Ma X, Macdonald PE, Pfeffer S, Tuschl T, Rajewsky N, Rorsman P, Stoffel M: A pancreatic islet-specific microRNA regulates insulin secretion. Nature. 2004, 432: 226-230. 10.1038/nature03076.View ArticlePubMedGoogle Scholar
- Lecellier CH, Dunoyer P, Arar K, Lehmann-Che J, Eyquem S, Himber C, Saib A, Voinnet O: A cellular microRNA mediates antiviral defense in human cells. Science. 2005, 308: 557-560. 10.1126/science.1108784.View ArticlePubMedGoogle Scholar
- Ambros V: The functions of animal microRNAs. Nature. 2004, 431: 350-355. 10.1038/nature02871.View ArticlePubMedGoogle Scholar
- Zhao Y, Samal E, Srivastava D: Serum response factor regulates amuscle-specific microRNA that targets Hand2 during cardiogenesis. Nature. 2005, 436: 214-220. 10.1038/nature03817.View ArticlePubMedGoogle Scholar
- Jian-Fu C, Mandel EM, Michael Thomson J, Qiulian W, Callis TE, Hammond SM, Conlon FL, Da-Zhi W: The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation. Nat Genet. 2005, 38: 228-233.Google Scholar
- McCarthy JJ: MicroRNA-206: the skeletal muscle-specific myomiR. Biochim Biophys Acta. 2008, 1779: 682-691. 10.1016/j.bbagrm.2008.03.001.PubMed CentralView ArticlePubMedGoogle Scholar
- Missiaglia E, Shepherd CJ, Patel S, Thway K, Pierron G, Pritchard-Jones K, Renard M, Sciot R, Rao P, Oberlin O, Delattre O, Shipley J: MicroRNA-206 expression levels correlate with clinical behaviour of rhabdomyosarcomas. Br J Cancer. 2010, 102: 1769-1777. 10.1038/sj.bjc.6605684.PubMed CentralView ArticlePubMedGoogle Scholar
- Sato MM, Nashimoto M, Katagiri T, Yawaka Y, Tamura M: Bone morphogenetic protein-2 down-regulates miR-206 expression by blocking its maturation process. Biochem Biophys Res Commun. 2009, 383: 125-129. 10.1016/j.bbrc.2009.03.142.View ArticlePubMedGoogle Scholar
- Naguibneva I, Polesskaya A, Ameyar-Zazoua M, Souidi M, Groisman R, Cuvellier S, Ait-Si-Ali S, Pritchard LL, Harel-Bellan A: Micro-RNAs and muscle differentiation. J Soc Biol. 2007, 201: 367-376. 10.1051/jbio:2007902.View ArticlePubMedGoogle Scholar
- Dey BK, Gagan J, Dutta A: miR-206 and −486 induce myoblast differentiation by downregulating Pax7. Mol Cell Biol. 2011, 31: 203-214. 10.1128/MCB.01009-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Crist CG, Montarras D, Pallafacchina G, Rocancourt D, Cumano A, Conway SJ, Buckingham M: Muscle stem cell behavior is modified by microRNA-27 regulation of Pax3 expression. Proc Natl Acad Sci USA. 2009, 106: 13383-13387. 10.1073/pnas.0900210106.PubMed CentralView ArticlePubMedGoogle Scholar
- Feng Y, Cao JH, Li XY, Zhao SH: Inhibition of miR-214 expression represses proliferation and differentiation of C2C12 myoblasts. Cell Biochem Funct. 2011, 29: 378-383. 10.1002/cbf.1760.View ArticlePubMedGoogle Scholar
- Gibs R, Weinstock G, Kappes S, Schook L, Skow L, Womack J: Bovine Genomic Sequencing Initiative. 2006, 1-12. http://www.genome.gov/Pages/Research/Sequencing/SeqProposals/BovineSEQ.pdf,Google Scholar
- Zhiliang G, Satyanaryana E, Honglin J: Identification and characterization of microRNAs from the bovine adipose tissue and mammary gland. FEBS Lett. 2007, 581: 981-988. 10.1016/j.febslet.2007.01.081.View ArticleGoogle Scholar
- Tripurani SK, Xiao C, Salem M, Yao J: Cloning and analysis of fetal ovary microRNAs in cattle. Anim Reprod Sci. 2010, 120: 16-22. 10.1016/j.anireprosci.2010.03.001.View ArticlePubMedGoogle Scholar
- Tripurani SK, Lee KB, Wee G, Smith GW, Yao J: MicroRNA-196a regulates bovine newborn ovary homebox gene (NOBOX) expression during early embryogenesis. BMC Dev Biol. 2011, 11: 25-10.1186/1471-213X-11-25.PubMed CentralView ArticlePubMedGoogle Scholar
- Lingenfelter BM, Tripurani SK, Tejomurtula J, Smith GW, Yao J: Molecular cloning and expression of bovine nucleoplasmin 2 (NPM2): a maternal effect gene regulated by miR-181a. Reprod Biol Endocrinol. 2011, 9: 40-10.1186/1477-7827-9-40.PubMed CentralView ArticlePubMedGoogle Scholar
- Jin W, Dodson MV, Moore SS, Basarab JA, Guan LL: Characterization of microRNA expression in bovine adipose tissues: a potential regulatory mechanism of subcutaneous adipose tissue development. BMC Mol Biol. 2010, 11: 29-10.1186/1471-2199-11-29.PubMed CentralView ArticlePubMedGoogle Scholar
- Winter J, Jung S, Keller S, Gregory RI, Diederichs S: Many roads to maturity: microRNA biogenesis pathways and their regulation. Nat Cell Biol. 2009, 11: 228-234. 10.1038/ncb0309-228.View ArticlePubMedGoogle Scholar
- Sciote JJ, Morris TJ: Skeletal muscle function and fibre types: the relationship between occlusal function and the phenotype of jaw-closing muscles in human. J Orthod. 2000, 27: 15-30. 10.1093/ortho/27.1.15.View ArticlePubMedGoogle Scholar
- Hocquette JF: Endocrine and metabolic regulation of muscle growth and body composition in cattle. Animal. 2010, 4: 1797-1809. 10.1017/S1751731110001448.View ArticlePubMedGoogle Scholar
- Brameld JM, Mostyn A, Dandrea J, Stephenson TJ, Dawson JM, Buttery PJ, Symonds ME: Maternal nutrition alters the expression of insulin-like growth factors in fetal sheep liver and skeletal muscle. J Endocrinol. 2000, 167: 429-437. 10.1677/joe.0.1670429.View ArticlePubMedGoogle Scholar
- Li R, Li Y, Kristiansen K, Wang J: SOAP: short oligonucleotide alignment program. Bioinformatics. 2008, 24: 713-714. 10.1093/bioinformatics/btn025.View ArticlePubMedGoogle Scholar
- Du M, Yan X, Tong JF, Zhao J, Zhu MJ: Maternal obesity, inflammation, and fetal skeletal muscle development. Biol Reprod. 2010, 82: 4-12. 10.1095/biolreprod.109.077099.PubMed CentralView 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
- Evgeny AG, Pauline AC, Wesley CB, Robert JM, Brian PD, Mark LT: A microRNA catalog of the developing chicken embryo identified by a deep sequencing approach. Genome Res. 2008, 18: 957-964. 10.1101/gr.074740.107.View ArticleGoogle Scholar
- Aravin AA, Lagos QM, Yalcin A, Zavolan M, Marks D, Ben S, Terry G, Jutta M, Thomas T: The small RNA profile during Drosophila melanogaster development. Dev Cell. 2003, 5: 337-350. 10.1016/S1534-5807(03)00228-4.View ArticlePubMedGoogle Scholar
- Evgeny AG, Kritaya K, Wanchai A, Paul FH, Neena M, Timothy JM: Repertoire of Bovine miRNA and miRNA-Like Small Regulatory RNAs Expressed upon Viral Infection. PLoS One. 2009, 4 (7): e6349-10.1371/journal.pone.0006349.View ArticleGoogle Scholar
- Zhang B, Stellwag EJ, Pan X: Large-scale genome analysis reveals unique features of microRNAs. Gene. 2009, 443: 100-109. 10.1016/j.gene.2009.04.027.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Blow M, Grocock R, Van Dongen S, Enright A: RNA editing of human microRNAs. Genome Biol. 2006, 7: R27-10.1186/gb-2006-7-4-r27.PubMed CentralView ArticlePubMedGoogle Scholar
- Xie Z, Allen E, Fahlgren N, Calamar A, Givan SA, Carrington JC: Expression of Arabidopsis miRNA Genes. Plant Physiol. 2005, 138: 2145-2154. 10.1104/pp.105.062943.PubMed CentralView ArticlePubMedGoogle Scholar
- Miska EA, Alvarez-Saavedra E, Abbott AL, Lau NC, Hellman AB, McGonagle SM, Bartel DP, Ambros VR, Horvitz HR: Most caenorhabditis elegans microRNAs are individually not essential for development or viability. PLoS Genet. 2007, 3: e215-10.1371/journal.pgen.0030215.PubMed CentralView ArticlePubMedGoogle Scholar
- Xie FL, Huang SQ, Guo K, Xiang AL, Zhu YY, Nie L, Yang ZM: Computational identification of novel microRNAs and targets in Brassica napus. FEBS Lett. 2008, 581: 1464-1474.View ArticleGoogle Scholar
- Zhang BH, Pan XP, Cox SB, Cobb GP, Anderson TA: Evidence that miRNAs are different from other RNAs. Cell Mol Life Sci. 2006, 63: 246-254. 10.1007/s00018-005-5467-7.View ArticlePubMedGoogle Scholar
- Chen C, Deng B, Qiao M, Zheng R, Chai J, Ding Y, Peng J, Jiang S: Solexa sequencing identification of conserved and novel microRNAs in Backfat of Large White and Chinese Meishan Pigs. PLoS One. 2012, 7 (2): e31426-10.1371/journal.pone.0031426.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang L, Liu H, Li D, Chen H: Identification and characterization of maize microRNAs involved in the very early stage of seed germination. BMC Genomics. 2011, 12: 154-10.1186/1471-2164-12-154.PubMed CentralView ArticlePubMedGoogle Scholar
- Jin W, Grant JR, Stothard P, Moore SS, Guan LL: Characterization of bovine miRNAs by sequencing and bioinformatics analysis. BMC Mol Bio. 2009, 10: 90-View ArticleGoogle Scholar
- Wang XW: A PCR-based platform for microRNA expression profiling studies. RNA. 2009, 15: 716-723. 10.1261/rna.1460509.PubMed CentralView ArticlePubMedGoogle Scholar
- Rosenberg MI, Georges SA, Asawachaicharn A, Analau E, Tapscott SJ: MyoD inhibits Fstl1 and Utrn expression by inducing transcription of miR-206. J Cell Biol. 2006, 175: 77-85. 10.1083/jcb.200603039.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen JF, Tao Y, Li J, Deng ZL, Yan Z, Xiao X, Wang DZ: microRNA-1 and microRNA-206 regulate skeletal muscle satellite cell proliferation and differentiation by repressing Pax7. J Cell Biol. 2010, 190: 867-879. 10.1083/jcb.200911036.PubMed CentralView ArticlePubMedGoogle Scholar
- Richardson C, Jones PC, Barnard V, Hebert CN, Terlecki S, Wijeratne WV: Estimation of the developmental age of the bovine fetus and newborn calf. Vet Rec. 1990, 126: 279-284.PubMedGoogle Scholar
- Ambros V, Bartel B, Bartel DP, Burge CB, Carrington JC, Chen XM, Dreyfuss G, Eddy SR, Griffiths-Jones S, Marshall M: A uniform system for microRNA annotation. RNA. 2003, 9 (3): 277-279. 10.1261/rna.2183803.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR, Lao KQ, Livak KJ, Guegler KJ: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33: e179-10.1093/nar/gni178.PubMed CentralView ArticlePubMedGoogle Scholar