- Open Access
MiRNA sequencing of Embryonic Myogenesis in Chengkou Mountain Chicken
BMC Genomics volume 23, Article number: 571 (2022)
Skeletal muscle tissue is among the largest organ systems in mammals, essential for survival and movement. Embryonic muscle development determines the quantity and quality of muscles after the birth of an individual. MicroRNAs (miRNAs) are a significant class of non-coding RNAs that bind to the 3’UTR region of mRNA to regulate gene function. Total RNA was extracted from the leg muscles of chicken embryos in different developmental stages of Chengkou Mountain Chicken and used to generate 171,407,341 clean small RNA reads. Target prediction, GO, and KEGG enrichment analyses determined the significantly enriched genes and pathways. Differential analysis determined the significantly different miRNAs between chicken embryo leg muscles at different developmental stages. Meanwhile, the weighted correlation network analysis (WGCNA) identified key modules in different developmental stages, and the hub miRNAs were screened following the KME value.
The clean reads contained 2047 miRNAs, including 721 existing miRNAs, 1059 known miRNAs, and 267 novel miRNAs. Many genes and pathways related to muscle development were identified, including ERBB4, MEF2C, FZD4, the Wnt, Notch, and MAPK signaling pathways. The WGCNA established the greenyellow module and gga-miR-130b-5p for E12, magenta module and gga-miR-1643-5p for E16, purple module and gga-miR-12218-5p for E19, cyan module and gga-miR-132b-5p for E21.
These results lay a foundation for further research on the molecular regulatory mechanism of embryonic muscle development in Chengkou mountain chicken and provide a reference for other poultry and livestock muscle development studies.
Skeletal muscles develop from the mesoderm during the embryonic development of vertebrates such as birds and mammals . Poultry has two kinds of fine muscle fiber; red and white. The leg muscle mainly consists of red muscle fiber. Poultry muscle is an important protein source for humans. The main edible parts are the pectoral and leg muscles, and all are skeletal muscles. Skeletal muscle is an important vertebrate tissue accounting for 40% of body weight  and mainly functions in movement, postural support, breathing, and thermogenesis .
Moreover, skeletal muscle myogenesis is a complex biological process affected by various regulators [4,5,6,7]. During embryogenesis, the skeletal muscle forms in the vertebrate limb from progenitor cells originating in the somites . The process takes four stages to develop into mature muscle fibers. In the first stage, the mesenchymal stem cells from the mesoderm undergo terminal differentiation to form mononuclear myoblasts. The second stage involves the fusion of the mononuclear myoblasts to form a fusiform multinucleated myotube. Then, the third stage involves further differentiation of muscle tubes into muscle fibers. The last stage involves the growth and eventual maturation of muscle fibers [9, 10]. Thus, microRNAs can be employed to study the four stages of skeletal muscle development.
MicroRNAs are a class of endogenous small non-coding RNAs, approximately 19–22 nucleotides long, playing various important regulatory roles in cells, such as regulating post-transcriptional gene expression in plants and animals. Approximately 70% of mammalian miRNAs are located in transcription units (TUs) , but most miRNAs are located in introns. Each miRNA can target multiple genes, and several miRNAs can regulate the same gene. MiRNA maturation involves several processes , including 1) Primary miRNA formation by RNA polymerase II; 2) pre-miRNA generation by the nuclear RNase III enzyme, Drosha and its cofactor DGCR8 cleavaged. Next, 3) the pre-miRNA form miR/miR* duplex (an siRNA-like duplex) via the RNase III enzyme, Dicer; and lastly, 4) the mature single-stranded miRNA from the duplex is incorpoeated into the RNA-induced silencing complex (RISC). MiRNA functions by binding to AGO protein to form RISC, which then binds to the 3 ‘UTR region of the functional gene .
This study used transcriptomics, the study of gene expression at the RNA level, to analyze the skeletal muscles of Chengkou Mountain Chicken. The transcriptomics technology, also known as RNA-Seq, is an important method for studying cell phenotype and function. The transcriptome is the sum of transcription products of all genes in a cell, including mRNA, rRNA, tRNA, and non-coding RNA at a specific state or physiological condition of an organism. Therefore, a major feature of the transcriptome is its spatiotemporal specificity. The recent development of the next-generation high-throughput sequencing technology has tremendously updated the transcriptome sequencing technology. New technologies such as single-cell transcriptome sequencing and spatial transcriptome sequencing are discovered, has extended transcriptome research extends to the translatome and structurome , greatly enriching scientific output.
China is a vast country with diverse terrain and species. Chongqing, located in southwest China is the hub city of southwest China and the economic center of the upper reaches of the Yangtze River. The Chongqing climate is within the northern subtropical mountain area, characterized by many mountains and rivers. The climate is mild, with abundant rainfall, sufficient sunshine, and four distinct seasons, making it very suitable for the growth and breeding of animals and plants. The area has abundant high-quality local livestock and poultry varieties. The Chengkou mountain chicken is a local poultry variety with excellent product features, including resistance to coarse feeding, delicious meat, strong resistance, and high nutritional value. However, it has similar challenges (slow growth rate and low meat yield) as other local breeds . Therefore, studying and clarifying the biological mechanism of muscle development is necessary for improving the production performance of Chengkou mountain chicken and retaining its advantages.
This study analyzed Chengkou mountain chicken to explore the superiority of local chicken species in heritage performance by investigating four-stage chicken leg muscles: 12-day (E12), 16-day (E16), 19-day (E19), and 21-day (E21) embryos. Small-RNA sequencing unraveled differentially expressed miRNAs involved in embryo development. The study further characterized the differentially expressed miRNA in muscle development and established the enrichment functions and structure of miRNAs.
Overview of small-RNA sequencing
We constructed 12 cDNA libraries (E12-1, E12-2, E12-3, E16-1, E16-2, E16-3, E19-1, E19-2, E19-3, E21-1, E21-2, and E21-3) from embryo leg muscle to obtain complete miRNA transcripts of the chicken embryo. A total of 171,407,341 clean reads were generated from 12 cDNA libraries after dropping low-quality reads. Thus, reads containing over one low-quality base or unknown nucleotides (N); without 3’ adaptors; containing 5’ adaptors; containing 3’ and 5’ adaptors but no small RNA fragment between them; containing polyA in small RNA fragments < 18 nt were removed. The remaining high-quality reads of each duplicate were approximately 99%, and the proportion of clean tags was > 94% (Table 1). Transcripts per million(TPM) showed that miRNAs had different expressions (Fig. 1A), and samples correlation heatmaps showed high reproducibility between all samples (Fig. 1B). The length distribution of small RNA sequences was approximately 22 bp, consistent with conventional animal samples (Fig. S1). Nearly 5% of the tags aligned to non-coding RNAs (including rRNA, scRNA, snRNA, snoRNA, and tRNA) based on the GenBank (Release 209.0, Table S1) and the Rfam (version 11.0, Table S2) databases. The additional 95% of the tags were used for follow-up analysis. Moreover, over 85% of the transcripts had a high genome match (Fig. 1C). The reference area statistics showed consistent proportions of sense and antisense tags in the exon and intron regions (Fig. S2). The repeat alignment results are shown in Table S3.
Identification of microRNA
The tag abundance identified as miRNA was > 70%, and the miRNAs identified in this study were divided into three categories. (1) Existing miRNA; the miRNA obtained by comparing existing miRNAs of Chicken in the miRBase database. (2) Known miRNA; the miRNA obtained by comparing miRNA of other species in the miRBase database. (3) Novel miRNA; the new miRNA obtained by hairpin structure prediction based on comparing small RNA with reference sequences. About 7% of the miRNA was base-edited for each sample (Table S5), and 1,059 known miRNA were identified (Table S6). The first nucleotide bias within the existing miRNA tag sequences was U (Fig. S3), and the first nucleotide bias with known miRNA sequences was A and U (Fig. S4). In summary, 2047 miRNAs were identified by classifying 266,267 tags, including 721 existing miRNAs (Table S4), 1059 known miRNAs, and 267 novel miRNAs (Table S7, Fig. 1D). Fig. S5 shows the tag annotations for different samples.
MicroRNA different expression analysis
The PCA analysis of all the miRNA with < 1 TPM revealed 12 samples divided into four groups by time point (Fig. 2A). Samples E19 and E21 were very close, probably because both are in late embryonic development, indicating the reliability of the sequence data. Meanwhile, the cluster analysis showed that most miRNAs are expressed in the early stage embryos, indicating the importance of miRNA in early embryonic development (Fig. 2B). The edgeR software identified 196 differentially expressed miRNAs, including 27 in E12_vs_E16, 151 in E12_vs_E19, 171 in E12_vs_E21, 13 in E16_vs_E19, 32 in E16_vs_E21, and 1 in E19_vs_E21 (Fig. 2C). We performed an Upset plot on miRNAs from different stages to identify key miRNAs in muscle development. The results showed that the larger the time span, the more differentially expressed genes (Fig. 2D).
Weighted correlation network analysis (WGCNA) of miRNAs
The “WGCNA” R package  identified the key module miRNAs associated with and their regulatory roles in the different stages of embryonic muscle development. Essentially, 12 soft thresholds were used to ensure that the module conforms to scale-free distribution (Fig. 3A). Therefore, 15 modules (excluding unclassified miRNAs) were identified in the different module colors (Fig. 3B). The brighter the color of the intersection between the row and column, the closer the gene connection between the corresponding row and column. The Pearson correlation was stronger (Fig. 3C). Association analysis revealed significant correlations between E12 and greenyellow (r = 0.92), E16 and magenta (r = 0.58), E19 and purple (r = 0.58), and E21 and cyan (r = 0.59) (Fig. 3D). High KME (eigengene connectivity) values indicated hub genes with the most connections. The top3 miRNAs with the highest KME in each module were chosen as hub miRNAs of the corresponding modules (Table 2).
Functional analysis of miRNAs and co-expressed genes
The main function of miRNA is to bind mRNA regulate the expression of target genes. Top1 hub miRNA in four modules with different muscle development time points were selected for subsequent analysis of enriched miRNAs. Subsequently, target genes were predicted (total predicted miRNA target genes were listed in Table S9). The GO enrichment analysis of the miRNAs that target genes showed several muscle development-related GO terms in E12, including the regulation of muscle tissue development, muscle organ development, and muscle tissue development (Fig. 4A). In E16, the regulation of muscle cell differentiation and the regulation of vascular associated smooth muscle cell migration were enriched (Fig. 4B). Nonetheless, the regulation of vascular smooth muscle cell differentiation and the vascular smooth muscle cell differentiation were enriched in E19 (Fig. 4C). The skeletal muscle satellite cell differentiation and skeletal muscle cell differentiation were enriched in E21 (Fig. 4D).
A Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/) revealed the miRNA enriched pathways. Consequently, the mucin-type O-glycan biosynthesis and Notch signaling pathways were the most significantly enriched in E12 (Fig. 5A). The phosphatidylinositol signaling system, toll-like, and ErbB signaling pathways were the most significantly enriched in E16 (Fig. 5B). In E19, the adipocytokine signaling pathway was the most significantly enriched (Fig. 5C), while the circadian rhythm-fly and the ErbB signaling pathways were the most significantly enriched in E21 (Fig. 5D). Meanwhile, some star signaling pathways were significantly enriched during different developmental periods, such as the Wnt (E12), GnRH (E12, E19, E21), MAPK (E16, E21), and PPAR (E21) signaling pathways.
Co-expression network establishment
Based on previous analyses, we focused on a few pathways related to embryo muscle development combined with a previous mRNA study . The key miRNA-mRNA-pathway regulatory networks for different embryo development stages were built via Cytoscape 3.9.1. Thus, the Wnt signaling pathway was involved in development at E12, and the key genes included DAAM1, WNT16, PPP3R1, PRICKLE1, FZD4, PSEN1, PLCB1, CAMK2G, and RSPO1 (Fig. 6A). MAPK signaling was the most significant pathway for gga-miR-1643-5p at E16, and the important genes included MAP3K13, SOS2, BRAF, PAK1, RAP1B, RAC3, TRAF2, MYD88, NTRK2, NFATC1, RPS6KA5, MAP3K8, and MEF2C (Fig. 6B). Similarly, the most important pathway for E19 and E21 were adipocytokine and MAPK signaling pathways, respectively. The key genes included SOCS3, ACSL3, IRS1, STAT3, NFKBIA, PPARGC1A, and PPM1A for the adipocytokine signaling pathway. In contrast, ERBB4, MAP3K13, MAPK3, RASGRF2, PAK2, MAP2K2, NF1, MAP2K5, KDR, MAPKAPK2, DUSP6, NTRK2, CACNA2D1, PDGFRB, CACNB2 were key in the MAPK signaling pathway (Fig. 6C-D).
Validation of candidate miRNAs and miRNA-mRNA relationship
The miRNAs with high KME values and high expression in the key modules corresponding to different developmental stages of Chengkou Mountain Chicken embryos were selected for RT-qPCR. They included gga-miR-130b-5p, gga-miR-363-5p, gga-miR-338-5p, gga-miR-499-5p, gga-miR-1729-5p, gga-miR-26a-5p, gga-miR-30e-3p, and gga-miR-10b-5p. The RT-qPCR results and small RNA-Seq results were highly correlated (Fig. 7), confirming the accuracy of the sequencing results. Four pairs of miRNA-mRNA were randomly selected for RT-qPCR based on candidate miRNAs and target gene prediction results, and found that although there was a significant negative correlation of miRNA-mRNA expression relationship (Fig. 8).
Muscle development, growth, and regeneration occur throughout the life cycle of vertebrates. Myogenesis occurs in four consecutive, time-distinct but overlapping stages in amniotes, including embryo, fetus, neonate, and adult . Fetal and neonatal myogenesis is key for muscle growth and maturation. Adult myogenesis is necessary for postpartum growth and repairing damaged muscles . Primary and secondary fibers are produced during poultry embryonic and fetal development; after that, the number of myofibers remains stable  except during damage repair.
The small RNA sequence detection range is 18-30nt endogenous RNA, including miRNA, siRNA, and piRNA. However, the main objective of this study was to detect miRNAs related to muscle development of the chicken embryo to promote the genetic improvement of Chengkou mountain chickens. Therefore, 2047 miRNAs were detected by small RNA sequencing of the embryonic leg muscle of Chengkou Mountain chickens at different developmental time points. The differential analysis identified 196 differentially expressed miRNAs, indicating the significance of the 196 miRNAs in the muscle development of the chicken embryo.
A WGCNA systematic biological method described the patterns of gene association among the different samples and identified highly covariant gene sets, candidate biomarker genes, or therapeutic targets for use in medical and biological fields. WGCNA revealed the TERThigh-specific miR-17–92 cluster can targets biological processes enriched in the TERTlow cancer in a pan-cancer analysis . Elsewhere, RNA sequencing and WGCNA of cord blood samples from fetal growth-restricted cases and controls combined with maternal peripheral blood quantification revealed miR-42-5p and miR-1306-3p as potential fetal growth restrictors . In calves, WGCNA identified bta-miR-145 and bta-miR-199a-3p as important hub miRNAs regulating rumen development, immune system, and protein digestion . This technique also established that DYNLL2 and its target miR-148-3p are important regulators of chicken myogenesis . In this study, the WGCNA technique revealed key modules and hub genes related to leg muscles of embryos at four different development stages. Thus, many interesting genes for muscle development-related biological process and signaling pathways were identified through GO and KEGG hub gene target analysis. The WGCNA analysis identified 537 miRNAs in the greenyellow module, highly correlated (r = 0.92) with E12, 78 miRNAs in the magenta module for E16, 78 miRNAs in the purple module for E19, and 49 miRNAs in the cyan module for E21. These results account for a quarter of the total muscle miRNA, indicating that the progression of muscle development in the embryo is extremely important during this period. The results are consistent with an earlier study of E12 to E21, the period of rapid maturity of muscle fiber during the embryonic development of Chengkou mountain chicken .
Previous studies have shown that miRNAs are involved in muscle development; miR-222a and miR-126-5p significantly reduced the CPEB3 and FGFR3 mRNA levels in chicken embryo fibroblasts , emphasizing the significance of miRNA-target interactions in embryonic muscle regulation. A luciferase reporter gene assay showed that miR-1 targets chicken ACVR2B UTR directly, but network analysis predicted that ACVR2B targets gga-miR-101, gga-miR-1a and gga-miR-499 . Although gga-miRNA-454-3p does not affect primary myoblast differentiation, it inhibits differentiation by targeting the myotube-associated protein SBF2 . A study focused on the chicken skeletal muscle indicated that miR-29b-1-5p inhibits the proliferation of chicken primary myoblasts and promotes the differentiation of myoblasts via an effective target gene, ANKRD9 . Meanwhile, miR-29b-1-5p and miR-133a-5p are sponges for circFGFR2 in skeletal muscle proliferation and differentiation . MiR-133 was earlier established as a miRNA specifically expressed in skeletal muscles . High-throughput sequencing revealed novel_miR_158, novel_miR_144, novel_miR_291, and miR-205a as crucial miRNAs for skeletal muscle development in Bian chicken, suggesting their vital function in chicken growth .
The proliferation and differentiation of skeletal muscle satellite cells are crucial in skeletal muscle development, especially during repair after muscle injury. MiR-21-5p targets KLF3 and regulates skeletal muscle satellite cell proliferation and differentiation . MicroR-27b-3p regulates the proliferation and differentiation of chicken primary myoblasts by targeting MSTN . Moreover, miRNA-214 regulates chicken myoblast proliferation and differentiation by targeting TRMT61A . The miRNA-gene pairs gga-miR-499-5p/SOX6 and gga-miR-196-5p/CALM1 might affect muscle fiber performance using a miRNA-mRNA integrated analysis . All these miRNAs, except novel_miRNA (naming rules were different), were detected in this sequencing data, and gga-miR-499-5p was the candidate hub miRNA within the purple module for E19, suggesting its role in chicken muscle development.
RT-qPCR results by randomly selecting miRNAs and their target genes showed that there was extremely significant correlation between some miRNAs and mRNA, however, there were some miRNAs and mRNA with no obvious correlation. It confirmed that the targeting relationship between miRNA and mRNA is not unique. Besides, gga-miR-130b-5p, gga-miR-1643-5p, gga-miR-12218-5p, and gga-miR-132b-5p were critical for muscle development at E12, E16, E19, and E21 embryo stages in Chengkou mountain chickens. Multiple genes related to muscle development were identified by predicting the targets of these miRNAs. MEF2C, a member of the myocyte enhancer factor 2 family of MADS (MCM1, agamous, deficiens, serum response factor), is an important regulator of cardiac myogenesis and right ventricular development. This gene is mainly expressed in cardiac precursor cells before linear cardiac tube formation in mice, and MEF2C mutation prevents right ventricle formation in mice .
Meanwhile, MEF2C probably synergizes with MyoD through amplification to establish skeletal muscle commitment during cardiac and skeletal myogenesis . Activation of satellite cells regulates the repair of injured human skeletal muscles. However, after knocking out the MEF2A, MEF2C, and MEF2D genes, satellite cells only proliferated but failed to differentiate, showing the dependence of skeletal muscle regeneration on MEF2 . ERBB4 (Erb-B2 receptor tyrosine kinase 4), a Tyr protein kinase family member and the epidermal growth factor receptor subfamily, regulates muscle differentiation . The Janus kinase/signal transducer and activator of transcription (JAK/STAT) signaling cascade has also been identified as a key factor in myogenesis. Nonetheless, the STAT3 isoform is critical for satellite cell migration and myogenic differentiation because it mediates the expression of muscle-specific myogenic factors . The SOCS (Suppressor of Cytokine Signaling) family of proteins down-regulates STAT activation . Other target genes such as BRAF , DAAM1 , FZD4 [44, 45], and NF1 [46, 47] affect muscle development in the same or different ways.
GO, and KEGG enrichment results showed that some muscle development entries were significantly enriched, including the Wnt, ErbB, MAPK, and Notch signaling pathways. As previously predicted, DAAM1, FZD4, and WNT16 constitute the Wnt signaling pathway that regulates the critical ability of muscles to break down and reorganize fibers during development. Wnt signaling is involved in muscle remodeling . Consequently, correct activation of the Wnt signaling pathway is essential during the various steps of muscle formation . Therefore, deficiency of Wnt signaling effectors during pregnancy leads to marked tissue damage and muscle dysplasia . ERBB4 represents the ErbB signaling pathway, whose inhibition leads to non-denervated skeletal muscle growth in mice, but activation causes an opposite outcome .
In this study, the predicted members of the MAPK signaling pathway were PPM1A, NF1, MEF2C, and TRAF2. Ras-MAPK signaling promotes neuroactivity-dependent differentiation of slow muscle fibers in vivo . Besides, the p38 MAPK is activated during myoblast differentiation, and it also affects the activity of the MEF2 family of transcription factors, suggesting that p38 may be involved in the myogenic program . Early morphogenesis of skeletal muscles during chicken embryo development requires transient activation of the Notch signaling pathway to drive terminal differentiation of muscle progenitors . Notch and NRG signaling antagonistically regulate the synthesis and degradation of the cardiomyocyte extracellular matrix in a mouse trabecular model, which is critical for the individualization and rearrangement growth of trabecular units . Generally, the growth and development of organisms is a complex process, often regulated by several signaling pathways. Studies have shown that myogenic progenitor cell differentiation transitions from Notch to Wnt signaling. The temporal balance between Notch and Wnt signaling coordinates the precise progression of muscle precursor cells along the myogenic lineage pathway .
This work constructed the miRNA sequencing library of Chengkou mountain chicken, generating 2047 miRNAs and 196 differentially expressed miRNAs. Key modules, hub miRNAs, and targets corresponding to different chicken embryo developmental stages were identified through WGCNA and functional enrichment analysis. GO and KEGG enrichment analysis of target genes revealed several significantly enriched signaling pathways during embryonic muscle development, including the Wnt, ErbB, MAPK, and Notch signaling pathways. This report is highly consistent with previous mRNA sequencing results . Combining these reports can provide a more accurate molecular basis for exploring the embryonic muscle development of Chengkou Mountain Chicken and guide the genetic improvement of local breeds.
Chicken embryo incubation and tissue collection
This study used the Chengkou mountain chicken as the experimental animal. Chengkou mountain chicken breeding eggs were obtained from the Chongqing Xuanpeng Agricultural Development Co. Ltd Chongqing, China. The eggs were incubated at 37.8 °C and 55% humidity. Twelve chicken embryos were obtained from four time points (12, 16, 19, and 21 embryonic ages), with three replicates at each time point. The embryos were euthanized via cervical spine dislocations, and leg muscles were collected from the same sampling sites. The 12 samples were stored at -80 ℃ (wrapped in RNA protective solution (QIAGEN, Hilden, Germany)) for RNA extraction.
cDNA library construction and sequencing
The Trizol reagent (Invitrogen, USA) was used to extract total RNA from chicken embryo leg muscles during the four stages, following the manufacturer’s protocol. The RNA molecules within 18–30 nt were enriched by polyacrylamide gel electrophoresis (PAGE). Nucleic acid tests and gel electrophoresis assessed total RNA quality and purity. rRNA was removed from the total RNA using the Ribo-Zero rRNA removal kit (Epicentre, USA). The Illumina HiSeq™ 2500 (Illumina, CA, USA) was used for sequencing at the GENE DENOVO Biotechnology co. LTD (Guangzhou, China). The original data were filtered as follows to ensure quality. Reads containing: (1) > 1low-qualityy base (Q-value ≤ 20) or unknown nucleotides (N), (2) without 3’ adaptors, (3) containing 5’ adaptors, (4) containing 3’ and 5’ adaptors but no small RNA fragment in between, (5) containing polyA in small RNA fragments and < 18 nt were excluded. The clean reads were compared with the GenBank and Rfam species databases using the Blastall tool. Meanwhile, the chicken genome short reads were aligned using the tool Bowtie. Reads within the databases were divided and compared to avoid no mismatches. Then, the reserved unmapped reads were used for subsequent transcriptome analysis. The TPM determined the sample expression. Sample repeatability was tested via principal component analysis (PCA).
Identification of miRNAs
All clean tags were searched against the miRbase database (Release 22) to identify existing miRNAs and known miRNAs via alignment with other species. The novel miRNAs were identified according to their genome positions and hairpin structures as predicted by the Mireap_v0.2 software. The default parameters of the Mireap_v0.2 software were as follows: (1) 18nt minimal and (2) 26nt maximal miRNA sequence length. (3) Minimal, 20nt, and (4) maximal miRNA sequence length, 24nt. (5) Minimal depth of Drosha/Dicer cutting site, 3, (6) maximal copy number of miRNAs on reference, 20, and (7) maximal free energy allowed for a miRNA precursor, 18 kcal/mol. (8) Maximal space between miRNA and miRNA*, 35nt, (9) minimal space between miRNA and miRNA*, 14nt, and (10) maximal bulge between miRNA and miRNA*, 4nt. (11) Maximal asymmetry of miRNA/miRNA* duplex, 5nt, and (12) flank sequence length of miRNA precursor, 10nt. The tag annotation results were determined in this priority order: rRNA etc. > existing miRNA > existing miRNA edit > known miRNA > repeat > exon > novel miRNA > intron. The tags that were not annotated at any of the above molecules were recorded as unannotated.
MiRNA expression analysis
The total miRNA consisted of existing miRNA, known miRNA, and novel miRNAs, based on their expression in each sample. The miRNA expression was calculated and normalized to TPM. In addition, the expression of existing miRNA, known miRNA, and novel miRNA was also analyzed individually. The edgeR tool revealed the significantly different miRNAs based on the P value < 0.05 and |log2FC|> 1 threshold. RNAhybrid (Version 2.1.2) + svm_light (Version 6.01), Miranda (Version 3.3a) and TargetScan (Version 7.0) were used to predict targets. The intersection of the results was more credible and chosen as predicted miRNA target genes.
Function enrichment analysis
The miRNA-mRNA regulatory relationship was constructed to analyze the function of target genes and clarify the mechanism of miRNA involvement in chicken embryo muscle development. The Genes to GO term mapping database calculated the number of genes in each GO term and GO functional statistics . A hypergeometric test identified the significantly enriched GO entries compared with the entire genome background. A hypergeometric test used the KEGG databases to identify significantly enriched pathways against the entire genome background . The most important biochemical metabolic and signal transduction pathways were determined through enrichment analysis. The calculated p-values were subjected to FDR correction, and pathways with FDR ≤ 0.05 were considered statistically significant.
Verification and statistical analysis
Herein, eight miRNAs and four mRNAs were used to verify the sequencing results and the relationship between miRNAs and mRNAs expression levels via RT-qPCR. The primers were designed by Primer Premier (Table S8). RNA reverse transcription and real-time fluorescence quantitative PCR were performed as previously described , U6 and ACTB were used as housekeeping genes for miRNA and mRNA, respectively. The relative miRNA and mRNA expression were calculated via the 2−△△CT method , and data were expressed as mean ± standard deviation of the mean. Duncan’s Multiple Range Test was used for two-group comparisons in SPSS 23.0 (SPSS Inc., IL, USA). Graphics were plotted using GraphPad Prism 9 (GraphPad Software, CA, USA). P < 0.05 and P < 0.01 were considered statistically significant and extremely significant, respectively .
Availability of data and materials
The raw data has been submitted to the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA), and the accession number is SRP290982.
Christ B, Ordahl CP. Early stages of chick somite development. Anat Embryol (Berl). 1995;191:381–96.
Nie M, Deng Z-L, Liu J, Wang D-Z. Noncoding RNAs, Emerging Regulators of Skeletal Muscle Development and Diseases. Biomed Res Int. 2015;2015:676575.
Dumont NA, Bentzinger CF, Sincennes M-C, Rudnicki MA. Satellite Cells and Skeletal Muscle Regeneration. Compr Physiol. 2015;5:1027–59.
Buckingham M. Gene regulatory networks and cell lineages that underlie the formation of skeletal muscle. Proc Natl Acad Sci USA. 2017;114:5830–7.
Xu M, Chen X, Chen D, Yu B, Li M, He J, et al. Regulation of skeletal myogenesis by microRNAs. J Cell Physiol. 2020;235:87–104.
Li Z, Cai B, Abdalla BA, Zhu X, Zheng M, Han P, et al. LncIRS1 controls muscle atrophy via sponging miR-15 family to activate IGF1-PI3K/AKT pathway. J Cachexia Sarcopenia Muscle. 2019;10:391–410.
Luo W, Lin Z, Chen J, Chen G, Zhang S, Liu M, et al. TMEM182 interacts with integrin beta 1 and regulates myoblast differentiation and muscle regeneration. J Cachexia Sarcopenia Muscle. 2021;12:1704–23.
Buckingham M, Bajard L, Chang T, Daubas P, Hadchouel J, Meilhac S, et al. The formation of skeletal muscle: from somite to limb. J Anat. 2003;202:59–68.
Allen RE, Merkel RA, Young RB. Cellular aspects of muscle growth: myogenic cell proliferation. J Anim Sci. 1979;49:115–27.
Buckingham M, Vincent SD. Distinct and dynamic myogenic populations in the vertebrate embryo. Curr Opin Genet Dev. 2009;19:444–53.
Rodriguez A, Griffiths-Jones S, Ashurst JL, Bradley A. Identification of mammalian microRNA host genes and transcription units. Genome Res. 2004;14:1902–10.
Krol J, Loedige I, Filipowicz W. The widespread regulation of microRNA biogenesis, function and decay. Nat Rev Genet. 2010;11:597–610.
Huang Y, Shen XJ, Zou Q, Wang SP, Tang SM, Zhang GZ. Biological functions of microRNAs: a review. J Physiol Biochem. 2011;67:129–39.
Stark R, Grzelak M, Hadfield J. RNA sequencing: the teenage years. Nat Rev Genet. 2019;20:631–56.
Liu A, Ao X, Ma X, Wang W, Wang X, Lv X, et al. Effect of Different Rearing Modes on the Growth Performance and Meat Quality of Chengkou Mountain Chicken. J Southwest Univ (Natural Science Edition). 2018;40:1–7.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
Ren L, Liu A, Wang Q, Wang H, Dong D, Liu L. Transcriptome analysis of embryonic muscle development in Chengkou Mountain Chicken. BMC Genomics. 2021;22:431.
Biressi S, Molinaro M, Cossu G. Cellular heterogeneity during vertebrate skeletal muscle development. Dev Biol. 2007;308:281–93.
Murphy M, Kardon G. Origin of vertebrate limb muscle: the role of progenitor and myoblast populations. Curr Top Dev Biol. 2011;96:1–32.
Picard B, Lefaucheur L, Berri C, Duclos MJ. Muscle fibre ontogenesis in farm animal species. Reprod Nutr Dev. 2002;42:415–31.
Luo Z, Wang W, Li F, Songyang Z, Feng X, Xin C, et al. Pan-cancer analysis identifies telomerase-associated signatures and cancer subtypes. Mol Cancer. 2019;18:106.
Wang G, Yu J, Yang Y, Liu X, Zhao X, Guo X, et al. Whole-transcriptome sequencing uncovers core regulatory modules and gene signatures of human fetal growth restriction. Clin Transl Med. 2020;9:9.
Do DN, Dudemaine P-L, Fomenky BE, Ibeagha-Awemu EM. Integration of miRNA weighted gene co-expression network and miRNA-mRNA co-expression analyses reveals potential regulatory functions of miRNAs in calf rumen development. Genomics. 2019;111:849–59.
Li Y, Yuan P, Fan S, Zhai B, Jin W, Li D, et al. Weighted gene co-expression network indicates that the DYNLL2 is an important regulator of chicken breast muscle development and is regulated by miR-148a-3p. BMC Genomics. 2022;23:258.
Jebessa E, Ouyang H, Abdalla BA, Li Z, Abdullahi AY, Liu Q, et al. Characterization of miRNA and their target gene during chicken embryo skeletal muscle development. Oncotarget. 2018;9:17309–24.
Li T, Wu R, Zhang Y, Zhu D. A systematic analysis of the skeletal muscle miRNA transcriptome of chicken varieties with divergent skeletal muscle growth identifies novel miRNAs and differentially expressed miRNAs. BMC Genomics. 2011;12:186.
Chen M, Zhang S, Xu Z, Gao J, Mishra SK, Zhu Q, et al. MiRNA Profiling in Pectoral Muscle Throughout Pre- to Post-Natal Stages of Chicken Development. Front Genet. 2020;11:570.
Li Y, Zhai B, Yuan P, Fan S, Jin W, Li W, et al. MiR-29b-1-5p regulates the proliferation and differentiation of chicken primary myoblasts and analysis of its effective targets. Poult Sci. 2022;101:101557.
Chen X, Ouyang H, Wang Z, Chen B, Nie Q. A Novel Circular RNA Generated by FGFR2 Gene Promotes Myoblast Proliferation and Differentiation by Sponging miR-133a-5p and miR-29b-1-5p. Cells. 2018;7:E199.
Sweetman D, Goljanek K, Rathjen T, Oustanina S, Braun T, Dalmay T, et al. Specific requirements of MRFs for the expression of muscle specific microRNAs, miR-1, miR-206 and miR-133. Dev Biol. 2008;321:491–9.
Wu P, Zhou K, Zhang L, Li P, He M, Zhang X, et al. High-throughput sequencing reveals crucial miRNAs in skeletal muscle development of Bian chicken. Br Poult Sci. 2021;62:658–65.
Zhang D, Ran J, Li J, Yu C, Cui Z, Amevor FK, et al. miR-21-5p Regulates the Proliferation and Differentiation of Skeletal Muscle Satellite Cells by Targeting KLF3 in Chicken. Genes (Basel). 2021;12:814.
Zhang G, He M, Wu P, Zhang X, Zhou K, Li T, et al. MicroRNA-27b-3p Targets the Myostatin Gene to Regulate Myoblast Proliferation and Is Involved in Myoblast Differentiation. Cells. 2021;10:423.
Duan Y, Wu Y, Yin X, Li T, Chen F, Wu P, et al. MicroRNA-214 Inhibits Chicken Myoblasts Proliferation, Promotes Their Differentiation, and Targets the TRMT61A Gene. Genes (Basel). 2020;11:E1400.
Liu Y, Zhang M, Shan Y, Ji G, Ju X, Tu Y, et al. miRNA-mRNA network regulation in the skeletal muscle fiber phenotype of chickens revealed by integrated analysis of miRNAome and transcriptome. Sci Rep. 2020;10:10619.
Lin Q, Schwarz J, Bucana C, Olson EN. Control of mouse cardiac morphogenesis and myogenesis by transcription factor MEF2C. Science. 1997;276:1404–7.
Al Madhoun AS, Mehta V, Li G, Figeys D, Wiper-Bergeron N, Skerjanc IS. Skeletal myosin light chain kinase regulates skeletal myogenesis by phosphorylation of MEF2C. EMBO J. 2011;30:2477–89.
Liu N, Nelson BR, Bezprozvannaya S, Shelton JM, Richardson JA, Bassel-Duby R, et al. Requirement of MEF2A, C, and D for skeletal muscle regeneration. Proc Natl Acad Sci U S A. 2014;111:4109–14.
Paatero I, Veikkolainen V, Mäenpää M, Schmelzer E, Belting H-G, Pelliniemi LJ, et al. ErbB4 tyrosine kinase inhibition impairs neuromuscular development in zebrafish embryos. Mol Biol Cell. 2019;30:209–18.
Tierney MT, Aydogdu T, Sala D, Malecova B, Gatto S, Puri PL, et al. STAT3 signaling controls satellite cell expansion and skeletal muscle repair. Nat Med. 2014;20:1182–6.
Caldow MK, Steinberg GR, Cameron-Smith D. Impact of SOCS3 overexpression on human skeletal muscle development in vitro. Cytokine. 2011;55:104–9.
Inoue S-I, Takahara S, Yoshikawa T, Niihori T, Yanai K, Matsubara Y, et al. Activated Braf induces esophageal dilation and gastric epithelial hyperplasia in mice. Hum Mol Genet. 2017;26:4715–27.
Rosado M, Barber CF, Berciu C, Feldman S, Birren SJ, Nicastro D, et al. Critical roles for multiple formins during cardiac myofibril development and repair. Mol Biol Cell. 2014;25:811–27.
Rana K, Chiu MWS, Russell PK, Skinner JP, Lee NKL, Fam BC, et al. Muscle-specific androgen receptor deletion shows limited actions in myoblasts but not in myofibers in different muscles in vivo. J Mol Endocrinol. 2016;57:125–38.
Alexander MS, Kawahara G, Motohashi N, Casar JC, Eisenberg I, Myers JA, et al. MicroRNA-199a is induced in dystrophic muscle and affects WNT signaling, cell proliferation, and myogenic differentiation. Cell Death Differ. 2013;20:1194–208.
Wei X, Franke J, Ost M, Wardelmann K, Börno S, Timmermann B, et al. Cell autonomous requirement of neurofibromin (Nf1) for postnatal muscle hypertrophic growth and metabolic homeostasis. J Cachexia Sarcopenia Muscle. 2020;11:1758–78.
Sullivan K, El-Hoss J, Quinlan KGR, Deo N, Garton F, Seto JTC, et al. NF1 is a critical regulator of muscle development and metabolism. Hum Mol Genet. 2014;23:1250–9.
LeBoeuf B, Chen X, Garcia LR. WNT regulates programmed muscle remodeling through PLC-β and calcineurin in Caenorhabditis elegans males. Development. 2020;147:dev181305.
Girardi F, Le Grand F. Wnt Signaling in Skeletal Muscle Development and Regeneration. Prog Mol Biol Transl Sci. 2018;153:157–79.
Cisternas P, Henriquez JP, Brandan E, Inestrosa NC. Wnt signaling in skeletal muscle dynamics: myogenesis, neuromuscular synapse and fibrosis. Mol Neurobiol. 2014;49:574–89.
Ho BL, Goh Q, Nikolaou S, Hu L, Shay-Winkler K, Cornwall R. NRG/ErbB signaling regulates neonatal muscle growth but not neuromuscular contractures in neonatal brachial plexus injury. FEBS Lett. 2021;595:655–66.
Murgia M, Serrano AL, Calabria E, Pallafacchina G, Lomo T, Schiaffino S. Ras is involved in nerve-activity-dependent regulation of muscle genes. Nat Cell Biol. 2000;2:142–7.
Keren A, Tamir Y, Bengal E. The p38 MAPK signaling pathway: a major regulator of skeletal muscle development. Mol Cell Endocrinol. 2006;252:224–30.
Rios AC, Serralbo O, Salgado D, Marcelle C. Neural crest regulates myogenesis through the transient activation of NOTCH. Nature. 2011;473:532–5.
Del Monte-Nieto G, Ramialison M, Adam AAS, Wu B, Aharonov A, D’Uva G, et al. Control of cardiac jelly dynamics by NOTCH1 and NRG1 defines the building plan for trabeculation. Nature. 2018;557:439–45.
Brack AS, Conboy IM, Conboy MJ, Shen J, Rando TA. A temporal switch from notch to Wnt signaling in muscle stem cells is necessary for normal adult myogenesis. Cell Stem Cell. 2008;2:50–9.
Liu L, Yi J, Ray WK, Vu LT, Helm RF, Siegel PB, et al. Fasting differentially alters the hypothalamic proteome of chickens from lines with the propensity to be anorexic or obese. Nutr Diabetes. 2019;9:13.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36(Database issue):D480-484.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402–8.
Liu L, Xiao Q, Gilbert ER, Cui Z, Zhao X, Wang Y, et al. Whole-transcriptome analysis of atrophic ovaries in broody chickens reveals regulatory pathways associated with proliferation and apoptosis. Sci Rep. 2018;8:7231.
We are grateful to the editors and reviewers for their comments and suggestions. The authors would like to acknowledge Dr. Haoxiang Zhu for helping to operate experimental equipment in this research.
This work was financially supported by the National Natural Science Foundation of China (32002155), the Undergraduate Innovation and Entrepreneurship Training Program of Southwest University (202110635029), China Postdoctoral Science Foundation (2020M673157), Fundamental Research Funds for the Central Universities (SWU-KT22010), the Graduate Research&Innovation Project in Chongqing (CYS21126).
Ethics approval and consent to participate
The Animal Ethics Committee of Southwest University approved all the animal experiment procedures, which were performed humanly following the Management and Use of Laboratory Animals (2016) guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Alignment non-coding RNA in GenBank. Table S2. Alignment non-coding RNA in Rfam. Table S3. Alignment the repeat area. Table S4. Number of identified miRNA and tag abundance statistics. Table S5. Base-edited for each sample. Table S6. Statistics for identifying the known miRNA. Table S7. Statistics of novel miRNA. Table S8. The primers for the RT-qPCR amplification. Figure S1. Sample tags length distribution. (1) E12-1; (2) E12-2; (3) E12-3; (4) E16-1; (5) E16-2; (6); E16-3; (7) E19-1; (8) E19-2; (9) E19-3; (10) E21-1; (11) E21-2; (12) E21-3. Figure S2. Comparison of reference area statistics. Sample reference ratios distribution (1) E12-1; (2) E12-2; (3) E12-3; (4) E16-1; (5) E16-2; (6); E16-3; (7) E19-1; (8) E19-2; (9) E19-3; (10) E21-1; (11) E21-2; (12) E21-3. Replace different components ratios with different colors. Figure S3. The first nucleotide bias distribution. Replace different nucleotides with different colors. Figure S4. The first nucleotide bias with known miRNAs. Replace different nucleotides with different colors. Figure S5. Tags annotation for all samples. (1) E12-1; (2) E12-2; (3) E12-3; (4) E16-1; (5) E16-2; (6); E16-3; (7) E19-1; (8) E19-2; (9) E19-3; (10) E21-1; (11) E21-2; (12) E21-3. Replace different components ratios with different colors.
The predicted miRNA target genes.
About this article
Cite this article
Shi, J., Li, W., Liu, A. et al. MiRNA sequencing of Embryonic Myogenesis in Chengkou Mountain Chicken. BMC Genomics 23, 571 (2022). https://doi.org/10.1186/s12864-022-08795-z
- Chengkou mountain chicken
- Embryonic Myogenesis