当前位置:首页 >> 畜牧兽医 >>

BMC Genomics 2009 Ralph MicroRNA transcriptome profiles during swine skeletal muscle development


BMC Genomics
Research article

BioMed Central

Open Access

MicroRNA transcriptome profiles during swine skeletal muscle development
Tara G McDaneld1, Timothy PL Smith*1, Matthew E Doumit2, Jeremy R Miles1, Luiz L Coutinho3, Tad S Sonstegard4, Lakshmi K Matukumalli5, Dan J Nonneman1 and Ralph T Wiedmann1
Address: 1USDA/ARS Meat Animal Research Center, Clay Center, NE, USA, 2Department of Animal Science, Michigan State University, East Lansing, MI, USA, 3University of Sao Paulo-ESALQ, Piracicaba, SP, Brazil, 4USDA/ARS Beltsville Area Research Center, Beltsville, MD, USA and 5Bioinformatics and Computational Biology, George Mason University, Manassas, VA, USA Email: Tara G McDaneld - tara.mcdaneld@ars.usda.gov; Timothy PL Smith* - tim.smith@ars.usda.gov; Matthew E Doumit - mdoumit@uidaho.edu; Jeremy R Miles - jeremy.miles@ars.usda.gov; Luiz L Coutinho - llcoutin@esalq.usp.br; Tad S Sonstegard - tads@anri.barc.usda.gov; Lakshmi K Matukumalli - lakshmi.matukumalli@ars.usda.gov; Dan J Nonneman - dan.nonneman@ars.usda.gov; Ralph T Wiedmann - ralph.wiedmann@ars.usda.gov * Corresponding author

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

Received: 16 October 2008 Accepted: 10 February 2009

This article is available from: http://www.biomedcentral.com/1471-2164/10/77 ? 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. Musclespecific 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 dayold 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 miR368, 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.

Page 1 of 11
(page number not for citation purposes)

BMC Genomics 2009, 10:77

http://www.biomedcentral.com/1471-2164/10/77

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

Page 2 of 11
(page number not for citation purposes)

BMC Genomics 2009, 10:77

http://www.biomedcentral.com/1471-2164/10/77

Figure 1 Saturation plot of microRNA libraries 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.

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 esti-

mates 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. 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

Page 3 of 11
(page number not for citation purposes)

BMC Genomics 2009, 10:77

http://www.biomedcentral.com/1471-2164/10/77

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

A
103

