Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Identification and Profiling of MicroRNAs from Skeletal Muscle of the Common Carp

  • Xuechun Yan ,

    Contributed equally to this work with: Xuechun Yan, Lei Ding

    Affiliation Key Laboratory of Aquatic Biotechnology, Heilongjiang River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Harbin, China

  • Lei Ding ,

    Contributed equally to this work with: Xuechun Yan, Lei Ding

    Affiliation College of Life Sciences, Northeast Forestry University, Harbin, China

  • Yunchao Li,

    Affiliation College of Life Sciences, Northeast Forestry University, Harbin, China

  • Xiaofeng Zhang,

    Affiliation Key Laboratory of Aquatic Biotechnology, Heilongjiang River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Harbin, China

  • Yang Liang,

    Affiliation College of Life Sciences, Northeast Forestry University, Harbin, China

  • Xiaowen Sun ,

    sxw2002@163.com (XS); chunboteng@nefu.edu.cn (CBT)

    Affiliation Key Laboratory of Aquatic Biotechnology, Heilongjiang River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Harbin, China

  • Chun-Bo Teng

    sxw2002@163.com (XS); chunboteng@nefu.edu.cn (CBT)

    Affiliation College of Life Sciences, Northeast Forestry University, Harbin, China

Abstract

The common carp is one of the most important cultivated species in the world of freshwater aquaculture. The cultivation of this species is particularly productive due to its high skeletal muscle mass; however, the molecular mechanisms of skeletal muscle development in the common carp remain unknown. It has been shown that a class of non-coding ∼22 nucleotide RNAs called microRNAs (miRNAs) play important roles in vertebrate development. They regulate gene expression through sequence-specific interactions with the 3′ untranslated regions (UTRs) of target mRNAs and thereby cause translational repression or mRNA destabilization. Intriguingly, the role of miRNAs in the skeletal muscle development of the common carp remains unknown. In this study, a small-RNA cDNA library was constructed from the skeletal muscle of the common carp, and Solexa sequencing technology was used to perform high throughput sequencing of the library. Subsequent bioinformatics analysis identified 188 conserved miRNAs and 7 novel miRNAs in the carp skeletal muscle. The miRNA expression profiling showed that, miR-1, miR-133a-3p, and miR-206 were specifically expressed in muscle-containing organs, and that miR-1, miR-21, miR-26a, miR-27a, miR-133a-3p, miR-206, miR-214 and miR-222 were differentially expressed in the process of skeletal muscle development of the common carp. This study provides a first identification and profiling of miRNAs related to the muscle biology of the common carp. Their identification could provide clues leading towards a better understanding of the molecular mechanisms of carp skeletal muscle development.

Introduction

The common carp (Cyprinus Carpio L.) occupies a prominent position in the world of freshwater aquaculture and serves as a major source of animal protein for millions of people [1]. In spite of its aquacultural importance and a production increase in carp breeding over the last decade, the current productive species are mainly selected through traditional breeding methods, including the use of mutual hybrids and backcrosses to domesticated ancestors [2]. Molecular techniques are still not widely used in the breeding of the common carp due to a lack of genetic and genomic information. In recent years, several carp traits, including cold tolerance, body weight and plasma glucose contents, have been studied through development of molecular markers and construction of genetic linkage maps [3][5]. Although it is an important economic trait, the skeletal muscle mass of the common carp has not been well defined by molecular markers or control mechanisms to date.

Several lines of evidence in mouse and zebrafish models have suggested that the development of vertebrate skeletal muscle is orchestrated by evolutionarily conserved gene expression networks in which transcription factors play important roles in regulating muscle growth and differentiation [6]. An essential myogenic transcription factor is the MADS box transcription factor myocyte enhancer factor-2 (MEF2), which activates the myogenic differentiation program in conjunction with the basic-helix-loop-helix (bHLH) transcription factors MyoD and myogenin in skeletal muscle [7]. Recent studies have shown that, in addition to transcriptional factors involved in muscle proliferation and differentiation, a set of microRNAs (miRNAs) also play important roles in skeletal muscle development in vertebrate animals [8], [9].

miRNAs are a class of non-coding RNAs of approximately 22 nucleotides in length, which are processed from stem-loop RNA precursors (pre-miRNAs, ∼65 nucleotides) by the ribonuclease III (RNase III) enzyme Dicer [10]. The pre-miRNAs are excised from the longer primary transcripts (pri-miRNAs) in the nucleus by the RNAse III enzyme Drosha [11]. The mature miRNAs exert their functions by incorporating one of the duplex strands into the RNA-induced silencing complexes (RISCs) and then binding the 3′ untranslated region (3′UTR) of the target mRNAs to degrade or repress their translation [12]. Many studies have demonstrated that miRNAs exist in a wide range of invertebrates and vertebrates, and that they are involved in the regulation of animal development, physiological functions and pathological processes [13][16].

