miR-460b-5p promotes proliferation and differentiation of chicken myoblasts and targets RBM19 gene.

The meat production of broilers is crucial to economic benefits of broiler industries, while the slaughter performance of broilers is directly determined by skeletal muscle development. Hence, the broiler breeding for growth traits shows a great importance. As a kind of small noncoding RNA, microRNA (miRNA) can regulate the expression of multiple genes and perform a wide range of regulation in organisms. Currently, more and more studies have confirmed that miRNAs are closely associated with skeletal muscle development of chickens. Based on our previous miR-seq analysis (accession number: PRJNA668199), miR-460b-5p was screened as one of the key miRNAs probably involved in the growth regulation of chickens. However, the regulatory effect of miR-460b-5p on the development of chicken skeletal muscles is still unclear. Therefore, miR-460b-5p was further used for functional validation at the cellular level in this study. The expression pattern of miR-460b-5p was investigated in proliferation and differentiation stages of chicken primary myoblasts. It was showed that the expression level of miR-460b-5p gradually decreased from the proliferation stage (GM 50%) to the lowest at 24 h of differentiation. As differentiation proceeded, miR-460b-5p expression increased significantly, reaching the highest and stabilizing at 72 h and 96 h of differentiation. Through mRNA quantitative analysis of proliferation marker genes, CCK-8 and Edu assays, miR-460b-5p was found to significantly facilitate the transition of myoblasts from G1 to S phase and promote chicken myoblast proliferation. mRNA and protein quantitative analysis of differentiation marker genes, as well as the indirect immunofluorescence results of myotubes, revealed that miR-460b-5p significantly stimulated myotube development and promote chicken myoblast differentiation. In addition, the target relationship was validated for miR-460b-5p according to the dual-luciferase reporter assay and mRNA quantitative analysis, which indicates that miR-460b-5p was able to regulate RBM19 expression by specifically binding to the 3' UTR of RBM19. In summary, miR-460b-5p has positive regulatory effects on the proliferation and differentiation of chicken myoblasts, and RBM19 is a target gene of miR-460b-5p.

Exploration of Potential Target Genes of miR-24-3p in Chicken Myoblasts by Transcriptome Sequencing Analysis.

Broiler skeletal muscle growth is significantly influenced by miRNAs. Our earlier research demonstrated that miR-24-3p significantly suppressed the proliferation of chicken myoblasts while promoting their differentiation. The purpose of this study is to investigate miR-24-3p potential target genes in chickens. We collected myoblasts of Jinghai yellow chicken and transfected four samples with mimics of miR-24-3p and another four samples with mimic NC (negative control) for RNA-seq. We obtained 54.34 Gb of raw data in total and 50.79 Gb of clean data remained after filtering. Moreover, 11,635 genes were found to be co-expressed in these two groups. The mimic vs. NC comparison group contained 189 DEGs in total, 119 of which were significantly up-regulated and 70 of which were significantly down-regulated. Important biological process (BP) terminology such as nuclear chromosomal segregation, reproduction, and nuclear division were discovered by GO enrichment analysis for DEGs in the mimic vs. NC comparison group. KEGG pathway analysis showed that focal adhesion, cytokine-cytokine receptor interaction, the TGF-ß signaling pathway, and the MAPK signaling pathway were enriched in the top 20. Variation site analysis illustrated the SNP (single nucleotide polymorphisms) and INDEL (insertion-deletion) in the tested samples. By comparing the target genes predicted by miRDB (MicroRNA target prediction database) and TargetScan with the 189 DEGs found by the transcriptome sequencing, we discovered two up-regulated DEGs (NEURL1 and IQSEC3) and two down-regulated DEGs (REEP1 and ST6GAL1). Finally, we carried out qPCR experiments on eight DEGs and discovered that the qPCR results matched the sequencing outcomes. These findings will aid in identifying potential miR-24-3p target genes in chicken skeletal muscle and offer some new directions for upcoming research on broiler breeding.

Genome-Wide Analysis of the KLF Gene Family in Chicken: Characterization and Expression Profile.

The kruppel-like factor (KLF) gene family is a group of transcription factors containing highly conserved zinc-finger motifs, which play a crucial role in cell proliferation and differentiation. Chicken has been widely used as a model animal for analyzing gene function, however, little is known about the function of the KLF gene family in chickens. In this study, we performed genome-wide studies of chicken KLF genes and analyzed their biological and expression characteristics. We identified 13 KLF genes from chickens. Our phylogenetic, motif, and conserved domain analyses indicate that the KLF gene family has remained conserved through evolution. Synteny analysis showed the collinear relationship among KLFs, which indicated that they had related biomolecular functions. Interaction network analysis revealed that KLFs worked with 20 genes in biological processes. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that KLF2 was involved in Apelin and Forkhead Box O (FOXO) signaling pathways. Moreover, qPCR showed that 13 KLF genes were expressed in the nine selected tissues and displayed various gene expression patterns in chickens. RNA-seq showed that KLF3 and KLF10 genes were differentially expressed in the normal and high-fat diet fed groups, and KLF4, KLF5, KLF6, KLF7, KLF9, KLF12, and KLF13 genes were differentially expressed between undifferentiated and differentiated chicken preadipocytes. Besides, RNA-seq also showed that KLF genes displayed different expression patterns in muscle at 11 and 16 embryonic days old, and in 1-day-old chickens. These results indicated that the KLF genes were involved in the development of muscle and fat in chickens. Our findings provide some valuable reference points for the subsequent study of the function of KLF genes.

Transcriptome Sequencing Analysis of circRNA in Skeletal Muscle between Fast- and Slow-Growing Chickens at Embryonic Stages.

Skeletal muscle growth has always been the focus of the broiler industry, and circRNAs play a significant role in this process. We collected leg muscles of slow- and fast-growing Bian chicken embryos in the study at 14 (S14 and F14) and 20 (S20 and F20) days for RNA-seq. Finally, 123 and 121 differentially expressed circRNAs (DECs) were identified in S14 vs. F14 and S20 vs. F20, respectively. GO enrichment analysis for DECs obtained important biological process (BP) terms including nicotinate nucleotide biosynthetic process, nicotinate nucleotide salvage, and NAD salvage in S20 vs. F20 and protein mannosylation in S14 vs. F14. KEGG pathway analysis showed Wnt signaling pathway, Tight junction, Ubiquitin mediated proteolysis, and Notch signaling pathway were enriched in the top 20. Based on the GO and KEGG analysis results, we found some significant host genes and circRNAs such as NAPRT and novel_circ_0004547, DVL1 and novel_circ_0003578, JAK2 and novel_circ_0010289, DERA and novel_circ_0003082, etc. Further analysis found 19 co-differentially expressed circRNAs between the two comparison groups. We next constructed a circRNA-miRNA network for them, and some candidate circRNA-miRNA pairs related to skeletal muscle were obtained, such as novel_circ_0002153-miR-12219-5p, novel_circ_0003578-miR-3064-3p, and novel_circ_0010661-miR-12260-3p. These results would help to reveal the mechanism for circRNAs in skeletal muscle and also provide some guidance for the breeding of broilers.

Identification of miRNA-mRNA Networks Associated with Pigeon Skeletal Muscle Development and Growth.

The growth and development of skeletal muscle determine the productivity of pigeon meat production, and miRNA plays an important role in the growth and development of this type of muscle. However, there are few reports regarding miRNA regulating the growth and development of skeletal muscle in pigeons. To explore the function of miRNA in regulating the growth and development of pigeon skeletal muscle, we used RNA sequencing technology to study the transcriptome of pigeons at two embryonic stages (E8 and E13) and two growth stages (D1 and D10). A total of 32,527 mRNAs were identified in pigeon skeletal muscles, including 14,378 novel mRNAs and 18,149 known mRNAs. A total of 2362 miRNAs were identified, including 1758 known miRNAs and 624 novel miRNAs. In total, 839 differentially expressed miRNAs (DEmiRNAs) and 11,311 differentially expressed mRNAs (DEGs) were identified. STEM clustering analysis assigned DEmiRNAs to 20 profiles, of which 7 were significantly enriched (p-value < 0.05). These seven significantly enriched profiles can be classified into two categories. The first category represents DEmiRNAs continuously downregulated from the developmental stage to the growth stage of pigeon skeletal muscle, and the second category represents DEmiRNAs with low expression at the development and early growth stage, and significant upregulation at the high growth stage. We then constructed an miRNA−mRNA network based on target relationships between DEmiRNAs and DEGs belonging to the seven significantly enriched profiles. Based on the connectivity degree, 20 hub miRNAs responsible for pigeon skeletal muscle development and growth were identified, including cli-miR-20b-5p, miR-130-y, cli-miR-106-5p, cli-miR-181b-5p, miR-1-z, cli-miR-1a-3p, miR-23-y, cli-miR-30d-5p, miR-1-y, etc. The hub miRNAs involved in the miRNA−mRNA regulatory networks and their expression patterns during the development and growth of pigeon skeletal muscle were visualized. GO and KEGG enrichment analysis found potential biological processes and pathways related to muscle growth and development. Our findings expand the knowledge of miRNA expression in pigeons and provide a database for further investigation of the miRNA−mRNA regulatory mechanism underlying pigeon skeletal muscle development and growth.