log (# observed per 1000)

4th and 5th passage 6th passage, #1 6th passage, #2
102

101

PN1

Let-7

miR-29

miR-34

miR-93

miR-450

miR503

miR-26

miR-15

miR-16

miR-18

miR-20

miR-21

miR-22

miR-24

miR-27

miR-103

miR-148

miR-185

miR-424

miR-543

miR-106

miR-130

miR-143

miR-151*

miR-199*

B
103 Adult BF #1

C
103 d90 F BF #1

log (# observed per 1000)

Adult BF #2 102

log (# observed per 1000)

miR-423-5p

102

101

101

miR-133

miR-143

miR-151

miR-206

miR-22

miR-24

miR-26

miR-27

Let-7

miR-1

miR-423-5p

Let-7

miR-199

miR-206

Figure 2 Variation between microRNA libraries 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.

(page number not for citation purposes)

miR-423-5p

Page 4 of 11

miR-432

miR-1

0

0

miR-744

d90 F BF #2

PN2

0

BMC Genomics 2009, 10:77

http://www.biomedcentral.com/1471-2164/10/77

result of this analysis, we concluded that changes in abundance are likely to be significant if they are greater than 6fold. 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. 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. 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, miR206 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. 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-

Page 5 of 11
(page number not for citation purposes)

BMC Genomics 2009, 10:77

http://www.biomedcentral.com/1471-2164/10/77

Figure 3 MicroRNA transcriptome profiles 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.

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.

Page 6 of 11
(page number not for citation purposes)

BMC Genomics 2009, 10:77

http://www.biomedcentral.com/1471-2164/10/77

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 miR1 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 miR381, 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: miR206, 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.

Table 1: miR-206 sequence variation

Sequence UGGAAUGUAAGGAAGUGUGUGA UGGAAUGUAAGGAAGUGUGUGAA UGGAAUGUAAGGAAGUGUGUGU UGGAAUGUAAGGAAGUGUGU UGGAAUGUAAGGAAGUGUGUGAU UGGAAUGUAAGGAAGUGUGUGA UGGAAUGUAAGGAAGUGUGUG UGGAAUGUAAGGAAGUGUGU UGGAAUGUAAGGAAGUGUGUGG CGGAAUGUAAGGAAGUGUGUGA

Quantity observed 2286 79 78 41 34 15 13 9 8 7

Mir-206 differed in length or sequence at the 5' and 3' ends. MiR-206 sequences of this type were clustered into a single miR sequence identified as the predominant sequence. Mir-206 sequences with the greatest observations are listed.

sequence of the cluster as the defining sequence for comparisons to known miR. 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.
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 (miR206, 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

In addition to identification of mismatches in the miR206 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

Page 7 of 11
(page number not for citation purposes)

BMC Genomics 2009, 10:77

http://www.biomedcentral.com/1471-2164/10/77

Table 2: Novel porcine miR

Temporary identifiers PN1 PN2 PN3 PN4 PN5 PN6 PN7 PN8 PN9 PN10 PN11 PN12

Sequence CCGCAGGUGCGGCCACUUGUUU GUGUUGGUGUGCACUUAUUU CGAACCGAAUCCCUCACUAAA AGGGGAGUGGUGGGGGGAG CGAACCGAACUCCUCACUAAA AGGGUUGGGCGGAGGCUUUCC CCACGAGGAGGAGACGCAGUG UGGCACAGGGUCCAGCUGUCGGC GGGGUGGGGGUCUGGGGGGUGU GAGAGAUCAGAGGCGCAGAGU GUGUGGGACGGUGGGGUGGGUU GUCGGGGAGGUUCCAGCUCUCAUUU

Twelve novel miR were detected that were not closely related to previously reported miR in the database. The miR were given temporary identifiers PN (porcine new) 1 to PN12.

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, miR27, 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.

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 accord-

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

Page 8 of 11
(page number not for citation purposes)

BMC Genomics 2009, 10:77

http://www.biomedcentral.com/1471-2164/10/77

ance 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.

Additional material 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. Click here for file [http://www.biomedcentral.com/content/supplementary/14712164-10-77-S1.docx]

Additional file 2
miR-206 sequence variation. The data provided represent the variation of miR-206 sequences identified in the current study. Click here for file [http://www.biomedcentral.com/content/supplementary/14712164-10-77-S2.docx]

Additional file 3
Predicted targets of up-regulated miR. The data provided represent the predicted gene targets for the up-regulated miR. Click here for file [http://www.biomedcentral.com/content/supplementary/14712164-10-77-S3.docx]

Page 9 of 11
(page number not for citation purposes)

BMC Genomics 2009, 10:77

http://www.biomedcentral.com/1471-2164/10/77

20.

Additional file 4
Predicted targets of up-regulated miR. The data provided represent the predicted gene targets for the down-regulated miR. Click here for file [http://www.biomedcentral.com/content/supplementary/14712164-10-77-S4.docx]

21. 22. 23.

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

24.

References
1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 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. 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. 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. Lee RC, Ambros V: An extensive class of small RNAs in Caenorhabditis elegans. Science 2001, 294(5543):862-864. Hutvagner G: Small RNA asymmetry in RNAi: function in RISC assembly and gene regulation. FEBS Lett 2005, 579(26):5850-5857. Lai EC: miRNAs: whys and wherefores of miRNA-mediated regulation. Curr Biol 2005, 15(12):R458-R460. Chan SP, Slack FJ: microRNA-mediated silencing inside P-bodies. RNA Biol 2006, 3(3):97-100. Carthew RW: Gene regulation by microRNAs. Curr Opin Genet Dev 2006, 16(2):203-208. Garzon R, Fabbri M, Cimmino A, Calin GA, Croce CM: MicroRNA expression and function in cancer. Trends Mol Med 2006, 12(12):580-587. Brennecke J, Stark A, Cohen SM: Not miR-ly muscular: microRNAs and muscle development. Genes Dev 2005, 19(19):2261-2264. Kajimoto K, Naraba H, Iwai N: MicroRNA and 3T3-L1 pre-adipocyte differentiation. RNA 2006, 12(9):1626-1632. Anderson C, Catoe H, Werner R: MIR-206 regulates connexin43 expression during skeletal muscle development. Nucleic Acids Res 2006, 34(20):5863-5871. 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. Callis TE, Chen JF, Wang DZ: MicroRNAs in skeletal and cardiac muscle development. DNA Cell Biol 2007, 26(4):219-225. McCarthy JJ, Esser KA: MicroRNA-1 and microRNA-133a expression are decreased during skeletal muscle hypertrophy. J Appl Physiol 2007, 102(1):306-313. Chen JF, Mandel EM, Thomson JM, Wu Q, Callis TE, Hammond SM, Conlon FL, Wang DZ: The role of microRNA-1 and microRNA133 in skeletal muscle proliferation and differentiation. Nat Genet 2006, 38(2):228-233. Nguyen HT, Frasch M: MicroRNA in muscle differentiation: lessons from Drosophila and beyond. Curr Opin Genet Dev 2006, 16:533-539. 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 miR133. Dev Biol 2008, 321(2):491-499. 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.

25. 26.

27.

28.

29. 30. 31. 32. 33.

34.

35.

36.

14. 15. 16.

37. 38. 39. 40.

17. 18.

19.

41.

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. Griffiths-Jones S: The microRNA registry. NAR 2004, 32:D109-D111. Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. NAR 2006, 34:D140-D144. Griffiths-Jones S, Saii HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. NAR 2008, 36:D154-D158. 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. Wienholds E, Plasterk RH: MicroRNA function in animal development. FEBS Lett 2005, 579(26):5911-5922. 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. 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. 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. 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. 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. Pette D, Staron RS: Mammalian skeletal muscle fiber type transitions. Int Rev Cytol 1997, 170:143-223. 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. 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. 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. 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. 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. Takada S, H Mano: Profiling of microRNA expression by mRAP. Nat Protoc 2007, 2(12):3136-3145. 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. 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. 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. 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.

Page 10 of 11
(page number not for citation purposes)

BMC Genomics 2009, 10:77

http://www.biomedcentral.com/1471-2164/10/77

42.

43. 44.

45. 46.

47.

48.

49. 50. 51. 52. 53.

54. 55. 56. 57. 58. 59.

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. 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. 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. Mott JL, Kobayashi S, Bronk SF, Gores GJ: mir-29 regulates Mcl-1 protein expression and apoptosis. Oncogene 2007, 26(42):6133-6140. 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. 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. 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. Lewis BP, Shih I, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell 2003, 115:787-798. Li X, Zhang YZ: Computational detection of microRNAs targeting transcription factor genes in Arabidopsis thaliana. Comput Biol Chem 2005, 29(5):360-367. Rehmsmeier M, Steffen P, H?chsmann M, Giegerich R: Fast and effective prediction of microRNA/target duplexes. RNA 2004, 10:1507-1517. Maglott D, Ostell J, Pruitt KD, Tatusova T: Entrez Gene: geneNucleic Acids Res 2005, centered information at NCBI. 33:D54-D58. 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. Tajbakhsh S, Rocancourt D, Cossu G, Buckingham M: Redefining the genetic hierarchies controlling skeletal myogenesis: Pax3 and Myf-5 act upstream of MyoD. Cell 1997, 89(1):127-138. Ueda H, Ohno S, Kobayashi T: Myotonic dystrophy and myotonic dystrophy protein kinase. Prog Histochem Cytochem 2000, 35(3):187-251. Doumit ME, Merkel RA: Conditions for isolation and culture of porcine myogenic satellite cells. Tissue Cell 1992, 24(2):253-262. 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. 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. 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.

Publish with Bio Med Central and every scientist can read your work free of charge
"BioMed Central will be the most significant development for disseminating the results of biomedical researc h in our lifetime."
Sir Paul Nurse, Cancer Research UK

Your research papers will be:
available free of charge to the entire biomedical community peer reviewed and published immediately upon acceptance cited in PubMed and archived on PubMed Central yours — you keep the copyright
Submit your manuscript here:
http://www.biomedcentral.com/info/publishing_adv.asp

BioMedcentral

Page 11 of 11
(page number not for citation purposes)


赞助商链接
相关文章:
更多相关标签: