Identification and characterization of microRNAs expressed in chicken skeletal muscle

Integrating miRNA and full-length transcriptome profiling to elucidate the mechanism of muscle growth in Muscovy ducks reveals key roles for miR-301a-3p/ANKRD1

BACKGROUND

The popularity of Muscovy ducks is attributed not only to their conformation traits but also to their slightly higher content of breast and leg meat, as well as their stronger-tasting meat compared to that of typical domestic ducks. However, there is a lack of comprehensive systematic research on the development of breast muscle in Muscovy ducks. In addition, since the number of skeletal muscle myofibers is established during the embryonic period, this study conducted a full-length transcriptome sequencing and microRNA sequencing of the breast muscle. Muscovy ducks at four developmental stages, namely Embryonic Day 21 (E21), Embryonic Day 27 (E27), Hatching Day (D0), and Post-hatching Day 7 (D7), were used to isolate total RNA for analysis.

RESULTS

A total of 68,161 genes and 472 mature microRNAs were identified. In order to uncover deeper insights into the regulation of mRNA by miRNAs, we conducted an integration of the differentially expressed miRNAs (known as DEMs) with the differentially expressed genes (referred to as DEGs) across various developmental stages. This integration allowed us to make predictions regarding the interactions between miRNAs and mRNA. Through this analysis, we identified a total of 274 DEGs that may serve as potential targets for the 68 DEMs. In the predicted miRNA‒mRNA interaction networks, let-7b, miR-133a-3p, miR-301a-3p, and miR-338-3p were the hub miRNAs. In addition, multiple DEMs also showed predicted target relationships with the DEGs associated with skeletal system development. These identified DEGs and DEMs as well as their predicted interaction networks involved in the regulation of energy homeostasis and muscle development were most likely to play critical roles in facilitating the embryo-to-hatchling transition. A candidate miRNA, miR-301a-3p, exhibited increased expression during the differentiation of satellite cells and was downregulated in the breast muscle tissues of Muscovy ducks at E21 compared to E27. A dual-luciferase reporter assay suggested that the ANKRD1 gene, which encodes a transcription factor, is a direct target of miR-301a-3p.

CONCLUSIONS

miR-301a-3p suppressed the posttranscriptional activity of ANKRD1, which is an activator of satellite cell proliferation, as determined with gain- and loss-of-function experiments. miR-301a-3p functions as an inducer of myogenesis by targeting the ANKRD1 gene in Muscovy ducks. These results provide novel insights into the early developmental process of black Muscovy breast muscles and will improve understanding of the underlying molecular mechanisms.

Genetic background and production periods shape the microRNA profiles of the gut in laying hens