Prediction of the Effect of Methylation in the Promoter Region of ZP2 Gene on Egg Production in Jinghai Yellow Chickens.

Egg production in chickens is a quantitative trait. The aim of this study was to investigate the effect of promoter methylation of the Zona pellucida 2 (ZP2) gene on egg production. Real-time fluorescence quantification showed that the expression of the ZP2 gene in the ovaries of 300-day-old Jinghai yellow chickens in the high-laying group was significantly higher than that in the low-laying group (p < 0.01). A series of deletion fragments of the ZP2 gene promoter in Jinghai yellow chickens had different promoter activities in DF-1 cells, and the core region of the ZP2 gene promoter was found to be between −1552 and −1348. Four CpG islands in the promoter region of the ZP2 gene were detected by software prediction. The overall degree of methylation of the ZP2-1 amplified fragment was negatively correlated with mRNA expression to some extent (R = −0.197); the overall degree of methylation of the ZP2-2 amplified fragment was also negatively correlated with mRNA expression to some extent (R = −0.264), in which the methylation of methylcytosine (mC)-9, mC-20, and mC-21 sites was significantly negatively correlated with mRNA expression (p < 0.05). In addition, the mC-20 and mC-21 sites are located on the Sp1 transcription factor binding site, and it is speculated that these two sites may be the main sites for regulating transcription. In summary, the methylation sites mC-20 and mC-21 of the ZP2 gene may inhibit the binding of Sp1 and DNA, affect the transcription of the ZP2 gene, and then affect the number of eggs produced by the Jinghai yellow chickens.

miRNA-10a-5p Targeting the BCL6 Gene Regulates Proliferation, Differentiation and Apoptosis of Chicken Myoblasts.

Proliferation, differentiation, and apoptosis are three essential stages in cell development, and miRNAs can achieve extensive regulation of cellular developmental processes by repressing the expression of target genes. According to our previous RNA-seq results, miRNA-10a-5p was differentially expressed at different periods in chicken myoblasts, revealing a possible association with muscle development. In this study, we concluded that miRNA-10a-5p inhibited chicken myoblasts' proliferation and differentiation and promoted chicken myoblasts' apoptosis by directly targeting BCL6, a critical transcription factor involved in muscle development and regeneration. Overexpression of BCL6 significantly facilitated myoblasts' proliferation and differentiation and suppressed myoblasts' apoptosis. On the contrary, knockdown of BCL6 significantly repressed myoblasts' proliferation and differentiation and induced myoblasts' apoptosis. The results above suggest that miRNA-10a-5p plays a potential role in skeletal muscle growth, development and autophagy by targeting the BCL6 gene. We first revealed the functions of miRNA-10a-5p and BCL6 in the proliferation, differentiation, and apoptosis of chicken myoblasts.

Transcriptome Integration Analysis at Different Embryonic Ages Reveals Key lncRNAs and mRNAs for Chicken Skeletal Muscle.

The growth and development of skeletal muscle at embryonic stages are vital and it directly affects the growth performance of chickens. Long non-coding RNA (lncRNA) plays an important role in this process. In the experiment, we collected the leg muscles of fast- and slow-growing Bian chickens both at 14- and 20-day embryo ages (14E and 20E) for RNA-seq. Finally, 292 and 347 differentially expressed (DE) lncRNAs were identified in F14vsF20 and S14vsS20, and 1,295 and 1,560 DE mRNAs were also screened, respectively. Then we constructed lncRNA-mRNA networks for the two groups, respectively, and found that 6 of the top 10 lncRNAs ranked with degree are same. GO analysis showed that 12 of the top 20 terms were same in the two comparison groups and most of them were related to energy metabolisms, such as cellular respiration and aerobic respiration. KEGG enrichment revealed that up to 16 pathways of the top 20 in F14vsF20 were same as that of S14vsS20 and most of them were related to growth, including citrate cycle (TCA cycle) and oxidative phosphorylation. Further analysis showed that there were 602 and 102 same DE mRNAs and DE lncRNAs between the two comparison groups. We then identified 442 lncRNA-mRNA pairs, including 201 mRNAs and 32 lncRNAs. Protein-Protein Interactions (PPI) network was predicted for the 201 mRNAs and three core networks were obtained using the plug-in MCODE of Cytoscape. Then the function of genes in the three core networks was further analyzed with ClueGo and they were mainly enriched in six groups of biological processes. On this basis, combined with KEGG pathways and lncRNA-mRNA networks, we identified several candidate lncRNAs and mRNAs. Among them, lncRNAs mainly include TCONS_00061389, TCONS_00025495, TCONS_00017622, TCONS_00216258 and TCONS_00084223, and mRNAs include PLK1, BUB1, TTK, NDUFS7 NDUFAB1, PDHA1, CDK1, SDHA, ACO2 and MDH1. The results would provide a foundation for further experiments on the role of lncRNAs in the regulation of muscle development. And it could also contribute to further clarify the regulatory mechanism of chicken skeletal muscle.

