BackgroundMicroRNAs (miRNAs) are a family of ~22 nucleotide small RNA molecules that regulate gene expression by fully or partially binding to their complementary sequences. Recently, a large number of miRNAs and their expression patterns have been identified in various species. However, to date no miRNAs have been reported to modulate muscle development in beef cattle.ResultsTotal RNAs from the Chinese Qinchuan bovine longissimus thoracis at fetal and adult stages were used to construct small RNA libraries for Solexa SBS technology sequencing. A total of 15,454,182 clean reads were obtained from the fetal bovine library and 13,558,164 clean reads from the adult bovine library. In total, 521 miRNAs including 104 novel miRNA candidates were identified. Furthermore, the nucleotide bias, base edit and family of the known miRNAs were also analyzed. Based on stem-loop qPCR, 25 high-read miRNAs were detected, and the results showed that bta-miRNA-206, miRNA-1, miRNA-133, miRNAn12, and miRNAn17 were highly expressed in muscle-related tissue or organs, suggesting that these miRNAs may play a role in the development of bovine muscle tissues.ConclusionsThis study confirmed the authenticity of 417 known miRNAs, discovered 104 novel miRNAs in bos taurus, and identified five muscle-specific miRNAs. The identification of novel miRNAs significantly expanded the repertoire of bovine miRNAs and could contribute to further studies on the muscle development of cattle.

DNA methylation is a key epigenetic modification in mammals and plays important roles in muscle development. We sampled longissimus dorsi muscle (LDM) from a well-known elite native breed of Chinese Qinchuan cattle living within the same environment but displaying distinct skeletal muscle at the fetal and adult stages. We generated and provided a genome-wide landscape of DNA methylomes and their relationship with mRNA and miRNA for fetal and adult muscle studies. Integration analysis revealed a total of 77 and 1,054 negatively correlated genes with methylation in the promoter and gene body regions, respectively, in both the fetal and adult bovine libraries. Furthermore, we identified expression patterns of high-read genes that exhibit a negative correlation between methylation and expression from nine different tissues at multiple developmental stages of bovine muscle-related tissue or organs. In addition, we validated the MeDIP-Seq results by bisulfite sequencing PCR (BSP) in some of the differentially methylated promoters. Together, these results provide valuable data for future biomedical research and genomic and epigenomic studies of bovine skeletal muscle that may help uncover the molecular basis underlying economically valuable traits in cattle. This comprehensive map also provides a solid basis for exploring the epigenetic mechanisms of muscle growth and development.

Studies of the miRNA expression profiles associated with the postnatal late development of skeletal muscle and IMF deposition are lacking in chicken. Here, we evaluated the patterns of muscle fiber growth and IMF deposition in breast muscle in the Chinese domestic breed called Gushi chicken, where we constructed four small RNA libraries from breast muscle tissues at 6, 14, 22, and 30 weeks. A total of 388 known miRNAs and 31 novel miRNAs were identified based on four small RNA libraries. Comparative analysis identified 92 significant differentially expressed (SDE) miRNAs based on six combinations. KEGG pathway analysis for the SDE miRNAs showed that metabolic pathways such as glycolysis and biosynthesis of amino acids were significantly enriched before 22 weeks, and pathways such as biosynthesis of unsaturated fatty acids and fatty acid elongation were significantly enriched after 22 weeks. This trend was consistent with the patterns of breast muscle fiber growth and IMF deposition in Gushi chickens. We also constructed miRNA-mRNA interaction networks related to breast muscle development and IMF deposition. The results showed that miRNAs such as gga-miR-1a-3p, and gga-miR-133a-5p may play important roles in breast muscle development, and miRNAs such as gga-miR-103-3p, and gga-miR-138-2-3p may have key roles in IMF deposition. This study determined the dynamic miRNA transcriptome in breast muscle tissue for the first time in Gushi chickens. The results provide a valuable resource for investigating the post-transcriptional regulation mechanisms during postnatal late development of breast muscle and IMF deposition and for evaluating the muscular disease.

BackgroundMicroRNAs (miRNAs) are small noncoding RNA molecules that serve as important post-transcriptional gene expression regulators by targeting messenger RNAs for post-transcriptional endonucleolytic cleavage or translational inhibition. miRNAs play important roles in many biological processes. Extensive high-throughput sequencing studies of miRNAs have been performed in several animal models. However, little is known about the diversity of these regulatory RNAs in goat (Capra hircus), which is one of the most important agricultural animals and the oldest domesticated species raised worldwide. Goats have long been used for their milk, meat, hair (including cashmere), and skins throughout much of the world.ResultsIn this study, two small RNA libraries were constructed based on dry period and peak lactation dairy goat mammary gland tissues and sequenced using the Illumina-Solexa high-throughput sequencing technology. A total of 346 conserved and 95 novel miRNAs were identified in the dairy goat. miRNAs expression was confirmed by qRT-PCR in nine tissues and in the mammary gland during different stages of lactation. In addition, several candidate miRNAs that may be involved in mammary gland development and lactation were found by comparing the miRNA expression profiles in different tissues and developmental stages of the mammary gland.ConclusionsThis study reveals the first miRNAs profile related to the biology of the mammary gland in the dairy goat. The characterization of these miRNAs could contribute to a better understanding of the molecular mechanisms of lactation physiology and mammary gland development in the dairy goat.

Laying performance is an important economic trait in hens, and this physiological process is largely influenced by the liver function. The livers of hens at 20- and 30-week-old stages were investigated using the next generation sequencing to identify the differences of microRNA expression profiles. Compared with the 20-week-old hens, 67 down- and 13 up-regulated microRNAs were verified to be significant differentially expressed (false discovery rate, FDR ≤ 0.05) (SDE) in the 30-week-old. We also identified 13 down- and 6 up-regulated novel differentially expressed (DE) microRNAs. miR-22-3p and miR-146b-5p, which exhibit critical roles in mammalian lipid metabolism, showed the most abundant expression and the highest fold-change, respectively. A total of 648 potential target genes of the SDE microRNAs were identified through an integrated analysis of microRNAs and the DE genes obtained in previous RNA-sequencing, including FADS1, FADS2, ELOVL6 and ACSL5, which are critical lipid metabolism-related regulators. Bioinformatic analyses revealed that target genes were mainly enriched in lipid-related metabolism processes. This work provides the first study of the expression patterns of hepatic microRNAs between 20- and 30-week old hens. The findings may serve as a fundamental resource for understanding the detailed functions of microRNAs in the molecular regulatory systems of lipid metabolism.

Domestication and breeding have reshaped the genomic architecture of chicken, but the retention and loss of genomic elements during these evolutionary processes remain unclear. We present the first chicken pan-genome constructed using 664 individuals, which identified an additional ∼66.5 Mb sequences that are absent from the reference genome (GRCg6a). The constructed pan-genome encoded 20,491 predicated protein-coding genes, of which higher expression level are observed in conserved genes relative to dispensable genes. Presence/absence variation (PAV) analyses demonstrated that gene PAV in chicken was shaped by selection, genetic drift, and hybridization. PAV-based GWAS identified numerous candidate mutations related to growth, carcass composition, meat quality, or physiological traits. Among them, a deletion in the promoter region of IGF2BP1 affecting chicken body size is reported, which is supported by functional studies and extra samples. This is the first time to report the causal variant of chicken body size QTL located at chromosome 27 which was repeatedly reported. Therefore, the chicken pan-genome is a useful resource for biological discovery and breeding. It improves our understanding of chicken genome diversity and provides materials to unveil the evolution history of chicken domestication.

Background/Aims: Poultry meat quality is affected by many factors, among which intramuscular fat (IMF) is predominant. IMF content affects the tenderness, juiciness, and favor of chicken. An increasing number of studies are focusing on the functions of microRNAs (miRNAs) during the adipogenic process. However, little is known about miRNAs associated with poultry IMF deposition, especially intramuscular adipocyte differentiation. Methods: The IMF content of two physiological stages was measured, and miRNA-Seq and RNA-Seq data were integrated and analyzed. A chicken intramuscular adipocyte cell differentiation model was constructed. A luciferase reporter assay, miRNA overexpression, and Oil Red O staining were used to confirm the targets of gga-miR-140-5p. Results: Our results showed that late-laying-period hens, which had a higher IMF content, exhibited lower global expression levels of miRNAs than juvenile hens. A total of 104 differentially expressed (DE) miRNAs were identified between the two groups. Integrated analysis of differentially expressed genes and DE miRNAs identified a total of 378 miRNA-mRNA pairs. Functional enrichment analysis revealed that these intersecting genes are involved in ubiquitin-mediated proteolysis, the peroxisome proliferator-activated receptor signaling pathway, glycerophospholipid metabolism, and fatty acid elongation and degradation pathways. Furthermore, we demonstrated that gga-miR-140-5p promoted intramuscular adipocyte differentiation via targeting retinoid X receptor gamma. Conclusion: Our findings may contribute to a more thorough understanding of chicken IMF deposition and the improvement of poultry meat quality.

There is a lack of understanding surrounding the molecular mechanisms involved in the development of chicken skeletal muscle in the late postnatal stage, especially in the regulation of breast muscle development related genes, pathways, miRNAs and other factors. In this study, 12 cDNA libraries and 4 small RNA libraries were constructed from Gushi chicken breast muscle samples from 6, 14, 22, and 30 weeks. A total of 15,508 known transcripts, 25,718 novel transcripts, 388 known miRNAs and 31 novel miRNAs were identified by RNA-seq in breast muscle at the four developmental stages. Through correlation analysis of miRNA and mRNA expression profiles, it was found that 417, 370, 240, 1,418, 496, and 363 negatively correlated miRNA–mRNA pairs of W14 vs. W6, W22 vs. W6, W22 vs. W14, W30 vs. W6, W30 vs. W14, and W30 vs. W22 comparisons, respectively. Based on the annotation analysis of these miRNA–mRNA pairs, we constructed the miRNA–mRNA interaction network related to biological processes, such as muscle cell differentiation, striated muscle tissue development and skeletal muscle cell differentiation. The interaction networks for signaling pathways related to five KEGG pathways (the focal adhesion, ECM-receptor interaction, FoxO signaling, cell cycle, and p53 signaling pathways) and PPI networks were also constructed. We found that ANKRD1, EYA2, JSC, AGT, MYBPC3, MYH11, ACTC1, FHL2, RCAN1, FOS, EGR1, and FOXO3, PTEN, AKT1, GADD45, PLK1, CCNB2, CCNB3 and other genes were the key core nodes of these networks, most of which are targets of miRNAs. The FoxO signaling pathway was in the center of the five pathway-related networks. In the PPI network, there was a clear interaction among PLK1 and CDK1, CCNB2, CDK1, and GADD45B, and CDC45, ORC1 and MCM3 genes. These results increase the understanding for the molecular mechanisms of chicken breast muscle development, and also provide a basis for studying the interactions between genes and miRNAs, as well as the functions of the pathways involved in postnatal developmental regulation of chicken breast muscle.