There is growing evidence of the importance of miRNAs for intestinal functional properties and nutritional uptake. Comparative miRNAs profiles of the jejunal mucosa were established against two genetic backgrounds (Lohmann Brown-Classic (LB) and Lohmann LSL-Classic (LSL), which are similar in egg production but differ in physiological traits including mineral utilization, along the production periods of laying hens. The target genes of miRNAs higher expressed in LB vs. LSL (miR-126-3p, miR-214, miR-24-3p, miR-726-5p, miR-29b-3p) were enriched for energy pathways at all stages. The target genes of the miRNAs higher in LSL (miR-1788-5p, miR-103-3p, miR-22-5p, miR-221-3p, miR-375) were more enriched for immune and the bone signalling pathways. The most prominent expression differences were between 16 and 24 weeks of age before and after onset of laying. Our results evidence central roles of intestinal miRNAs as regulators of gene expression, influencing intestinal homeostasis and adaptation to environment in different strains and production phases.

Mechanism and Functions of Identified miRNAs in Poultry Skeletal Muscle Development – A Review

Development of the skeletal muscle goes through several complex processes regulated by numerous genetic factors. Although much efforts have been made to understand the mechanisms involved in increased muscle yield, little work is done about the miRNAs and candidate genes that are involved in the skeletal muscle development in poultry. Comprehensive research of candidate genes and single nucleotide related to poultry muscle growth is yet to be experimentally unraveled. However, over a few periods, studies in miRNA have disclosed that they actively participate in muscle formation, differentiation, and determination in poultry. Specifically, miR-1, miR-133, and miR-206 influence tissue development, and they are highly expressed in the skeletal muscles. Candidate genes such as CEBPB, MUSTN1, MSTN, IGF1, FOXO3, mTOR, and NFKB1, have also been identified to express in the poultry skeletal muscles development. However, further researches, analysis, and comprehensive studies should be made on the various miRNAs and gene regulatory factors that influence the skeletal muscle development in poultry. The objective of this review is to summarize recent knowledge in miRNAs and their mode of action as well as transcription and candidate genes identified to regulate poultry skeletal muscle development.

Gga-let-7f-3p promotes apoptosis in selenium deficiency-induced skeletal muscle by targeting selenoprotein K

Selenoprotein K (SELENOK) is primarily observed in the endoplasmic reticulum, and serves to maintain the normal physiological functions of skeletal muscle. Skeletal muscle development and regeneration are associated with significant changes in the expression of specific microRNAs (miRNAs). Downregulated SELENOK expression is observed in chicken muscles deficient of Se. However, the mechanisms of miRNA regulation of SELENOK expression remain elusive. Here, deep sequencing was used to detect the miRNA profiles of muscle in Se deficient (-Se group) and normal (C group) chickens. A dual-luciferase reporter assay was adopted to verify the relationship between SELENOK and gga-let-7f-3p. In addition, gga-let-7f-3p was either overexpressed or knocked-down in chicken myoblasts. Furthermore, the cells were treated with N-acetyl-l-cysteine (NAC) or hydrogen peroxide (H2O2) in order to probe the factors involved in oxidative stress, endoplasmic reticulum stress (ERS) and apoptosis, respectively. Relative to the C group, there were 132 differentially expressed miRNAs (including 57 upregulated and 75 downregulated) in the muscles of the -Se group. The dual-luciferase reporter assay showed that SELENOK was a primary target of gga-let-7f-3p. It was also observed that the overexpression or knock-down of gga-let-7f-3p significantly influenced the SELENOK expression. Moreover, NAC blocked mimics of ga-let-7f-3p, thus inducing oxidative stress, ERS and apoptosis. Simultaneously, gga-let-7f-3p inhibitors blocked the stimulant effects caused by H2O2 in chicken myoblasts. Furthermore, Se deficiency downregulated the SELENOK protein expression and induced oxidative stress, ERS and apoptosis in chicken muscles. In conclusion, the gga-let-7f-3p-SELENOK pathway played a pivotal role in Se deficiency mediated muscle injuries through the induction of oxidative stress and ERS, ultimately promoting apoptosis.

Characterization of miRNA transcriptome profiles related to breast muscle development and intramuscular fat deposition in chickens

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.

Sequencing and characterization of lncRNAs in the breast muscle of Gushi and Arbor Acres chickens

Chicken muscle quality is one of the most important factors determining the economic value of poultry, and muscle development and growth are affected by genetics, environment, and nutrition. However, little is known about the molecular regulatory mechanisms of long non-coding RNAs (lncRNAs) in chicken skeletal muscle development. Our study aimed to better understand muscle development in chickens and thereby improve meat quality. In this study, Ribo-Zero RNA-Seq was used to investigate differences in the expression profiles of muscle development related genes and associated pathways between Gushi (GS) and Arbor Acres (AA) chickens. We identified two muscle tissue specific expression lncRNAs. In addition, the target genes of these lncRNAs were significantly enriched in certain biological processes and molecular functions, as demonstrated by Gene Ontology (GO) analysis, and these target genes participate in five signaling pathway, as revealed by an analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Taken together, these data suggest that different lncRNAs might be involved in regulating chicken muscle development and growth and provide new insight into the molecular mechanisms of lncRNAs.

Genome-wide DNA methylation profiles reveal novel candidate genes associated with meat quality at different age stages in hens

Poultry meat quality is associated with breed, age, tissue and other factors. Many previous studies have focused on distinct breeds; however, little is known regarding the epigenetic regulatory mechanisms in different age stages, such as DNA methylation. Here, we compared the global DNA methylation profiles between juvenile (20 weeks old) and later laying-period (55 weeks old) hens and identified candidate genes related to the development and meat quality of breast muscle using whole-genome bisulfite sequencing. The results showed that the later laying-period hens, which had a higher intramuscular fat (IMF) deposition capacity and water holding capacity (WHC) and less tenderness, exhibited higher global DNA methylation levels than the juvenile hens. A total of 2,714 differentially methylated regions were identified in the present study, which corresponded to 378 differentially methylated genes, mainly affecting muscle development, lipid metabolism, and the ageing process. Hypermethylation of the promoters of the genes ABCA1, COL6A1 and GSTT1L and the resulting transcriptional down-regulation in the later laying-period hens may be the reason for the significant difference in the meat quality between the juvenile and later laying-period hens. These findings contribute to a better understanding of epigenetic regulation in the skeletal muscle development and meat quality of chicken.

Characterization of miR-206 Promoter and Its Association with Birthweight in Chicken

miRNAs have been widely investigated in terms of cell proliferation and differentiation. However, little is known about their effects on bird growth. Here we characterized the promoter of miR-206 in chicken and found that the preferable promoter was located in 1200 bp upstream of pri-miR-206. In this region, many key transcription factors, including MyoD, c-Myb, CEBPα/β, AP-4, RAP1, Brn2, GATA-1/2/3, E47, Sn, upstream stimulatory factor (USF) and CdxA, were predicted to bind and interact with miR-206 promoter. Overexpression of MyoD sharply increased miR-206 expression in both fibroblast and myoblast cells, and also the regulation in the myoblast cells was much stronger, indicating that miR-206 was regulated by MyoD combined with other muscle specific transcriptional factors. Aiming to further investigate the relationship between miR-206 mutation and transcriptional expression, total of 23 SNPs were identified in the two distinct bird lines by sequencing. Interestingly, the motif bound by MyoD was individually destroyed by G-to-C mutation located at 419 bp upstream of miR-206 precursor. Co-transfecting MyoD and miR-206 promoter in DF-1 cells, the luciferase activity of promoter containing homozygous GG types was significantly higher than CC ones (p < 0.05). Thus, this mutation caused low expression of miR-206. Consistently, eight variants including G-419C mutation exhibited a great effect on birthweight through maker-trait association analysis in F2 population (p < 0.05). Additionally, the regulation of miR-206 on embryo muscle mass mainly by increasing MyoG and muscle creatine kinase (MCK) expression (p < 0.05) with little change in MyoD, TMEM8C and myosin heavy chain (MHC). In conclusion, our findings provide a novel mutation destroying the promoter activity of miR-206 in birds and shed new light to understand the regulation mechanism of miR-206 on the embryonic muscle growth.

Human miR-221/222 in Physiological and Atherosclerotic Vascular Remodeling

A cluster of miR-221/222 is a key player in vascular biology through exhibiting its effects on vascular smooth muscle cells (VSMCs) and endothelial cells (ECs). These miRNAs contribute to vascular remodeling, an adaptive process involving phenotypic and behavioral changes in vascular cells in response to vascular injury. In proliferative vascular diseases such as atherosclerosis, pathological vascular remodeling plays a prominent role. The miR-221/222 cluster controls development and differentiation of ECs but inhibits their proangiogenic activation, proliferation, and migration. miR-221/222 are primarily implicated in maintaining endothelial integrity and supporting quiescent EC phenotype. Vascular expression of miR-221/222 is upregulated in initial atherogenic stages causing inhibition of angiogenic recruitment of ECs and increasing endothelial dysfunction and EC apoptosis. In contrast, these miRNAs stimulate VSMCs and switching from the VSMC "contractile" phenotype to the "synthetic" phenotype associated with induction of proliferation and motility. In atherosclerotic vessels, miR-221/222 drive neointima formation. Both miRNAs contribute to atherogenic calcification of VSMCs. In advanced plaques, chronic inflammation downregulates miR-221/222 expression in ECs that in turn could activate intralesion neoangiogenesis. In addition, both miRNAs could contribute to cardiovascular pathology through their effects on fat and glucose metabolism in nonvascular tissues such as adipose tissue, liver, and skeletal muscles.