Identification of crucial circRNAs in skeletal muscle during chicken embryonic development.

BACKGROUND: Chicken provides humans with a large amount of animal protein every year, in which skeletal muscle plays a leading role. The embryonic skeletal muscle development determines the number of muscle fibers and will affect the muscle production of chickens. CircRNAs are involved in a variety of important biological processes, including muscle development. However, studies on circRNAs in the chicken embryo muscle development are still lacking. RESULTS: In the study, we collected chicken leg muscles at 14 and 20-day embryo ages both in the fast- and slow-growing groups for RNA-seq. We identified 245 and 440 differentially expressed (DE) circRNAs in the comparison group F14vsF20 and S14vsS20 respectively. GO enrichment analysis for the host genes of DE circRNAs showed that biological process (BP) terms in the top 20 related to growth in F14vsF20 were found such as positive regulation of transcription involved in G1/S phase of mitotic cell cycle, multicellular organismal macromolecule metabolic process, and multicellular organismal metabolic process. In group S14vsS20, we also found some BP terms associated with growth in the top 20 including actomyosin structure organization, actin cytoskeleton organization and myofibril assembly. A total of 7 significantly enriched pathways were obtained, containing Adherens junction and Tight junction. Further analysis of those pathways found three crucial host genes MYH9, YBX3, IGF1R in both fast- and slow-growing groups, three important host genes CTNNA3, AFDN and CREBBP only in the fast-growing group, and six host genes FGFR2, ACTN2, COL1A2, CDC42, DOCK1 and MYL3 only in the slow-growing group. In addition, circRNA-miRNA network also revealed some key regulation pairs such as novel_circ_0007646-miR-1625-5p, novel_circ_0007646-miR-1680-5p, novel_circ_0008913-miR-148b-5p, novel_circ_0008906-miR-148b-5p and novel_circ_0001640-miR-1759-3p. CONCLUSIONS: Comprehensive analysis of circRNAs and their targets would contribute to a better understanding of the molecular mechanisms in poultry skeletal muscle and it also plays an important guiding role in the next research.

Comparative Analysis of miRNA Expression Profiles in Skeletal Muscle of Bian Chickens at Different Embryonic Ages.

MicroRNAs (miRNAs) are widely involved in the growth and development of skeletal muscle through the negative regulation of target genes. In order to screen out the differentially expressed miRNAs (DEMs) associated with skeletal muscle development of Bian chickens at different embryonic ages, we used the leg muscles of fast-growing and slow-growing Bian chickens at the 14th and 20th embryonic ages (F14, F20, S14 and S20) for RNA-seq. A total of 836 known miRNAs were identified, and 121 novel miRNAs were predicted. In the F14 vs. F20 comparison group, 127 DEMs were screened, targeting a total of 2871 genes, with 61 miRNAs significantly upregulated and 66 miRNAs significantly downregulated. In the S14 vs. S20 comparison group, 131 DEMs were screened, targeting a total of 3236 genes, with 60 miRNAs significantly upregulated and 71 miRNAs significantly downregulated. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the predicted target genes were significantly enriched in 706 GO terms and 6 KEGG pathways in the F14 vs. F20 group and 677 GO terms and 5 KEGG pathways in the S14 vs. S20 group. According to the interaction network analysis, we screened five coexpressed DEMs (gga-miR-146a-3p, gga-miR-2954, gga-miR-34a-5p, gga-miR-1625-5p and gga-miR-18b-3p) with the highest connectivity degree with predicted target genes between the two comparison groups, and five hub genes (HSPA5, PKM2, Notch1, Notch2 and RBPJ) related to muscle development were obtained as well. Subsequently, we further identified nine DEMs (gga-let-7g-3p, gga-miR-490-3p, gga-miR-6660-3p, gga-miR-12223-5p, novel-miR-327, gga-miR-18a-5p, gga-miR-18b-5p, gga-miR-34a-5p and gga-miR-1677-3p) with a targeting relationship to the hub genes, suggesting that they may play important roles in the muscle development of Bian chickens. This study reveals the miRNA differences in skeletal muscle development between 14- and 20-day embryos of Bian chickens from fast- and slow-growing groups and provides a miRNA database for further studies on the molecular mechanisms of the skeletal muscle development in Bian chickens.