miRNAs can be identified through manual cloning and sequencing, miRNA array screening, and bioinformatic approaches [17], [18]. Recently, high throughput sequencing technology has been widely used to facilitate the identification and detection of miRNAs in multiple species [19][23]. For this study, we constructed a small-RNA cDNA library from the skeletal muscle of the common carp. Through high throughput sequencing of the small RNA library and subsequent bioinformatic analysis, miRNAs in skeletal muscle tissue of the common carp were identified. Finally, using expression profiling analysis of the miRNAs from different tissues and multiple development stages, miRNAs likely to play a role in carp muscle development were identified.

Results

Construction of cDNA library for sequencing and small RNA discovery

The cDNA library of small RNAs was constructed with pooled total RNAs from skeletal muscle tissues collected from three different 1-year-old common carps. Through high throughput Solexa sequencing, 21335605 total reads were obtained; 93.49% of the reads were 21–23 nucleotides in length (Figure 1A). After removal of the 5′ and 3′ adapters, pollution reads and reads smaller than 18 nucleotides, 20842912 high-quality clean reads were extracted, and 25966 of the resulting unique small RNAs were annotated as either rRNA, snRNA, snoRNA, or tRNA using BLAST searches against the NCBI Genebank and Rfam database (Figure 1B, 1C) and then removed from the following analysis. The remainder of the unannotated small RNAs was mapped to the zebrafish genome, which is the only completely sequenced genome of the Cyprinidae fish family, to perform distribution analysis using SOAP (http://soap.genomics.org.cn). 35815 unique sequences (16965286 reads) were perfectly mapped to the zebrafish genome.

thumbnail
Figure 1. Anaylsis of small RNAs derived from the skeletal muscle of the common carp.

A: Length distribution of sequenced small RNAs; B: The clean reads were blasted against the GenBank database to annotate rRNA, RNA, snRNA, or snoRNA; C: The clean reads were blasted against the Rfam database noncoding RNA to annotate rRNA, RNA, snRNA, or snoRNA.

https://doi.org/10.1371/journal.pone.0030925.g001

Identification of miRNAs in the skeletal muscle of common carp

To identify conserved miRNAs in the skeletal muscle tissue of the common carp, the small RNAs with a length of 18–23 nucleotides were searched against the latest miRBase release version 17.0 (miRBase V17.0). 2422 unique sequences (16522982 reads) were annotated as miRNA candidates, while the rest (3916545 reads) were unannotated. The 2422 miRNA candidates were then clustered into 188 categories according to sequence similarity. Each category included multiple homologs, which differed in sequence length by only 1–5 nucleotides. Such homologous sequences with different lengths are thought to be variants produced by various biochemical modifications and by imprecise processing of primary or precursor miRNAs by Drosha and Dicer enzymes. The sequence with the highest read in a category was considered to be the representative sequence for the following analysis.

Apart from the characteristic length (∼22 nucleotides) and conservation, an important criteria for annotating a cloned sequence as a miRNA is the presence of a compact pre-miRNA fold-back structure [24]. Since we are sequencing the genome of the common carp ourselves and aligning the sequence to the chromosomes, we selected 20 candidate miRNA sequences to compare to the sequenced ESTs and genome sequence of the common carp and predict possible precursor sequences. Subsequently, we predicted the miRNA precursor structure using mfold software and found that all the 20 sequences could form miRNA precursors with stem-loop structure (Figure 2, Figure S1), indicating that annotating the cloned sequences as miRNAs through conservation analysis is quite reliable. The 188 conserved sequences were annotated as common carp miRNAs (Table S1).

thumbnail
Figure 2. Prediction of the fold-back structure of 10 carp miRNA precursors.

The precursor sequences were obtained by sequence alignment with the sequences of the common carp genome and ESTs. The mature miRNA sequences in the precursors are indicated in red.

https://doi.org/10.1371/journal.pone.0030925.g002

The unannotated small RNAs that could be mapped to the zebrafish genome were subjected to novel miRNA prediction analysis of their secondary structure, the Dicer enzyme cleavage site and the minimum free energy using Mireap software (Beijing Genomics Institute). 15 specific reads of ∼22 nucleotides in length were considered miRNA candidates. These sequences were further searched against the ESTs and genome sequence of the common carp in an attempt to predict possible precursor sequences. 7 miRNA precursors with stem-loop structure were found using mfold software (Figure S2); thus 7 novel and specific miRNAs were identified from the skeletal muscle tissue of the common carp (Table 1).

Conservation analysis of the identified common carp miRNAs

To date, although the miRBase V17.0 has included miRNAs from 84 different animal species, miRNA members from non-model animals have not been fully identified. To analyze the conservation of common carp miRNAs, we compared them to those identified in model animals, including worm (Caenorhabditis elegans), fly (Drosophila melanogaster), frog (Xenopus tropocalis), zebrafish (Danio rerio), chicken (Gallus gallus), mouse (Mus musculus) and human (Homo sapiens). Only let-7, miR-1, and miR-139 are highly conserved from worm to human, but 96 common carp miRNAs are conserved in all vertebrate model animals (Table S1), although there exist base edition and length change of sequences across the different animal species. Moreover, 31 common carp miRNAs are only identified in teleosts (Table 2). The comparison was performed between the miRNAs discovered in the common carp with those identified in the bighead carp and silver carp [23], and the result shows that 80% common carp miRNAs are discovered in the bighead carp and silver carp (Table S3).

thumbnail
Table 2. Common carp miRNAs only conserved in the teleosts.

https://doi.org/10.1371/journal.pone.0030925.t002

Let-7, one of the most conserved miRNAs, has evolved into a miRNA family that includes 8 members in the common carp. Of these, 5 members possess base or length edition between different vertebrate model animals but no changes in the seed sequence (Figure 3). Apart from let-7, 32 other miRNA families were found to be included in the 188 conserved common carp miRNAs. Of these, only cca-miR-16c, cca-miR-193b, cca-miR-200c, and cca-miR-301c have differences in the seed sequence from their respective family members (Table S2).

thumbnail
Figure 3. Comparation of the member sequences of the let-7 family in model vertebrate animals.

Bases that differ between model vertebrate animals are indicated in red.

https://doi.org/10.1371/journal.pone.0030925.g003

Through literature mining, we found that 27 miRNAs are known to be involved in animal muscle development [8], [25][47], among them 16 miRNAs were found in the common carp (Table 3). Conservation analysis showed that the 16 miRNAs were highly conserved in multiple animal species (Table S1).

thumbnail
Table 3. The miRNAs both identified in the common carp and demonstrated to be muscle-related in several species.

https://doi.org/10.1371/journal.pone.0030925.t003

Validation and expression analysis of the identified common carp miRNAs

For validation and identification of muscle-related miRNAs in the common carp, RT-PCR analysis of miRNA expression was performed in skeletal muscle, heart, gut, liver, kidney, skin, gill, eye, and brain of the common carp. We first randomly picked out several high-read miRNAs, including let-7a, miR-10a-5p, miR-143, and miR-146a. The results showed that all 4 miRNAs could be amplified from the skeletal muscle, heart, gut, liver, kidney, skin, gill, eye, and brain, and that the expression levels of these miRNAs were relatively consistent across all 9 tissue types (Fig. 4A). These results indicate that let-7a, miR-10a-5p, miR-143, and miR-146a are ubiquitously expressed in multiple tissues and organs.

thumbnail
Figure 4. The expression of miRNAs in tissues and organs of the common carp were detected by RT-PCR.

A: Several high-read miRNAs were detected in skeletal muscle (m), heart (h), gut (gu), liver (l), kidney (k), skin (s), gill (gi), eye (e), and brain (b) of the common carp; B: Eleven muscle-related miRNAs were detected in skeletal muscle (m), heart (h), gut (gu), liver (l), kidney (k), skin (s), gill (gi), eye (e), and brain (b) of the common carp. Expression of miRNAs was normalized to U6.

https://doi.org/10.1371/journal.pone.0030925.g004

We then selected miR-1, miR-21, miR-23a, miR-24, miR-26a, miR-27a, miR-133a-3p, miR-181a-5p, miR-206, miR-214, and miR-222, which have been reported to play roles in the development of vertebrate skeletal muscle (Table 3). Of these, miR-1, miR-133a-3p, and miR-206, were found to be specifically expressed in heart, skeletal muscle, gut, eye, and skin (Fig. 4B), whereas miR-21, miR-26a, miR-181a-5p, and miR-222 were detected in all the samples with little expression level change. Importantly, miR-23a, miR-24, miR-27a, and miR-214 were highly expressed in muscle-related tissue or organs (skeletal muscle, heart, gut), although they could also be detected in non-muscle organs (Fig. 4B). These muscle-specific and highly expressed miRNAs may therefore play roles in the skeletal muscle development of the common carp.

Profiling of the identified common carp miRNAs at different developmental stages

To further explore the miRNAs involved in the skeletal muscle development of the common carp, we performed quantitative analysis of the miRNAs in skeletal muscle tissues of 30 day post-hatching (dph) larva, 1-year-old common carps, and 2-year-old common carps. The results showed that the expression of miR-1, miR-133a-3p, miR-206, and miR-21 was increased in the skeletal muscles from 30 day post-hatching to 2-year-old (Figure 5). Among four miRNAs, miR-1 and miR-133a-3p were expressed at significantly higher levels in the 2-year skeletal muscles than 1-year skeletal muscles (Figure 5). MiR-26a expression was not changed in the 30 dph and 1-year skeletal muscles, but increased in the 2-year skeletal muscles (Figure 5). These results suggest that miR-1, miR-133a-3p, miR-21, miR-26, and miR-206 play roles in the growth and functional maintenance of carp skeletal muscles. miR-27a, miR-214, and miR-222 are expressed at relatively high levels in the 30 dph skeletal muscles, and then decreased in the 1-year and 2-year skeletal muscles (Figure 5). These results indicate that miR-27a, miR-214, and miR-222 are related to the process of skeletal muscle proliferation and differentiation in the common carp.

thumbnail
Figure 5. The expression of muscle-related miRNAs during skeletal muscle development of the common carp was detected by qRT-PCR.

The profile of miRNAs was measured at the skeletal muscle tissues of 30 day post-hatching (dph) larva, 1-year old and 2-year old common carps. The level of miRNA expression in the skeletal muscles of 30 dph larva was defined as 1. 30 d, 30 dph skeletal muscles; 1 y, 1-year skeletal muscles; 2 y, 2-year skeletal muscles. At least 5 animals were used for each timepoint.

https://doi.org/10.1371/journal.pone.0030925.g005

The expression of miRNAs was also analyzed at different developmental stages, including fertilized oocytes, 24 hour post-fertilization (hpf) embryos, 48 hpf embryos, 1 dph larva, 7 dph larva, 15 dph larva, and 30 dph larva. The results showed that miR-1, miR-133a-3p, and miR-206 were increased expression 48 hour post fertilization, and the other miRNAs were increased expression 7 day post hatching (Figure S3). Of note, half of detected miRNAs (miR-21, miR-26a, miR-27a, and miR-214) were expressed at a higher level in the fertilized oocytes than in embryos after fertilization (Figure S3); this is likely caused by the existence of maternal miRNAs or dicer.

Discussion

Vertebrate skeletal muscle develops through the proliferation and terminal differentiation of somatic cells committed to myogenic lineage [48]. Multiple transcriptional factors and miRNAs are known to be involved in skeletal myogenesis in mouse and zebrafish models [8], [9]. However, the molecular mechanisms of skeletal muscle development in the common carp are not well understood. In this study, we focused on the identification and profiling of miRNAs involved in the development of carp skeletal muscle. Using Solexa sequencing technology and bioinformatic analysis, 188 conserved miRNAs and 7 novel miRNAs were identified from carp skeletal muscle. Among the 188 conserved miRNAs, 97 are conserved in all vertebrate model animals. Moreover, the closer is the evolutionary link between two species, the more conserved their miRNA sequences are. Of the 27 miRNAs that have previously been reported to be involved in muscle development, 16 were found in the common carp and multiple animal species. These highly conserved miRNAs might be expected to perform similar functions in carp muscle development as they do in model animal species.

Profiling of miRNAs expression in the common carp indicate that miR-1, miR-133a-3p, and miR-206 are specifically expressed in the muscle-rich tissues and organs (skeletal muscle, heart, gut) in the common carps, and that their expression increases with the progress of common carp muscle development, ultimately reaching their highest levels in the skeletal tissues of mature carp. These results are consistent with previous findings in mouse, zebrafish and pig [29], [30], [32], [49]. It has been reported that a loss of Dicer function or a down-regulation of miR-1 and miR-133 in zebrafish alters muscle gene expression and disrupts actin organization during sarcomere assembly [32]. And studies in vitro have shown that miR-1, miR-133, and miR-206 can target multiple muscle-development-related genes, including Fstl1, Utrn [50], HDAC4 [29], Pax7 [30], c-Met [31], [37], SRF [29]. Although the carp-specific target genes of miR-1, miR-133, and miR-206 are not known, their specific expression pattern and high conservation indicate that they are also likely to play a role in the development of carp skeletal muscle.

Other muscle-related miRNAs in vertebrates, including miR-21, miR-214, miR-26, miR-27, miR-221, and miR-222, have also been identified in the common carp. miR-21 is a key determinant in the IL-11/Stat3 anti-apoptotic signaling pathway involved in the preconditioning of skeletal myoblasts [25], miR-26a targets Ezh2 during myogenesis [44], miR-27 regulates behaviors of muscle stem cells by repression of Pax3 expression [45], miR-214 promotes myogenic differentiation by facilitating exit from mitosis via down-regulation of proto-oncogene N-ras [38], and miR-221 and miR-222 modulate differentiation and maturation of skeletal muscle cells by targeting p27 [41]. However, whether these miRNAs influence carp muscle development by targeting similar genes requires further investigation.

Materials and Methods

Animals and Tissue collection

All animal experiments were performed using wild common carps bred in the Heilongjiang Fishery Institute in accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of China. The study was approved by the Committee on the Ethics of Animal Experiments of the Heilongjiang Fishery Institute (Permit Number: 11-1970).

Male and female adult common carp were separately kept in tanks with oxygen-saturated water at 23–25°C. Pituitary extract was used to induce and synchronize ovulation and spermiation. Artificial propagation was performed in hatcheries, and the incubation of early embryos was carried out in hatchery jars. Hatched fry were kept in a large glass tanks at room temperature. The experimental materials include fertilized oocytes, 24 hour post-fertilization (hpf) embryos, 48 hpf embryos, 1 day post-hatching (dph) larva, 7 dph larva, 15 dph larva, 30 dph larva, skeletal muscle tissues from 30 dph larva, 1-year-old and 2-year-old common carp. Tissues including skeletal muscle, skin, heart, liver, kidney, gill, gut, eye, and brain were collected by dissection from 1-year and 2-year-old common carp. The samples were snap-frozen in liquid nitrogen and stored at −80°C.

Small RNA isolation and cDNA library construction

Skeletal muscle tissues used for the generation of the small RNA library were obtained from 3 individual 1-year-old fish. Total RNA was isolated from skeletal muscle using TRIzol reagent (Invitrogen). Subsequently, the population of recovered small-RNAs, ranging in size from 18 to 30 nucleotides was purified from 15% polyacrylamide gel, and these small RNAs were ligated with the 3′ and 5′ adapter ligation. Reverse transcription reactions were performed using the RT primer, and PCR reactions were performed using the forward and reverse primers. The PCR product was purified via phenol/chloroform extraction and ethanol precipitation and was delivered to the Beijing Genomics Institute to be sequenced using Solexa technologies.

Sequence analysis and identification of miRNAs

Initial reads obtained from Solexa sequencing were processed by removing poor quality reads, 5′ adapter pollution reads, reads without 3′ adapter, reads without insert fragment, reads containing poly(A) stretches, and reads less than 18 nt. The clean reads were blasted against the Rfam database (http://www.sanger.ac.uk/software/Rfam) and the GenBank noncoding RNA database (http://blast.ncbi.nlm.nih.gov/) to annotate rRNA, RNA, snRNA and snoRNA. The other small RNAs were mapped to the zebrafish genome to perform distribution analysis and miRNA prediction using SOAP (http://soap.genomics.org.cn).

miRNA identification was performed by comparing the sequenced small RNAs with the known miRNAs of zebrafish or other animal species in miRBase v17.0. The miRNA precursors with a stem-loop structure were predicted by homologous comparison miRNA sequences to ESTs and to the genome sequence of the common carp with mfold software (http://mfold.rna.albany.edu/?q=mfold/RNA-Folding-Form).

RT-PCR and qRT-PCR

Total RNAs were isolated from embryonic and adult samples using TRIzol reagent (Invitrogen) and then treated with DNase I (RNase-free) (Takara). To detect miRNA expression, total RNAs were polyadenylated with ATP by E. coli poly (A) polymerase (Biolabs, New England). The polyadenylated total RNAs were reverse transcribed with M-MLV Reverse Transcription Reagents (Invitrogen) and a poly (T) primer ligated with a RACE adapter for miRNA quantitative assays. The RNAs without polyadenylation were reverse transcribed with oligo d (T) for the mRNA assay.

qRT-PCR was performed using the ABI 7500 sequence detection system (Applied Biosystems). U6 small nuclear RNA was used as an endogenous control for miRNAs. All reactions were run in triplicate and included no template controls for each gene. The relative amount of miRNA to U6 RNA was calculated using the 2−ΔΔCT method, and the level of significance was determined by one-way analysis of variance (ANOVA) with SPSS Statistics 17.0. All quantitative data presented were the mean ± SEM. All primers for RT-PCR and qPCR were shown in supplementary materials (Table S4).

Supporting Information

Figure S1.

Prediction of the fold-back structure of 10 conserved carp miRNA precursors. The precursor sequences were obtained by sequence alignment with the sequences of the common carp genome and ESTs. The mature miRNA sequences in the precursors are indicated in red.

https://doi.org/10.1371/journal.pone.0030925.s001

(TIF)

Figure S2.

Prediction of the precursor structure of 7 novel carp miRNAs. The precursor sequences were obtained by sequence alignment with the sequences of the common carp genome and ESTs. The mature miRNA sequences in the precursors are indicated in red.

https://doi.org/10.1371/journal.pone.0030925.s002

(TIF)

Figure S3.

The expression of muscle-related miRNAs during development of the common carp was detected by qRT-PCR. 0 h, fertilized oocytes; 24 h, 24 hpf embryos; 48 h, 48 hpf embryos; 1 d, 1 dph larva; 7 d, 7 dph larva; 15 d, 15 dph larva; 30 d, 30 dph larva. At least 5 animals were used for each timepoint. The level of miRNA expression in the 30 dph larva was defined as 1.

https://doi.org/10.1371/journal.pone.0030925.s003

(TIF)

Table S1.

Highly-conserved miRNAs identified in the common carp.

https://doi.org/10.1371/journal.pone.0030925.s004

(DOC)

Table S2.

miRNA families identified in the common carp.

https://doi.org/10.1371/journal.pone.0030925.s005

(DOC)

Table S3.

Comparison between the discovered miRNAs in the common carp and those found in the bighead carp and silver carp.

https://doi.org/10.1371/journal.pone.0030925.s006

(DOC)

Table S4.

Primer sequences for RT-PCR and qPCR assays.

https://doi.org/10.1371/journal.pone.0030925.s007

(DOC)

Acknowledgments

We thank American Journal Experts for editing the manuscript and Beijing Genomics Institute (Shenzhen, China) for high throughput sequencing of the small RNA library and helping partial bioinformatic analysis.

Author Contributions

Conceived and designed the experiments: CBT XS. Performed the experiments: XY LD Y. Li. Analyzed the data: CBT LD XZ. Contributed reagents/materials/analysis tools: XY Y. Lang. Wrote the paper: LD CBT XS.

References

  1. 1. Food and Agriculture Organization (FAO) (2002) Yearbook of fisheries statistics summary tables. Available: http://www.fao.org/fi/statist/summtab/default.asp.
  2. 2. Vandeputte M (2003) Selective breeding of quantitative traits in the common carp (Cyprinus carpio): a review. Aquat Living Resour 16: 399–407.
  3. 3. Sun X, Liang L (2004) A genetic linkage map of common carp (Cyprinus carpio L.) and mapping of a locus associated with cold tolerance. Aquaculture 38: 165–172.
  4. 4. Tanck MW, Palstra AP, van der Weerd M, Leffering CP, van der Poel JJ, et al. (2001) Segregation of microsatellite alleles and residual heterozygosity at single loci in homozygous androgenetic common carp (Cyprinus carpio L.). Genome 44: 743–751.
  5. 5. Cheng L, Liu L, Yu X, Wang D, Tong J (2010) A linkage map of common carp (Cyprinus carpio) based on AFLP and microsatellite markers. Anim Genet 41: 191–198.
  6. 6. Braun T, Gautel M (2011) Transcriptional mechanisms regulating skeletal muscle differentiation, growth and homeostasis. Nat Rev Mol Cell Biol 12: 349–361.
  7. 7. Black BL, Olson EN (1998) Transcriptional control of muscle development by myocyte enhancer factor-2 (MEF2) proteins. Annu Rev Cell Dev Biol 14: 167–196.
  8. 8. Ge Y, Chen J (2011) MicroRNAs in skeletal myogenesis. Cell Cycle 10: 441–448.
  9. 9. Güller I, Russell AP (2010) MicroRNAs in skeletal muscle: their role and regulation in development, disease and function. J Physiol 588: 4075–4087.
  10. 10. Ketting RF, Fischer SE, Bernstein E, Sijen T, Hannon GJ, et al. (2001) Dicer functions in RNA interference and in synthesis of small RNA involved in developmental timing in C. elegans. Genes Dev 15: 2654–2659.
  11. 11. Lee Y, Ahn C, Han J, Choi H, Kim J, et al. (2003) The nuclear RNase III Drosha initiates microRNA processing. Nature 425: 415–419.
  12. 12. Okamura K, Ishizuka A, Siomi H, Siomi MC (2004) Distinct roles for Argonaute proteins in small RNA-directed RNA cleavage pathways. Genes Dev 18: 1655–1666.
  13. 13. Du T, Zamore PD (2005) microPrimer: the biogenesis and function of microRNA. Development 132: 4645–4652.
  14. 14. Bartel DP (2004) MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116: 281–297.
  15. 15. Urbich C, Kuehbacher A, Dimmeler S (2008) Role of microRNAs in vascular diseases, inflammation, and angiogenesis. Cardiovasc Res 79: 581–588.
  16. 16. Blenkiron C, Miska EA (2007) miRNAs in cancer: approaches, aetiology, diagnostics and therapy. Hum Mol Genet 16: 106–113.
  17. 17. Bentwich I, Avniel A, Karov Y, Aharonov R, Gilad S, et al. (2005) Identification of hundreds of conserved and nonconserved human microRNAs. Nature Genet 37: 766–770.
  18. 18. Berezikov E, Guryev V, van de BJ, Wienholds E, Plasterk RH, et al. (2005) Phylogenetic shadowing and computational identification of human microRNA genes. Cell 120: 21–24.
  19. 19. Ruby JG, Jan C, Player C, Axtell MJ, Lee W, et al. (2006) Large-scale sequencing reveals 21U-RNAs and additional microRNAs and endogenous siRNAs in C. elegans. Cell 127: 1193–1207.
  20. 20. Wei Y, Chen S, Yang P, Ma Z, Kang L (2009) Characterization and comparative profiling of the small RNA transcriptomes in two phases of locust. Genome Bio 10: R6.
  21. 21. Fahlgren N, Howell MD, Kasschau KD, Chapman EJ, Sullivan CM, et al. (2007) High-throughput sequencing of Arabidopsis microRNAs: evidence for frequent birth and death of MIRNA genes. PLoS ONE 2: e219.
  22. 22. Cho IS, Kim J, Seo HY, Lim do H, Hong JS, et al. (2010) Cloning and characterization of microRNAs from porcine skeletal muscle and adipose tissue. Mol Biol Rep 37: 3567–3574.
  23. 23. Chi W, Tong C, Gan X, He S (2011) Characterization and comparative profiling of MiRNA transcriptomes in bighead carp and silver carp. PLoS One 6: e23549.
  24. 24. Ambros V, Bartel B, Bartel DP, Burge CB, Carrington JC, et al. (2003) A uniform system for microRNA annotation. RNA 9: 277–279.
  25. 25. Haider KH, Idris NM, Kim HW, Ahmed RP, Shujia J, et al. (2010) MicroRNA-21 is a key determinant in IL-11/Stat3 anti-apoptotic signalling pathway in preconditioning of skeletal myoblasts. Cardiovasc Res 88: 168–178.
  26. 26. Seok HY, Tatsuguchi M, Callis TE, He A, Pu W, et al. (2011) miR-155 inhibits the expression of MEF2A to repress skeletal muscle differentiation. J Biol Chem 286: 35339–35346.
  27. 27. Gagan J, Dey BK, Layer R, Yan Z, Dutta A (2011) MicroRNA-378 targets the myogenic repressor MyoR during myoblast differentiation. J Biol Chem 286: 19431–19438.
  28. 28. Crippa S, Cassano M, Messina G, Galli D, Galvez BG, et al. (2011) miR669a and miR669q prevent skeletal muscle differentiation in postnatal cardiac progenitors. J Cell Biol 193: 1197–1212.
  29. 29. Chen JF, Mandel EM, Thomson JM, Wu Q, Callis TE, et al. (2006) The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation. Nat Genet 38: 228–233.
  30. 30. Chen JF, Tao Y, Li J, Deng Z, Yan Z, et al. (2010) microRNA-1 and microRNA-206 regulate skeletal muscle satellite cell proliferation and differentiation by repressing Pax7. J Cell Biol 190: 867–879.
  31. 31. Yan D, Dong Xda E, Chen X, Wang L, Lu C, et al. (2009) MicroRNA-1/206 targets c-Met and inhibits rhabdomyosarcoma development. J Biol Chem 284: 29596–29604.
  32. 32. Mishima Y, Abreu-Goodger C, Staton AA, Stahlhut C, Shou C, et al. (2009) Zebrafish miR-1 and miR-133 shape muscle gene expression and regulate sarcomeric actin organization. Genes Dev 23: 619–632.
  33. 33. Ge Y, Sun Y, Chen J (2011) IGF-II is regulated by microRNA-125b in skeletal myogenesis. J Cell Biol 192: 69–81.
  34. 34. Naguibneva I, Ameyar-Zazoua M, Polesskaya A, Ait-Si-Ali S, Groisman R, et al. (2006) The microRNA miR-181 targets the homeobox protein Hox-A11 during mammalian myoblast differentiation. Nat Cell Biol 8: 278–284.
  35. 35. Anderson C, Catoe H, Werner R (2006) MIR-206 regulates connexin43 expression during skeletal muscle development. Nucleic Acids Res 34: 5863–5871.
  36. 36. Kim HK, Lee YS, Sivaprasad U, Malhotra A, Dutta A (2006) muscle-specific microRNA miR-206 promotes muscle differentiation. J Cell Biol 174(5): 677–668.
  37. 37. Taulli R, Bersani F, Foglizzo V, Linari A, Vigna E, et al. (2009) The muscle-specific microRNA miR-206 blocks human rhabdomyosarcoma growth in xenotransplanted mice by promoting myogenic differentiation. J Clin Invest 119: 2366–2378.
  38. 38. Liu J, Luo XJ, Xiong AW, Zhang ZD, Yue S, et al. (2010) MicroRNA-214 promotes myogenic differentiation by facilitating exit from mitosis via downregulation of proto-oncogene N-ras. J Biol Chem 285: 26599–26607.
  39. 39. Juan AH, Kumar RM, Marx JG, Young RA, Sartorelli V (2009) Mir-214-dependent regulation of the polycomb protein Ezh2 in skeletal muscle and embryonic stem cells. Mol Cell 36: 61–74.
  40. 40. Flynt AS, Li N, Thatcher EJ, Solnica-Krezel L, Patton JG (2007) Zebrafish miR-214 modulates Hedgehog signaling to specify muscle cell fate. Nat Genet 39: 259–263.
  41. 41. Cardinali B, Castellani L, Fasanaro P, Basso A, Alema S, et al. (2009) Microrna-221 and microrna-222 modulate differentiation and maturation of skeletal muscle cells. PLoS One 4: 7607.
  42. 42. Drummond MJ, Glynn EL, Fry CS, Dhanani S, Volpi E, et al. (2009) Essential amino acids increase microRNA-499, -208b, and -23a and downregulate myostatin and myocyte enhancer factor 2C mRNA expression in human skeletal muscle. 139: 2279–2284.
  43. 43. Sun Q, Zhang Y, Yang G, Chen X, Zhang Y, et al. (2008) Transforming growth factor-beta-regulated miR-24 promotes skeletal muscle differentiation. Nucleic Acids Res 36: 2690–2699.
  44. 44. Wong CF, Tellam RL (2008) MicroRNA-26a targets the histone methyltransferase Enhancer of Zeste homolog 2 during myogenesis. J Biol Chem 283: 9836–9843.
  45. 45. Crist CG, Montarras D, Pallafacchina G, Rocancourt D, Cumano A, et al. (2009) Muscle stem cell behavior is modified by microRNA-27 regulation of Pax3 expression. Proc Natl Acad Sci USA 106: 13383–13387.
  46. 46. Allen DL, Loh AS (2011) Posttranscriptional mechanisms involving microRNA-27a and b contribute to fast-specific and glucocorticoid-mediated myostatin expression in skeletal muscle. Am J Physiol Cell Physiol 300: C124–137.
  47. 47. Wang H, Garzon R, Sun H, Ladner KJ, Singh R, et al. (2008) NFkappaB-YY1-miR-29 regulatory circuitry in skeletal myogenesis and rhabdomyosarcoma. Cancer Cell 14: 369–381.
  48. 48. Buckingham M (2001) Skeletal muscle formation in vertebrates. Curr Opin Genet Dev 11: 440–448.
  49. 49. Huang TH, Zhu MJ, Li XY, Zhao SH (2008) Discovery of porcine microRNAs and profiling from skeletal muscle tissues during development. PLoS One 3: e3225.
  50. 50. Rosenberg MI, Georges SA, Asawachaicharn A, Analau E, Tapscott SJ (2006) MyoD inhibits Fstl1 and Utrn expression by inducing transcription of miR-206. J Cell Biol 175: 77–85.