Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Highly Accessed Research article

MicroRNA transcriptome profiles during swine skeletal muscle development

Tara G McDaneld1, Timothy PL Smith1*, Matthew E Doumit2, Jeremy R Miles1, Luiz L Coutinho3, Tad S Sonstegard4, Lakshmi K Matukumalli5, Dan J Nonneman1 and Ralph T Wiedmann1

Author Affiliations

1 USDA/ARS Meat Animal Research Center, Clay Center, NE, USA

2 Department of Animal Science, Michigan State University, East Lansing, MI, USA

3 University of Sao Paulo-ESALQ, Piracicaba, SP, Brazil

4 USDA/ARS Beltsville Area Research Center, Beltsville, MD, USA

5 Bioinformatics and Computational Biology, George Mason University, Manassas, VA, USA

For all author emails, please log on.

BMC Genomics 2009, 10:77  doi:10.1186/1471-2164-10-77


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/10/77


Received:16 October 2008
Accepted:10 February 2009
Published:10 February 2009

© 2009 McDaneld et al; licensee 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.

Abstract

Background

MicroRNA (miR) are a class of small RNAs that regulate gene expression by inhibiting translation of protein encoding transcripts. To evaluate the role of miR in skeletal muscle of swine, global microRNA abundance was measured at specific developmental stages including proliferating satellite cells, three stages of fetal growth, day-old neonate, and the adult.

Results

Twelve potential novel miR were detected that did not match previously reported sequences. In addition, a number of miR previously reported to be expressed in mammalian muscle were detected, having a variety of abundance patterns through muscle development. Muscle-specific miR-206 was nearly absent in proliferating satellite cells in culture, but was the highest abundant miR at other time points evaluated. In addition, miR-1 was moderately abundant throughout developmental stages with highest abundance in the adult. In contrast, miR-133 was moderately abundant in adult muscle and either not detectable or lowly abundant throughout fetal and neonate development. Changes in abundance of ubiquitously expressed miR were also observed. MiR-432 abundance was highest at the earliest stage of fetal development tested (60 day-old fetus) and decreased throughout development to the adult. Conversely, miR-24 and miR-27 exhibited greatest abundance in proliferating satellite cells and the adult, while abundance of miR-368, miR-376, and miR-423-5p was greatest in the neonate.

Conclusion

These data present a complete set of transcriptome profiles to evaluate miR abundance at specific stages of skeletal muscle growth in swine. Identification of these miR provides an initial group of miR that may play a vital role in muscle development and growth.

Background

Functionally important small RNAs were first described in nematodes in 1993 [1,2]. However, it was not until 2001 that researchers began to understand the function of this family of RNAs that includes microRNA (miR) and to recognize that their significance was not confined to lower order organisms [3,4]. The small RNA classified as miR are short sequences, 18–26 nucleotide long, encoded by nuclear genes that produce characteristic stem-loop RNA structures when transcribed. During processing from the primary transcript, the mature miR sequence is loaded into an RNA:protein complex known as the "RNA induced silencing complex" (RISC) [5,6]. The sequence of the miR loaded in the complex targets the RISC to specific binding sites in the 3' untranslated region of mRNA transcripts, resulting in either degradation of the miR:mRNA complex or translocation to processing bodies. In either case, association of RISC with mRNA causes decreased translation of the targeted gene product [6,7]. As a result of decreased translation of their cognate targets, miR have been reported to guide developmental decisions including cell fate, cell cycle progression, apoptosis, adipocyte differentiation, and processes that alter muscle development and growth including myoblast proliferation, differentiation, and skeletal muscle hypertrophy [8-15].

The objective of the current research was to evaluate miR transcriptome profiles during skeletal muscle development in swine. MicroRNA were initially reported to have a role in skeletal muscle development utilizing mouse, drosophila, and zebrafish models. Three muscle-specific miRNA (miR-1, miR-133, and miR-206) were identified to increase in abundance during muscle cell differentiation [10,16,17]. However, these miRNA have been reported to regulate different stages of myogenesis [12,13,15]. MiR-133 increases proliferation of C2C12 myoblasts, whereas miR-206 and miR-1 promote differentiation [16]. Research in livestock models has begun to evaluate the role these miRNA have in skeletal muscle development. Expression of the muscle regulatory factor, myogenic factor 5, has been reported to regulate miR-1 and miR-206 transcription level in a chicken cell culture model [18]. In addition, over-expression of fibroblast growth factor-4 has been reported to decrease miR-206 abundance, resulting in developmental changes in the somite of developing chicken embryos [19]. Muscle-specific miR have also been reported to regulate a gene that directly impacts economic traits in livestock [20]. A mutation in the myostatin gene of heavily muscled Belgian Texel sheep creates a target site for miR-1 and miR-206 containing RISC complexes in the 3' untranslated portion of the transcript, resulting in decreased translation of the myostatin protein and consequent increase in muscle mass.

With the dramatic increase in identified miR sequences for multiple species including livestock species, a public database dedicated to the cataloguing of predicted and experimentally observed miRs has been developed (miRbase) [21-23]. In human, 678 miRs have been described (miRbase release 11.0, April 2008). However, genomic sequence scans and miR cloning results indicate that the actual total number of human miRs may be closer to 800 [24]. Comparative analysis of these sequences indicates that they are highly conserved among species as diverse as nematodes and mammals, supporting the hypothesis that they are of central importance to biological processes. In addition, expression of miR genes is tightly regulated spatially among tissues and temporally within tissues during development in all species studied, indicating the importance of determining miR transcriptome profiles to fully understand their biological importance [25-27]. In order to identify miR and determine their role in skeletal muscle of livestock, we evaluated miR transcriptome profiles at specific stages of muscle development including proliferating satellite cells, three stages of fetal growth, day-old neonate, and the adult.

Results and Discussion

Validation of miR cDNA libraries

The method of miR identification and quantification by cloning and sequencing has been utilized in numerous reports of miR biology [28]. It is particularly advantageous to use this method when working in species with poorly characterized genomes, since variation in miR sequences between species has become apparent [27,29,30].

Our initial goal was to obtain complete transcriptome profiles of miR abundance in skeletal muscle, and to accurately measure changes in the level of abundance for these miR between developmental stages. Time points evaluated included proliferating satellite cells (4th, 5th, and 6th passage), three stages of fetal development (60, 90, and 105 day-old fetus), day-old neonate, and adult. After identification of miR, accurate quantification of changes in miR abundance level between libraries and stages of skeletal muscle development was imperative. For all species, the dynamic range and precision obtained with enumerating each sequence observed enables a more quantitative description of miR abundance levels. Converting this precision to accuracy, requires running enough samples to saturate the signal until the abundance levels stay constant as new sequences are added. Initially, data from the neonatal muscle profiles were used to determine miR signal saturation as an indicator that abundance level of miR would remain constant as additional data were added (Figure 1). Although the number of unique singletons was far from exhausted after 5,000 observed putative miR, the number of new sequences observed for known and unknown miR began to plateau, suggesting that the singletons were either rare miR, represent contamination of the tissue with trace amounts of other tissues, or are sequence artifacts. Therefore, the data sets consisting of several thousand miR sequences were determined to be extensive enough to capture the diversity of miR abundance and estimate relative steady state levels in the developmental stages examined. Similar plots (not shown) were created from proliferating satellite cell profiles of the 4th and 5th passage. As with the neonate, abundance levels of known and unknown miR were well defined once a few thousand miR had been sequenced, as estimated abundance levels based on the first 3,000 sequenced clones from the 4th and 5th passage satellite cells agreed with estimates for the full library of 8,832 clones within 10%. Together, these results provide confidence that the number of clones sequenced provide an accurate representation of the miR transcriptome profiles of porcine skeletal muscle.

thumbnailFigure 1. Saturation plot of microRNA libraries. Saturation plots were created from the neonatal muscle sample to determine saturation of the signal as an indicator that abundance levels would remain constant as new results were added. A total of 5,000 observed miR were evaluated. The supply of unique singletons is represented by the small dashed line, while the known and unknown miR are represented by the solid black line and large dashed line, respectively.

Comparison of miR expression profiles between libraries created from independent samples at the same developmental state was of interest to examine biological and technical replication of the clone-based approach to miR profiling. Previous studies in other species have not specifically addressed this issue. Therefore, two approaches were completed to evaluate possible variation between miR libraries. First, transcriptome profiles were compared between miR libraries that were replicated in the current experiment including the three proliferating satellite cell, two adult biceps femoris, and two d 90 fetal biceps femoris libraries (Figure 2). As detailed in the Methods, RNA was obtained from multiple cell culture isolations of proliferating satellite cells at different cell passages including 4th, 5th, and 6th passages. The 4th and 5th passages were pooled together and compared to two different cell culture isolations from the 6th passage cells. The miR trascriptome profiles from the two 6th passage samples differed from each other as much as from the 4th and 5th passage sample (Figure 2a). Therefore, variability in abundance was quantified using a histogram (data not shown) of relative changes in abundance. Data utilized in the histogram was restricted to miR observed at a minimum abundance ratio level of 5 per thousand miR observed. This analysis resulted in 30 observed changes between samples that were not expected to have differing transcriptome profiles. Of these observed changes, the majority of changes were small (17 of the 30 ratios were between 1 and 2 per thousand miR observed, with an additional 8 between 2 and 3 per thousand miR observed). However, four ratios (13%) were between 3 and 6 per thousand miR expressed, with one almost reaching a ratio of 10. As a result of this analysis, we concluded that changes in abundance are likely to be significant if they are greater than 6-fold. Upon final evaluation of the miR transcriptome profiles, we took a more conservative approach and restricted our discussion to changes in miR abundance levels that were 10-fold or greater. Second, two individual libraries from the same 6th passage RNA ligation template were created to determine if variability was introduced during the PCR amplification step (see Methods). The transcriptome profiles of the replicate libraries were identical (data not shown), suggesting that variation observed previously between the three satellite cell libraries was due to satellite cell populations.

thumbnailFigure 2. Variation between microRNA libraries. Transcriptome profiles from samples at the same developmental state were compared within the three satellite cell libraries (2a), two individual adult biceps femoris libraries (2b), and two d 90 fetal biceps femoris libraries (2c). MicroRNA cDNA libraries for satellite cells (2a) were created from the pool of the 4th and 5th passage satellite cells and the 6th passage (6th passage, #1). A second library was created from a second flask of 6th passage stellite cells (6th passage, #2) to evaluate variation between satellite cells at the same passage. MicroRNA abundance is represented as number of individual miR tags observed per thousand tags evaluated.

Further evaluation of miR tags from the adult biceps femoris libraries (Figure 2b) demonstrated reproducible abundance levels between libraries of the moderate to highly abundant miR including let-7 and the muscle-specific miR-1 and miR-206. However, clear differences between the two libraries were evident as seven additional miR were identified in the first library compared to the second. In contrast to the adult and satellite cell miR libraries, the fetal libraries exhibited similarity across all miR that were present (Figure 2c). This observed variation in miR abundance within the adult libraries versus the fetal libraries could be attributed to animal variation. For the fetal libraries, biceps femoris of four female fetuses was collected and pooled for RNA extraction separately from two individual sows. The pooling strategy necessary to obtain sufficient starting material for the fetal muscle samples may have reduced variation present between individual samples and present a less variable overall miR profile. In contrast, each adult library was created from a muscle sample of an individual sow. Additionally, difference between libraries could be a result of variation in muscle sample. Muscle is not a homogeneous tissue [31], and it has been reported that gene abundance levels change based on location of the sample and distribution of Type I and Type II muscle fibers throughout the sample [32-36]. For the fetal libraries, the entire biceps femoris was obtained from four fetuses. The pooled sample was then powdered to a homogenous mix for RNA extraction. However, for the adult libraries a two to four gram sample was obtained for RNA extraction due to the greater size of the adult muscle. Therefore, it is possible that the two biceps femoris samples for the adult libraries were heterogeneous and may have contributed to the variation in miR tags between adult libraries.

MicroRNA transcriptome profiles in skeletal muscle

A digital transcriptome profile approach [37] was applied to evaluate miR abundance based on cloning the miR population from each sample and evaluating abundance as the number of transcripts for a given miR gene per thousand transcripts observed (see Additional file 1). This cloning and sequencing-based approach was highly successful, as indicated by the high degree of homology between our results and those previously described in miRbase for other species.

Additional file 1. Adundance levels of miR in skeletal muscle at specific developmental states. The data provided represent the transcription profiles of miR during specific stages throughout skeletal muscle development in swine.

Format: DOCX Size: 27KB Download fileOpen Data

Sequence comparison to known miR identified the muscle-specific miR-206 as the highest abundant miR across all muscle samples, which represented greater than 60% of all miR present (Figure 3). This contrasts sharply with the abundance of miR-206 in proliferating satellite cells at 1.8 per thousand (or 0.18%). In mouse C2C12 cells, miR-206 is also lowly abundant in proliferating cells and has been reported to be induced during differentiation [38], suggesting that its presence is associated with the switch from precursor to mature muscle cell. In addition, mR-1 and miR-133 are lowly abundant in proliferating C2C12 myoblasts [38]. The data herein confirm these results in satellite cells and also demonstrate that miR-206 is present at a high level through most or all of fetal development, and continues through maturation of the adult pig. The constant high-level presence at all stages following early differentiation suggests that the role of miR-206 is to repress functions associated with muscle precursor cells.

thumbnailFigure 3. MicroRNA transcriptome profiles. MicroRNA cDNA clone libraries were created from skeletal muscle during specific stages of swine development including satellite cells, d 60 (primary fiber development) of fetal development, d 90 (secondary fiber development) of fetal development, d 105 of fetal development, one day-old neonate and the adult. For the fetal samples, biceps femoris (BF) and longissimus dorsi (LD) were collected from female and male fetuses. Biceps femoris (BF) samples were collected for the one day-old neonate and the adult. MicroRNA abundance is represented as number of individual miR tags observed per thousand tags evaluated. Data for the satellite cells and adult muscle is presented as the average of multiple transcriptome libraries presented in Additional file 1.

In comparison to the high abundance of miR-206 in muscle tissue, porcine miR-1 had relatively moderate abundance that increased throughout development, similar to the pattern observed for this miR in mouse muscle development [16]. These data suggest that miR-1 and miR-206 play different roles in muscle development, with miR-1 affecting regulation of genes that require inactivation in later fetal stages and miR-206 having a more constant role in repressing genes immediately after differentiation. In contrast to miR-1 and miR-206, miR-133 was detected only at low levels in fetal development yet increased in the neonate and adult (Figure 3). Based on abundance level, these data suggest that miR-206 and miR-1 may have a greater role in fetal muscle development than miR-133, or their targets are higher in abundance. In addition, both miR-206 and miR-1 promote differentiation [16], suggesting that these two miR may have a greater impact compared to miR-133 on increased differentiation, which characterizes fetal development [39].

In addition to muscle-specific miR, a larger number of ubiquitous miR were present across all libraries. While a greater percentage of ubiquitous miR were lowly abundant throughout development, a number increased in abundance at specific developmental stages. MiR present 10 fold greater in satellite cells compared to neonate muscle included miR-16, miR-18, miR-27, miR-29, miR-34, and miR-106. Slightly below the factor-of-ten cut-off included miR-24 at 9.2. These miR have been implicated in multiple cellular processes including cell growth (miR-24), apoptosis (miR-16, miR-24 and miR-29), and cell cycle regulation in normal (miR-16) and cancerous cells (miR-24, miR-27, miR-29 and miR-34) [9,40-46]. Conversely, research has reported low abundance during proliferation and differentiation in cell culture for these miR with the exception of miR-24, which has a moderately high abundance level during late differentiation [16]. This difference in abundance level may be due in part to the different cell culture models utilized in the experiments. Chen et al. [16] evaluated miR in an immortalized C2C12 cell line, while the experiment herein utilized a porcine primary muscle cell line (satellite cells) during proliferation.

The fetal time points examined in swine were selected to coincide with important events in muscle development, specifically the waves of primary and secondary fiber formation (See Methods). Overall, miR were lowly abundant throughout fetal development with the exception of let-7 and muscle specific miR-1 and miR-206. Differential abundance of lowly to moderately abundant miR was observed between time points for primary (d 60) and secondary fiber development (d 90 and d 105; Figure 3). MiR-432 was moderately abundant during early fetal development at d 60, while miR-424 abundance increased during d 90 and d 105, and miR-126 abundance increased during later stages of fetal development at d 90 and 105. Currently, function of miR-432 in muscle development has not been determined. However, previous research reports that miR-424 regulates monocyte and macrophage differentiation [47], while miR-126 expression alters cell cycle progression of cancer cells by decreasing tumor growth and proliferation [48]. Four miR, miR-338, miR-368, miR-376, and miR-381, were also identified to be specific to the one-day old neonate as demonstrated by a ten-fold increase in abundance level, compared to both the adult muscle and satellite cells, suggesting that these miR may have a role in muscle growth immediately following birth. As for the adult, miR abundance was greatest for miR-151 and the muscle specific miR-1 and miR-133 compared to fetal and neonate miR libraries. While miR-1 and miR-133 abundance increases during differentiation in cell culture [16], the role of these miR in adult skeletal muscle has yet to be fully determined.

Classification of potential novel miR

Evaluation of sequence clusters identified two different classifications of novel miR; sequence tags that differed at only one (highly conserved) position and sequence tags that had no match to miR in the database. Five observed tags, miR-168a, miR-206, miR-24a, miR-368, and miR-381, represented possible exceptions to the pattern of exact sequence conservation across positions 4–17 of the reference miR. These five observed sequences differed from known miR at only one position, and the miR were observed between 23 and 59 times, and are not likely to be attributable to experimental artifacts. Three of the five were observed more frequently without mismatches: miR-206, miR-24a, and miR-168a (a cross-contamination from the parallel control oligo processing, see Methods). Interestingly, all of the miR-381 related tags displayed mismatches with the reported human miR-381 sequence with a single mismatch occurring 59 times (G?A at position 10 in Sus scrofa miR sequence), indicating that the Sus scrofa version of miR-381 does indeed differ in sequence from human and mouse. MiR-368 was observed 23 times as a mismatch and only 12 times without a mismatch. These five examples of single base mismatches were included in the transcriptome profiles along with those that matched exactly.

In addition to identification of mismatches in the miR-206 sequence, length of sequence at the 5' and 3' ends also differed (Table 1, see Additional file 2). While there are known instances of highly similar or identical mature miR being produced from discrete genes, it seems more likely that these differences are attributed to minor alterations in processing, cloning or sequencing of miR originating from the same gene. Since the goal was to produce relative transcriptome profiles of potential miR during muscle development, we clustered sequences of this type into a single miR sequence, using the most commonly observed sequence of the cluster as the defining sequence for comparisons to known miR.

Additional file 2. miR-206 sequence variation. The data provided represent the variation of miR-206 sequences identified in the current study.

Format: DOCX Size: 118KB Download fileOpen Data

Table 1. miR-206 sequence variation

Combining the data from all transcriptome profiles identified a total of 94 distinct miR that matched reference sequences and 12 sequences that did not match, and are therefore potential novel or porcine-specific miR (given temporary identifiers PN (porcine new) 1 to PN12; Table 2). While a greater number of these novel miR were lowly abundant throughout development, PN1 abundance increased in the proliferating satellite cells, neonate, and adult. Further validation of these putative novel miR awaits development of the swine genome to look for hallmarks of microRNA genes related to these tag sequences.

Table 2. Novel porcine miR

Computational identification of miR targets

Relatively few miRNA targets have been identified experimentally, but numerous computational predictions are readily available including miRanda, RNAhybrid and TargetScan [49-51]. Initially, miR targets were predicted for a sub-set of the miR up-regulated in this experiment (miR-206, miR-338, miR-368, miR-376a, and miR-381). These miR were predicted by miRNA viewer to target 47 "common genes" and 864 genes when the "all genes" option was chosen. For the purposes of this study, the predicted target genes were narrowed down to include only those with "muscle" listed in the gene ontology, as noted by the Entrezgene project [52]. A total of 19 genes were selected based on these criteria (see Additional file 3). Of these, 7 targets were identified for the highest expressed miR during fetal development to the adulthood, miR-206, including genes that have been implicated in multiple myogenic processes (see Additional file 3) [50,53-56]. These predicted targets for miR-206 include dystophia myotonica protein kinase, which has been implicated in myotonic dystrophy [55] and the transcription factor paired box gene 3 that regulates myogenic cell fate through the myogenic transcription factor MyoD [54]. Secondly, targets for a sub-set of down-regulated miR (miR-15, miR-16, miR-27, miR-29, miR-34 and miR-106), of which 44 targets were identified based on our previous criteria were predicted (see Additional file 4). From these identified gene targets, the ability to link the specific miR to target transcripts will improve as computational methods evolve, and as the databases of expressed and genomic porcine sequences grow.

Additional file 3. Predicted targets of up-regulated miR. The data provided represent the predicted gene targets for the up-regulated miR.

Format: DOCX Size: 74KB Download fileOpen Data

Additional file 4. Predicted targets of up-regulated miR. The data provided represent the predicted gene targets for the down-regulated miR.

Format: DOCX Size: 120KB Download fileOpen Data

Conclusion

Together, these data suggest that miR have a role in progression of myogenesis throughout development and their function may be specific to different stages of skeletal muscle growth. In addition, the data reported herein are the most complete set of transcriptome profiles to evaluate miR abundance in skeletal muscle at specific time points during fetal development of swine. Identification of these miR provide an initial group of expressed miR that change in abundance during specific developmental stages and therefore may target genes that regulate this process.

Methods

Skeletal muscle collection and clone libraries

MicroRNA libraries for the satellite cells were created from cells cultured from semimebranosus of 8-week old piglets as previously described [57] and incubated at 37°C with 5% CO2. Satellite cells at passage four, five, and six were collected for RNA extraction. RNA from the 4th and 5th passages was combined for creation of the first satellite cell cDNA library. RNA from two sets of 6th passage satellite cells were used for the second cDNA library. Tissues during fetal development were collected at d 60, 90, and 105 of fetal development [57] by removing porcine fetuses immediately after sacrifice and dissecting longissimus dorsi and biceps femoris muscle. Four female and four male fetuses were obtained from a single sow at each time point, and samples were pooled by sex and muscle type before immersion in liquid nitrogen. Neonatal biceps femoris was obtained at day one after birth, while two to four grams of adult biceps femoris was also obtained from an adult sow at slaughter. The experimental procedures were approved and performed in accordance with U. S. Meat Animal Research Center Animal Care Guidelines and the Guide for the Care and Use of Agricultural Animals in Agricultural Research and Teaching (1999). RNA extraction was performed using TRIreagent following the manufacturer's recommended protocol (Ambion, Austin, TX). Concentration and quality of RNA was determined using an Agilent 2100 Bioanalyzer for RNA (Agilent Technologies, Santa Clara, CA). Single insert cDNA libraries were constructed as described previously by Lu et al. [58] with the following modifications. RNA fractions were isolated from denaturing acrylamide gels and eluted by FlashPAGE (Ambion, Austin, TX). First strand cDNA synthesis was performed utilizing a primer to the 3' adapter sequence and SUPERSCRIPT reverse transcriptase (Invitrogen, Carlsbad, CA). The microRNA were then amplified by PCR and digested with EcoRI for ligation into pBLUEscript and electroporation in EC100 electrocompetent cells. Individual colonies were transferred into 384-well plates and grown in ampicillin selective LB media for plasmid preparation and sequencing using an Applied Biotechnology 3730 sequencer.

Statistical analysis

Chromatograms were converted into sequences and scored using Phred [59]. Sequences were collected based on identification of flanking vector and linker sequences. The intervening sequences were kept as putative miR, as long as their length was between 16 and 27 bases. The putative miR were then clustered based on sequence similarity into characteristic consensus sequences, where each member was required to match 14 consecutive bases to the most common member of the cluster. This approach was used due to frequent observance of highly similar miR differing only at their 3' ends, which often varied only by length or in the sequence of the last base [27]. The characteristic sequences were then compared to all known miR from the miRBase [21-23]. The criterion used for a positive match was that the putative miR contained an exact match to positions 4–17 of a known miR, as this segment of the miR sequence are highly conserved and unique to each miR. Validity of our clustering approach was tested using 455 known human miR. This resulted in 401 miR matching only to their unique sequence. Of the remaining 54 sequences, none involved the human homologue of the porcine miR reported in our current experiment. While it is possible that a portion of the clusters may represent miR from more than one distinct gene, the homology of the core targeting sequence indicates they are likely to have similar targets. Those that were not identified this way were screened, using BLAST, against tRNA, rRNA, snoRNA and mitochondrial sequence. The remaining unidentified sequences were checked for single base mismatches within positions 4–17 of previously identified miR. Only five examples were found, and they were counted with the full matching sequences. The observed putative miR that had at least two mismatches and that were observed at least 20 times were given temporary labels PN (porcine new) 1 to PN12 in order of decreasing levels of abundance.

MicroRNA were considered expressed if the tag cluster had at least ten members per thousand tags observed. Low abundance was defined as 0 to 16 tags per thousand observed, moderate abundance was defined as 16 to 256 tags, and high abundance was defined as greater than 256 tags (Figure 3). A difference between tag counts was accounted if greater than 10 fold.

Abbreviations

miR: microRNA; MiRNA: microRNA; mRNA: messenger RNA; PCR: polymerase chain reaction; PN: Porcine new; RISC: RNA induced silencing complex; RNA: ribonucleic acid

Authors' contributions

TGM conducted experimental design, tissue collection, RNA extraction, cDNA library construction, and participated in drafting of the manuscript. TPLS conceived the project and participated in experimental design and drafting of the manuscript. TPLS, LLC and TSD participated in training for development of cDNA libraries. MED collected and provided satellite cells. JRM participated in experimental design and tissue collection. LKM participated in data analysis. DJN provided tissue for day-old neonate. RTW conducted data analysis and participated in drafting the manuscript. All authors read and approved of the final manuscript.

Acknowledgements

The authors would like to recognize Renee Godtel, Bob Lee and Steve Simcox for technical assistance and Sherry Kluver for secretarial assistance.

References

  1. 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(5):843-854. PubMed Abstract | Publisher Full Text OpenURL

  2. Wightman B, Ha I, Ruvkun G: Posttranscriptional regulation of the heterochronic gene lin-14 by lin-4 mediates temporal pattern formation in C. elegans.

    Cell 1993, 75(5):855-862. PubMed Abstract | Publisher Full Text OpenURL

  3. Lau NC, Lim LP, Weinstein EG, Bartel DP: An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditis elegans.

    Science 2001, 294(5543):858-862. PubMed Abstract | Publisher Full Text OpenURL

  4. Lee RC, Ambros V: An extensive class of small RNAs in Caenorhabditis elegans.

    Science 2001, 294(5543):862-864. PubMed Abstract | Publisher Full Text OpenURL

  5. Hutvagner G: Small RNA asymmetry in RNAi: function in RISC assembly and gene regulation.

    FEBS Lett 2005, 579(26):5850-5857. PubMed Abstract | Publisher Full Text OpenURL

  6. Lai EC: miRNAs: whys and wherefores of miRNA-mediated regulation.

    Curr Biol 2005, 15(12):R458-R460. PubMed Abstract | Publisher Full Text OpenURL

  7. Chan SP, Slack FJ: microRNA-mediated silencing inside P-bodies.

    RNA Biol 2006, 3(3):97-100. PubMed Abstract | Publisher Full Text OpenURL

  8. Carthew RW: Gene regulation by microRNAs.

    Curr Opin Genet Dev 2006, 16(2):203-208. PubMed Abstract | Publisher Full Text OpenURL

  9. Garzon R, Fabbri M, Cimmino A, Calin GA, Croce CM: MicroRNA expression and function in cancer.

    Trends Mol Med 2006, 12(12):580-587. PubMed Abstract | Publisher Full Text OpenURL

  10. Brennecke J, Stark A, Cohen SM: Not miR-ly muscular: microRNAs and muscle development.

    Genes Dev 2005, 19(19):2261-2264. PubMed Abstract | Publisher Full Text OpenURL

  11. Kajimoto K, Naraba H, Iwai N: MicroRNA and 3T3-L1 pre-adipocyte differentiation.

    RNA 2006, 12(9):1626-1632. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Anderson C, Catoe H, Werner R: MIR-206 regulates connexin43 expression during skeletal muscle development.

    Nucleic Acids Res 2006, 34(20):5863-5871. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Nakajima N, Takahashi T, Kitamura R, Isodono K, Asada S, Ueyama T, Matsubara H, Oh H: MicroRNA-1 facilitates skeletal myogenic differentiation without affecting osteoblastic and adipogenic differentiation.

    Biochem Biophys Res Commun 2006, 350(4):1006-1012. PubMed Abstract | Publisher Full Text OpenURL

  14. Callis TE, Chen JF, Wang DZ: MicroRNAs in skeletal and cardiac muscle development.

    DNA Cell Biol 2007, 26(4):219-225. PubMed Abstract | Publisher Full Text OpenURL

  15. McCarthy JJ, Esser KA: MicroRNA-1 and microRNA-133a expression are decreased during skeletal muscle hypertrophy.

    J Appl Physiol 2007, 102(1):306-313. PubMed Abstract | Publisher Full Text OpenURL

  16. Chen JF, Mandel EM, Thomson JM, Wu Q, Callis TE, Hammond SM, Conlon FL, Wang DZ: The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation.

    Nat Genet 2006, 38(2):228-233. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Nguyen HT, Frasch M: MicroRNA in muscle differentiation: lessons from Drosophila and beyond.

    Curr Opin Genet Dev 2006, 16:533-539. PubMed Abstract | Publisher Full Text OpenURL

  18. Sweetman D, Goljanek K, Rathjen T, Oustanina S, Braun T, Dalmay T, Münsterberg A: Specific requirements of MRFs for the expression of muscle specific microRNAs, miR-1, miR-206 and miR-133.

    Dev Biol 2008, 321(2):491-499. PubMed Abstract | Publisher Full Text OpenURL

  19. Sweetman D, Rathjen T, Jefferson M, Wheeler G, Smith TG, Wheeler GN, Münsterberg A, Dalmay T: FGF-4 signaling is involved in mir-206 expression in developing somites of chicken embryos.

    Dev Dyn 2006, 235(8):2185-2191. PubMed Abstract | Publisher Full Text OpenURL

  20. Clop A, Marcq F, Takeda H, Pirottin D, Tordoir X, Bibe B, Bouix J, Caiment F, Elsen JM, Eychenne F, Larzul C, Laville E, Meish F, Milenkovic D, Tobin J, Charlier C, Georges M: A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep.

    Nat Genet 2006, 38(7):813-818. PubMed Abstract | Publisher Full Text OpenURL

  21. Griffiths-Jones S: The microRNA registry.

    NAR 2004, 32:D109-D111. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature.

    NAR 2006, 34:D140-D144. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Griffiths-Jones S, Saii HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics.

    NAR 2008, 36:D154-D158. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Bentwich I, Avniel A, Karov Y, Aharonov R, Gilad S, Barad O, Barzilai A, Einat P, Einav U, Meiri E, Sharon E, Spector Y, Bentwich Z: Identification of hundreds of conserved and nonconserved human microRNAs.

    Nat Genet 2005, 37(7):766-770. PubMed Abstract | Publisher Full Text OpenURL

  25. Wienholds E, Plasterk RH: MicroRNA function in animal development.

    FEBS Lett 2005, 579(26):5911-5922. PubMed Abstract | Publisher Full Text OpenURL

  26. Aboobaker AA, Tomancak P, Patel N, Rubin GM, Lai EC: Drosophila microRNAs exhibit diverse spatial expression patterns during embryonic development.

    Proc Natl Acad Sci 2005, 102(50):18017-18022. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. 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 2006, 29(1):35-43. PubMed Abstract | Publisher Full Text OpenURL

  28. Hafner M, Landgraf P, Ludwig J, Rice A, Ojo T, Lin C, Holoch D, Lim C, Tuschl T: Identification of microRNAs and other small regulatory RNAs using cDNA library sequencing.

    Methods 2008, 44(1):3-12. PubMed Abstract | Publisher Full Text OpenURL

  29. 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. PubMed Abstract | Publisher Full Text OpenURL

  30. Kim J, Cho IS, Hong JS, Choi YK, Kim H, Lee YS: Identification and characterization of new microRNAs from pig.

    Mamm Genome 2008, 19(7-8):570-580. PubMed Abstract | Publisher Full Text OpenURL

  31. Pette D, Staron RS: Mammalian skeletal muscle fiber type transitions.

    Int Rev Cytol 1997, 170:143-223. PubMed Abstract OpenURL

  32. Gunawan AM, Richert BT, Schinckel AP, Grant AL, Gerrard DE: Ractopamine induces differential gene expression in porcine skeletal muscles.

    J Anim Sci 2007, 85(9):2115-2124. PubMed Abstract | Publisher Full Text OpenURL

  33. Bai Q, McGillivray C, da Costa N, Dornan S, Evans G, Stear MJ, Chang KC: Development of a porcine skeletal muscle cDNA microarray: analysis of differential transcript expression in phenotypically distinct muscles.

    BMC Genomics 2003, 4(1):8. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Kim NK, Joh JH, Park HR, Kim OH, Park BY, Lee CS: Differential expression profiling of the proteomes and their mRNAs in porcine white and red skeletal muscles.

    Proteomics 2004, 4(11):3422-3428. PubMed Abstract | Publisher Full Text OpenURL

  35. Hamelin M, Sayd T, Chambon C, Bouix J, Bibe B, Milenkovic D, Leveziel H, Georges M, Clop A, Marinova P, Laville E: Differential expression of sarcoplasmic proteins in four heterogeneous ovine skeletal muscles.

    Proteomics 2007, 7(2):271-280. PubMed Abstract | Publisher Full Text OpenURL

  36. Neufer PD, Carey JO, Dohm GL: Transcriptional regulation of the gene for glucose transporter GLUT4 in skeletal muscle. Effects of diabetes and fasting.

    J Biol Chem 1993, 268(19):13824-13829. PubMed Abstract | Publisher Full Text OpenURL

  37. Takada S, H Mano: Profiling of microRNA expression by mRAP.

    Nat Protoc 2007, 2(12):3136-3145. PubMed Abstract | Publisher Full Text OpenURL

  38. Rao PK, Kumar RM, Farkhondeh M, Baskerville S, Lodish HF: Myogenic factors that regulate expression of muscle-specific microRNAs.

    Proc Natl Acad Sci USA 2006, 103(23):8721-8726. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Beermann DH, RG Cassens, GJ Hausman: A second look at fiber type differentiation in porcine skeletal muscle.

    J Anim Sci 1978, 46(1):125-132. PubMed Abstract | Publisher Full Text OpenURL

  40. Cimmino A, Calin GA, Fabbri M, Iorio MV, Ferracin M, Shimizu M, Wojcik SE, Aqeilan RI, Zupo S, Dono M, Rassenti L, Alder H, Volinia S, Liu CG, Kipps TJ, Negrini M, Croce CM: miR-15 and miR-16 induce apoptosis by targeting BCL2.

    Proc Natl Acad Sci 2005, 102(39):13944-13949. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. Linsley PS, Schelter J, Burchard J, Kibukawa M, Martin MM, Bartz SR, Johnson JM, Cummins JM, Raymond CK, Dai H, Chau N, Cleary M, Jackson AL, Carleton M, Lim L: Transcripts targeted by the microRNA-16 family cooperatively regulate cell cycle progression.

    Mol Cell Biol 2007, 27(6):2240-2252. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Cheng AM, Byrom MW, Shelton J, Ford LP: Antisense inhibition of human miRNAs and indications for an involvement of miRNA in cell growth and apoptosis.

    Nucleic Acids Res 2005, 33(4):1290-1297. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. He L, He X, Lowe SW, Hannon GJ: microRNAs join the p53 network – another piece in the tumour-suppression puzzle.

    Nat Rev Cancer 2007, 7(11):819-822. PubMed Abstract | Publisher Full Text OpenURL

  44. Gottardo F, Liu CG, Ferracin M, Calin GA, Fassan M, Bassi P, Sevignani C, Byrne D, Negrini M, Pagano F, Gomella LG, Croce CM, Baffa R: Micro-RNA profiling in kidney and bladder cancers.

    Urol Oncol 2007, 25(5):387-392. PubMed Abstract | Publisher Full Text OpenURL

  45. Mott JL, Kobayashi S, Bronk SF, Gores GJ: mir-29 regulates Mcl-1 protein expression and apoptosis.

    Oncogene 2007, 26(42):6133-6140. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Pekarsky Y, Santanam U, Cimmino A, Palamarchuk A, Efanov A, Maximov V, Volinia S, Alder H, Liu CG, Rassenti L, Calin GA, Hagan JP, Kipps T, Croce CM: Tcl1 expression in chronic lymphocytic leukemia is regulated by miR-29 and miR-181.

    Cancer Res 2006, 66(24):11590-11593. PubMed Abstract | Publisher Full Text OpenURL

  47. Rosa A, Ballarino M, Sorrentino A, Sthandier O, De Angelis FG, Marchioni M, Masella B, Guarini A, Fatica A, Peschle C, Bozzoni I: The interplay between the master transcription factor PU.1 and miR-424 regulates human monocyte/macrophage differentiation.

    Proc Natl Acad Sci USA 2007, 104(50):19849-19854. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Tavazoie SF, Alarcón C, Oskarsson T, Padua D, Wang Q, Bos PD, Gerald WL, Massagué J: Endogenous human microRNAs that suppress breast cancer metastasis.

    Nature 2008, 451(7175):147-152. PubMed Abstract | Publisher Full Text OpenURL

  49. Lewis BP, Shih I, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets.

    Cell 2003, 115:787-798. PubMed Abstract | Publisher Full Text OpenURL

  50. Li X, Zhang YZ: Computational detection of microRNAs targeting transcription factor genes in Arabidopsis thaliana.

    Comput Biol Chem 2005, 29(5):360-367. PubMed Abstract | Publisher Full Text OpenURL

  51. Rehmsmeier M, Steffen P, Höchsmann M, Giegerich R: Fast and effective prediction of microRNA/target duplexes.

    RNA 2004, 10:1507-1517. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Maglott D, Ostell J, Pruitt KD, Tatusova T: Entrez Gene: gene-centered information at NCBI.

    Nucleic Acids Res 2005, 33:D54-D58. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  53. Winder SJ, Walsh MP, Vasulka C, Johnson JD: Calponin-calmodulin interaction: properties and effects on smooth and skeletal muscle actin binding and actomyosin ATPases.

    Biochemistry 1993, 32(48):13327-13333. PubMed Abstract OpenURL

  54. Tajbakhsh S, Rocancourt D, Cossu G, Buckingham M: Redefining the genetic hierarchies controlling skeletal myogenesis: Pax-3 and Myf-5 act upstream of MyoD.

    Cell 1997, 89(1):127-138. PubMed Abstract | Publisher Full Text OpenURL

  55. Ueda H, Ohno S, Kobayashi T: Myotonic dystrophy and myotonic dystrophy protein kinase.

    Prog Histochem Cytochem 2000, 35(3):187-251. PubMed Abstract OpenURL

  56. Doumit ME, Merkel RA: Conditions for isolation and culture of porcine myogenic satellite cells.

    Tissue Cell 1992, 24(2):253-262. PubMed Abstract | Publisher Full Text OpenURL

  57. Foxcroft GR, Dixon WT, Novak S, Putman CT, Town SC, Vinsky MD: The biological basis for prenatal programming of postnatal performance in pigs.

    J Anim Sci 2006, 84(Suppl):E105-E112. PubMed Abstract | Publisher Full Text OpenURL

  58. Lu C, Tej SS, Luo S, Haudenschild CD, Meyers BC, Green PJ: Elucidation of the small RNA component of the transcriptome.

    Science 2005, 309:1567-1569. PubMed Abstract | Publisher Full Text OpenURL

  59. Ewing B, Hillier L, Wendl MC, Green P: Base calling of automated sequencer traces using phred. I. Accuracy assessment.

    Genome Res 1998, 8:175-185. PubMed Abstract | Publisher Full Text OpenURL