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

Genomic Insights into Molecular Regulation Mechanisms of Intramuscular Fat Deposition in Chicken

Intramuscular fat (IMF) plays an important role in the tenderness, water-holding capacity, and flavor of chicken meat, which directly affect meat quality. In recent years, regulatory mechanisms underlying IMF deposition and the development of effective molecular markers have been hot topics in poultry genetic breeding. Therefore, this review focuses on the current understanding of regulatory mechanisms underlying IMF deposition in chickens, which were identified by multiple genomic approaches, including genome-wide association studies, whole transcriptome sequencing, proteome sequencing, single-cell RNA sequencing (scRNA-seq), high-throughput chromosome conformation capture (HiC), DNA methylation sequencing, and m6A methylation sequencing. This review comprehensively and systematically describes genetic and epigenetic factors associated with IMF deposition, which provides a fundamental resource for biomarkers of IMF deposition and provides promising applications for genetic improvement of meat quality in chicken.

Revelation of genes associated with energy generating metabolic pathways in the fighter type Aseel chicken of India through skeletal muscle transcriptome sequencing

In this study, comparative analysis of skeletal muscle transcriptome was carried out for four biological replicates of Aseel, a fighter type breed and Punjab Brown, a meat type breed of India. The profusely expressed genes in both breeds were related to muscle contraction and motor activity. Differential expression analysis identified 961 up-regulated and 979 down-regulated genes in Aseel at a threshold of log2 fold change ≥ ±2.0 (padj<0.05). Significantly enriched KEGG pathways in Aseel included metabolic pathways and oxidative phosphorylation, with higher expression of genes associated with fatty acid beta-oxidation, formation of ATP by chemiosmotic coupling, response to oxidative stress, and muscle contraction. The highly connected hub genes identified through gene network analysis in the Aseel gamecocks were HNF4A, APOA2, APOB, APOC3, AMBP, and ACOT13, which are primarily associated with energy generating metabolic pathways. The up-regulated genes in Punjab Brown chicken were found to be related to muscle growth and differentiation. There was enrichment of pathways such as focal adhesion, insulin signaling pathway and ECM receptor interaction in these birds. The results presented in this study help to improve our understanding of the molecular mechanisms associated with fighting ability and muscle growth in Aseel and Punjab Brown chicken, respectively.

Weighted gene co-expression network analysis identified hub genes critical to fatty acid composition in Gushi chicken breast muscle

BACKGROUND

The composition and content of fatty acids in the breast muscle are important factors influencing meat quality. In this study, we investigated the fatty acid composition and content in the breast muscle of Gushi chickens at different developmental stages (14 weeks, 22 weeks, and 30 weeks). Additionally, we utilized transcriptomic data from the same tissue and employed WGCNA and module identification methods to identify key genes associated with the fatty acid composition in Gushi chicken breast muscle and elucidate their regulatory networks.

RESULTS

Among them, six modules (blue, brown, green, light yellow, purple, and red modules) showed significant correlations with fatty acid content and metabolic characteristics. Enrichment analysis revealed that these modules were involved in multiple signaling pathways related to fatty acid metabolism, including fatty acid metabolism, PPAR signaling pathway, and fatty acid biosynthesis. Through analysis of key genes, we identified 136 genes significantly associated with fatty acid phenotypic traits. Protein-protein interaction network analysis revealed that nine of these genes were closely related to fatty acid metabolism. Additionally, through correlation analysis of transcriptome data, we identified 51 key ceRNA regulatory networks, including six central genes, 7 miRNAs, and 28 lncRNAs.

CONCLUSION

This study successfully identified key genes closely associated with the fatty acid composition in Gushi chicken breast muscle, as well as their post-transcriptional regulatory networks. These findings provide new insights into the molecular regulatory mechanisms underlying the flavor characteristics of chicken meat and the composition of fatty acids in the breast muscle.

Genetic architecture and key regulatory genes of fatty acid composition in Gushi chicken breast muscle determined by GWAS and WGCNA

BACKGROUND

Fatty acids composition in poultry muscle is directly related to its tenderness, flavour, and juiciness, whereas its genetic mechanisms have not been elucidated. In this study, the genetic structure and key regulatory genes of the breast muscle fatty acid composition of local Chinese chicken, Gushi-Anka F2 resource population by integrating genome-wide association study (GWAS) and weighted gene co-expression network analysis (WGCNA) strategies. GWAS was performed based on 323,306 single nucleotide polymorphisms (SNPs) obtained by genotyping by sequencing (GBS) method and 721 chickens from the Gushi-Anka F2 resource population with highly variable fatty acid composition traits in the breast muscle. And then, according to the transcriptome data of the candidate genes that were obtained and phenotypic data of fatty acid composition traits in breast muscle of Gushi chickens at 14, 22, and 30 weeks of age, we conducted a WGCNA.

RESULTS

A total of 128 suggestive significantly associated SNPs for 11 fatty acid composition traits were identified and mapped on chromosomes (Chr) 2, 3, 4, 5, 13, 17, 21, and 27. Of these, the two most significant SNPs were Chr13:5,100,140 (P = 4.56423e-10) and Chr13:5,100,173 (P = 4.56423e-10), which explained 5.6% of the phenotypic variation in polyunsaturated fatty acids (PUFA). In addition, six fatty acid composition traits, including C20:1, C22:6, saturated fatty acid (SFA), unsaturated fatty acids (UFA), PUFA, and average chain length (ACL), were located in the same QTL intervals on Chr13. We obtained 505 genes by scanning the linkage disequilibrium (LD) regions of all significant SNPs and performed a WGCNA based on the transcriptome data of the above 505 genes. Combining two strategies, 9 hub genes (ENO1, ADH1, ASAH1, ADH1C, PIK3CD, WISP1, AKT1, PANK3, and C1QTNF2) were finally identified, which could be the potential candidate genes regulating fatty acid composition traits in chicken breast muscle.

CONCLUSION

The results of this study deepen our understanding of the genetic mechanisms underlying the regulation of fatty acid composition traits, which is helpful in the design of breeding strategies for the subsequent improvement of fatty acid composition in poultry muscle.

miRNA-seq analysis in skeletal muscle of chicken and function exploration of miR-24-3p

The regulation of skeletal muscle growth and development in chicken is complex. MicroRNAs (miRNAs) have been found to play an important role in the process, and more research is needed to further understand the regulatory mechanism of miRNAs. In this study, leg muscles of Jinghai yellow chickens at 300 d with low body weight (slow-growing group) and high body weight (fast-growing group) were collected for miRNA sequencing (miRNA-seq) and Bioinformatics analysis revealed 12 differentially expressed miRNAs (DEMs) between the two groups. We predicted 150 target genes for the DEMs, and GO and KEGG pathway analysis showed the target genes of miR-24-3p and novel_miR_133 were most enriched in the terms related to growth and development. Moreover, networks of DEMs and target genes showed that miR-24-3p and novel_miR_133 were the 2 core miRNAs. Hence, miR-24-3p was selected for further functional exploration in chicken primary myoblasts (CPMs) with molecular biology technologies including qPCR, cell counting kit-8 (CCK-8), 5-ethynyl-2'-deoxyuridine (EdU) and immunofluorescence. When proliferating CPMs were transfected with miR-24-3p mimic, the expression of cyclin dependent kinase inhibitor 1A (P21) was up-regulated and both CCK-8 and EdU assays showed that the proliferation of CPMs was inhibited. However, when the inhibitor was transfected into the proliferating CPMs, the opposite results were found. In differentiated CPMs, transfection with miR-24-3p mimic resulted in up regulation of MYOD, MYOG and MYHC after 48 h. Myotube areas also increased significantly compared to the mimic negative control (NC) group. When treated with inhibitor, differentiation CPMs produced the opposite effects. Overall, we revealed 2 miRNAs (novel_miR_133 and miR-24-3p) significantly related with growth and development and further proved that miR-24-3p could suppress the proliferation and promote differentiation of CPMs. The results would facilitate understanding the effects of miRNAs on the growth and development of chickens at the post-transcriptional level and could also have an important guiding role in yellow-feathered chicken breeding.

CircRNAs Related to Breast Muscle Development and Their Interaction Regulatory Network in Gushi Chicken

Circular RNAs (circRNAs) play a significant regulatory role during skeletal muscle development. To identify circRNAs during postnatal skeletal muscle development in chickens, we constructed 12 cDNA libraries from breast muscle tissues of Chinese Gushi chickens at 6, 14, 22, and 30 weeks and performed RNA sequencing. In total, 2112 circRNAs were identified, and among them 79.92% were derived from exons. CircRNAs are distributed on all chromosomes of chickens, especially chromosomes 1-9 and Z. Bioinformatics analysis showed that each circRNA had an average of 38 miRNA binding sites, 61.32% of which have internal ribosomal entry site (IRES) elements. Furthermore, in total 543 differentially expressed circRNAs (DE-circRNAs) were identified. Functional enrichment analysis revealed that DE-circRNAs source genes are engaged in biological processes and muscle development-related pathways; for example, cell differentiation, sarcomere, and myofibril formation, mTOR signaling pathway, and TGF-β signaling pathway, etc. We also established a competitive endogenous RNA (ceRNA) regulatory network associated with skeletal muscle development. The results in this report indicate that circRNAs can mediate the development of chicken skeletal muscle by means of a complex ceRNA network among circRNAs, miRNAs, genes, and pathways. The findings of this study might help increase the number of known circRNAs in skeletal muscle tissue and offer a worthwhile resource to further investigate the function of circRNAs in chicken skeletal muscle development.

Whole transcriptome RNA Sequencing Reveals the Global Molecular Responses and circRNA/lncRNA-miRNA-mRNA ceRNA Regulatory Network in Chicken Fat Deposition

Fat deposition is a vital factor affecting the economics of poultry production. Numerous studies on fat deposition have been done. However, the molecular regulatory mechanism is still unclear. In the present study, the whole-transcriptome RNA sequencing in abdominal fat, back skin, and liver both high- and low-abdominal fat groups was used to uncover the competitive endogenous RNA (ceRNA) regulation network related to chicken fat deposition. The results showed that differentially expressed (DE) genes in abdominal fat, back skin, liver were 1207(784 mRNAs, 330 lncRNAs, 41 circRNAs, 52 miRNAs), 860 (607 mRNAs, 166 lncRNAs, 26 circRNAs, 61 miRNAs), and 923 (501 mRNAs, 262 lncRNAs, 15 circRNAs, 145 miRNAs), respectively. The ceRNA regulatory network analysis indicated that the fatty acid metabolic process, monocarboxylic acid metabolic process, carboxylic acid metabolic process, glycerolipid metabolism, fatty acid metabolism, and peroxisome proliferator-activated receptor (PPAR) signaling pathway took part in chicken fat deposition. Meanwhile, we scan the important genes, FADS2, HSD17B12, ELOVL5, AKR1E2, DGKQ, GPAM, PLIN2, which were regulated by gga-miR-460b-5p, gga-miR-199-5p, gga-miR-7470-3p, gga-miR-6595-5p, gga-miR-101-2-5p. While these miRNAs were competitive combined by lncRNAs including MSTRG.18043, MSTRG.7738, MSTRG.21310, MSTRG.19577, and circRNAs including novel_circ_PTPN2, novel_circ_CTNNA1, novel_circ_PTPRD. This finding provides new insights into the regulatory mechanism of mRNA, miRNA, lncRNA, and circRNA in chicken fat deposition.

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

Background

The characteristics of muscle fibers determine the growth and meat quality of poultry. In this study, we performed a weighted gene co-expression network analysis (WGCNA) on the muscle fiber characteristics and transcriptome profile of the breast muscle tissue of Gushi chicken at 6, 14, 22, and 30 weeks.

Results

A total of 27 coexpressed biological functional modules were identified, of which the midnight blue module had the strongest correlation with muscle fiber and diameter. In addition, 7 hub genes were found from the midnight blue module, including LC8 dynein light chain 2 (DYNLL2). Combined with miRNA transcriptome data, miR-148a-3p was found to be a potential target miRNA of DYNLL2. Experiments on chicken primary myoblasts (CPMs) demonstrated that miR-148a-3p promotes the expression of myosin heavy chain (MYHC) protein by targeting DYNLL2, proving that it can promote differentiation of myoblasts.

Conclusions

This study proved that the hub gene DYNLL2 and its target miR-148-3p are important regulators in chicken myogenesis. These results provide novel insights for understanding the molecular regulation mechanisms related to the development of chicken breast muscle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08522-8.

Evaluation of Commercial Meat Products of Red Chicken Reared under LED Lights

Simple Summary

The aim of the study was to investigate the role of three different light-color temperatures of Light-Emitting Diodes (LEDs) (Neutral, Cool and Warm) on some attributes of breast meat. Various changes were observed in the physical and chemical characteristics of breast meat samples and the results obtained in fatty and volatile profiles indicate that the LEDs change the characteristics of meat.

Abstract

The objective of our study was to investigate the role of three different light-color temperatures of Light-Emitting Diodes (LEDs) [Neutral (K=3300−3700); Warm (K=3000−2500) and Cool (K=5500−6000)] on the qualitative attributes of breast meat obtained from male AZ Extra Heavy Red chickens. The comparison was made with meat deriving from chickens reared in the presence of classic neon lighting (Control). The meat was analyzed for the determination of both physical and chemical properties (cooking loss, moisture, total lipids and fatty acid composition). Furthermore, meat samples subjected to cooking were also analyzed for the identification of volatile compounds produced during the process; such evaluation was performed both immediately after cooking (T0) and after 7 days (T7) of cooked-meat storage at 4 °C. Cooking-loss values were higher for samples from chickens raised with Neutral LED (p < 0.05) compared to the other groups. For the fatty acid profiles of the meat, higher values were found for monounsaturated fatty acids (MUFAs) such as C18:1, C9 and C16:1 in Cool LED compared to the Control. Regarding the volatile profile of cooked meat, compounds belonging to the families of aldehydes, alcohols, ketones, and aromatic compounds were identified. Compounds belonging to the aldehyde family, such as hexanal, increased in Cool LED samples at T0 in comparison to the Control. On the other hand, the amounts of 1-Pentanol, 1-Octanol and 2-Octen-1-ol, which belong to the alcohol family, increased at T7 in Cool LED samples compared to the Warm LED. In conclusion, LED lighting showed to be effective in inducing significant variations on chicken breast meat ready to be introduced to the market, in particular regarding fatty acid profiles and the accumulation of volatile compounds. However more in-depth evaluation is needed for the identification of modifications regarding the sensorial sphere, which could have an impact on the consumer acceptability of the product.

Identification and Characterization of Circular RNAs in Association With the Deposition of Intramuscular Fat in Aohan Fine-Wool Sheep

Aohan fine-wool sheep (AFWS) is a high-quality fine-wool sheep breed that supplies wool and meat. Research is needed on the molecular mechanism behind intramuscular fat (IMF) deposition that greatly improves mutton quality. The widely expressed non-coding RNA is physiologically used in roles such as competitive endogenous RNA (ceRNA) that includes circular RNAs (circRNAs). Although circRNAs were studied in many fields, little research was devoted to IMF in sheep. We used the longissimus dorsi muscle of 2 and 12-month-old AWFS as research material to identify circRNAs related to IMF deposition in these sheep by RNA-seq screening for differentially expressed circRNAs in the two age groups. A total of 11,565 candidate circRNAs were identified, of which the 104 differentially expressed circRNAs in the two age groups were analyzed. Enrichment analysis was performed using Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes. The enriched pathways included lipid transport (GO:0006869), negative regulation of canonical Wnt signaling pathway (GO:0090090), fat digestion and absorption (ko04975), and sphingolipid metabolism (ko00600). The differentially expressed circRNAs included ciRNA455, circRNA9086, circRNA7445, circRNA4557, and others. The source genes involved in these pathways might regulate IMF deposition. We used the TargetScan and miRanda software for interaction analysis, and a network diagram of circRNA-miRNA interactions was created. CircRNA455-miR-127, circRNA455-miR-29a, circRNA455-miR-103, circRNA4557-mir149-5p, and circRNA2440-mir-23a might be involved in the IMF deposition process. The targeting relationship of circRNA4557-miR-149-5p was verified by a dual-luciferase reporter assay. The RT-qPCR results of seven randomly selected circRNAs were consistent with the sequencing results. This study provides additional information on circRNA regulation of IMF deposition in AFWS and is a useful resource for future research on this sheep breed.

Construction of miRNA-mRNA network in the differentiation of chicken preadipocytes

1. MicroRNAs (miRNAs) play key roles in regulating lipid metabolism, adipogenesis and fat deposition in chicken. To date, there are only a few miRNAs that had been confirmed to be involved in chicken adipogenesis. The detailed mechanisms by which miRNAs regulate chicken adipogenesis remain largely unknown.

2. To identify candidate miRNAs involved in chicken preadipocyte differentiation and explore potential mechanisms behind their functions, the following study analysed and identified miRNA and mRNA expression levels in undifferentiated and differentiated preadipocytes. Hub miRNA-mRNA interactions were identified, and the degree of connectivity of DE miRNAs in the network was established.

3. A total of 145 DE miRNAs and 660 DE mRNAs were identified between undifferentiated and differentiated preadipocytes. An miRNA-mRNA network was constructed, including 29 DE miRNAs and 155 DE mRNAs, forming 470 miRNA-mRNA interactions. Functional enrichment analysis showed that DE mRNAs in the network were significantly enriched in 712 biological processes and 13 KEGG pathways. Based on the connectivity degree, five DE miRNAs with higher degrees miR-195-x, gga-miR-200a-3p, gga-miR-135a-5p, novel-m0067-5p and novel-m0270-5p were identified as hub miRNAs. Fifty-eight DE mRNAs interacted with these five hub miRNAs and formed 70 miRNA-mRNA interactions.

4. This study constructed a miRNA-mRNA network associated with chicken preadipocyte differentiation and identified five hub miRNAs in the network. The findings identified the number of chicken adipogenic miRNAs and laid the foundation for elucidating the miRNA-mediated regulatory mechanism in chicken adipogenesis.

MiR-29b-1-5p regulates the proliferation and differentiation of chicken primary myoblasts and analysis of its effective targets

Several recent studies investigated the role of the miR-29 family in muscle development. However, only a few studies focused on chicken skeletal muscle. In the present study, cell cycle, 5-ethynyl-2'-deoxyuridine (EdU), cell counting kit-8 (CCK-8), and other assays indicated that miR-29b-1-5p can inhibit the proliferation of chicken primary myoblasts (CPMs); the western blot assay and immunofluorescence detection of MYHC demonstrated that miR-29b-1-5p can promote the differentiation of myoblasts. The functional enrichment analysis revealed that the target genes of miR-29b-1-5p may be involved in muscle tissue development, muscle organ development, and striated muscle tissue development, which are biological processes related to muscle development. The correlation analysis showed that these 6 genes, that is, ankyrin repeat domain 9 (ANKRD9), lactate dehydrogenase A (LDHA), transcription factor 12 (TCF12), FAT atypical cadherin 1 (FAT1), lin-9 homolog (LIN9), and integrin beta 3 binding protein (ITGB3BP), can be used as effective candidate target genes of miR-29b-1-5p. Moreover, miR-29b-1-5p inhibits the expression of ANKRD9 by directly binding the 3'UTR of ANKRD9. Overall, these data indicate that miR-29b-1-5p inhibits the proliferation of primary chicken myoblasts, stimulates their differentiation, and is involved in the process of muscle development and that its effective target gene is ANKRD9. This study identified the molecular mechanism of miR-29b-1-5p in chicken muscle development.

Integrated mRNA and miRNA transcriptome analysis provides novel insights into the molecular mechanisms underlying goose pituitary development during the embryo-to-hatchling transition

It is well established that the endocrine system plays a pivotal role in preparing the avian embryos for the abrupt switch from chorioallantoic to pulmonary respiration during the critical embryo-to-hatchling transition. However, as the master gland of the endocrine system, there has been little research focusing on the molecular mechanisms controlling the development and function of the pituitary gland during the peri-hatch period in birds. In the present study, we aimed to determine the genome-wide mRNA and miRNA transcriptome profiles of the pituitary during the embryo-to-hatchling transition period from embryonic day 22 (E22) to post-hatching day 6 (P6) in the goose (Anser cygnoides). Of note, expression of Anser_cygnoides_newGene_32456 and LOC106031011 were significantly different among these 4 stages (i.e., E22, E26, P2, and P6). Meanwhile, the neuroactive ligand-receptor interaction pathway was significantly enriched by the DEGs commonly identified among three pairwise comparisons. At the miRNA transcriptome level, there were not commonly identified DE miRNAs among these 4 stages, while the 418 of their predicted target genes were mutually shared. Both the target genes of DE miRNAs in each comparison and these 418 shared target genes were significantly enriched in the ECM-receptor interaction and focal adhesion pathways. In the predicted miRNA-mRNA interaction networks of these 2 pathways, novel_miRNA_467, novel_miRNA_154, and novel_miRNA_340 were the hub miRNAs. In addition, multiple DE miRNAs also showed predicted target relationships with the DEGs associated with extracellular matrix (ECM) components. Among them, expression of novel_miR_120, tgu-miR-92-3p, and novel_miR_398 was significantly negatively correlated with that of LAMC3 (laminin subunit gamma3), suggesting that these miRNAs may regulate pituitary tissue remodeling and functional changes through targeting LAMC3 during development. These identified DE mRNAs and miRNAs as well as their predicted interaction networks involved in regulation of tissue remodeling and cellular functions were most likely to play critical roles in facilitating the embryo-to-hatchling transition. These results provide novel insights into the early developmental process of avian pituitary gland and will help better understand the underlying molecular mechanisms.

Effect of the ACAA1 Gene on Preadipocyte Differentiation in Sheep

Acetyl-CoA acyltransferase 1 (ACAA1) functions as a key regulator of fatty acid β-oxidation in peroxisomes by catalyzing the cleavage of 3-ketoacyl-CoA to acetyl-CoA and acyl-CoA, which participate in the extension and degradation of fatty acids. Thus, ACAA1 is an important regulator of lipid metabolism and plays an essential role in fatty acid oxidation and lipid metabolism. Our previous study findings revealed that ACAA1 is closely associated with the peroxisome proliferator-activated receptor (PPAR) signaling and fatty acid metabolism pathways, which are involved in fat deposition in sheep, leading to our hypothesis that ACAA1 may be involved in fat deposition by regulating lipid metabolism. However, the associated molecular mechanism remains unclear. In the present study, to assess the potential function of ACAA1 in sheep preadipocyte differentiation, we knocked down and overexpressed ACAA1 in sheep preadipocytes and evaluated the pattern of ACAA1 gene expression during preadipocyte differentiation by qRT-PCR. ACAA1 was significantly expressed in the early stage of adipocyte differentiation, and then its expression decreased. ACAA1 deficiency increased lipid accumulation and the triglyceride content and promoted sheep preadipocyte differentiation, whereas ACAA1 overexpression inhibited adipogenesis and decreased lipid accumulation and the triglyceride content. Simultaneously, we demonstrated that ACAA1 deficiency upregulated the expressions of the adipogenic marker genes PPARγ and C/EBPα in sheep preadipocytes, but ACAA1 overexpression inhibited the expressions of these markers, indicating that ACAA1 affects lipid metabolism by regulating adipogenic marker genes. Our results may promote a better understanding of the regulation of adipogenesis by ACAA1.

Characteristics and expression profiles of circRNAs during abdominal adipose tissue development in Chinese Gushi chickens

Circular RNAs (circRNAs) play important roles in adipogenesis. However, studies on circRNA expression profiles associated with the development of abdominal adipose tissue are lacking in chickens. In this study, 12 cDNA libraries were constructed from the abdominal adipose tissue of Chinese domestic Gushi chickens at 6, 14, 22, and 30 weeks. A total of 1,766 circRNAs were identified by Illumina HiSeq 2500 sequencing. These circRNAs were primarily distributed on chr1 through chr10 and sex chromosomes, and 84.95% of the circRNAs were from gene exons. Bioinformatic analysis showed that each circRNA has 35 miRNA binding sites on average, and 62.71% have internal ribosome entry site (IRES) elements. Meanwhile, these circRNAs were primarily concentrated in TPM < 0.1 and TPM > 60, and their numbers accounted for 18.90% and 80.51%, respectively, exhibiting specific expression patterns in chicken abdominal adipose tissue. In addition, 275 differentially expressed (DE) circRNAs were identified by comparison analysis. Functional enrichment analysis showed that the parental genes of DE circRNAs were primarily involved in biological processes and pathways related to lipid metabolism, such as regulation of fat cell differentiation, fatty acid homeostasis, and triglyceride homeostasis, as well as fatty acid biosynthesis, fatty acid metabolism, and glycerolipid metabolism. Furthermore, ceRNA regulatory networks related to abdominal adipose development were constructed. The results of this study indicated that circRNAs can regulate lipid metabolism, adipocyte proliferation and differentiation, and cell junctions during abdominal adipose tissue development in chickens through complex ceRNA networks between circRNAs, miRNAs, genes, and pathways. The results of this study may help to expand the number of known circRNAs in abdominal adipose tissue and provide a valuable resource for further research on the function of circRNAs in chicken abdominal adipose tissue.

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.

Molecular Regulation of Lipogenesis, Adipogenesis and Fat Deposition in Chicken

In the poultry industry, excessive fat deposition is considered an undesirable factor, affecting feed efficiency, meat production cost, meat quality, and consumer’s health. Efforts to reduce fat deposition in economically important animals, such as chicken, can be made through different strategies; including genetic selection, feeding strategies, housing, and environmental strategies, as well as hormone supplementation. Recent investigations at the molecular level have revealed the significant role of the transcriptional and post-transcriptional regulatory networks and their interaction on modulating fat metabolism in chickens. At the transcriptional level, different transcription factors are known to regulate the expression of lipogenic and adipogenic genes through various signaling pathways, affecting chicken fat metabolism. Alternatively, at the post-transcriptional level, the regulatory mechanism of microRNAs (miRNAs) on lipid metabolism and deposition has added a promising dimension to understand the structural and functional regulatory mechanism of lipid metabolism in chicken. Therefore, this review focuses on the progress made in unraveling the molecular function of genes, transcription factors, and more notably significant miRNAs responsible for regulating adipogenesis, lipogenesis, and fat deposition in chicken. Moreover, a better understanding of the molecular regulation of lipid metabolism will give researchers novel insights to use functional molecular markers, such as miRNAs, for selection against excessive fat deposition to improve chicken production efficiency and meat quality.

LncRNAs and their regulatory networks in breast muscle tissue of Chinese Gushi chickens during late postnatal development

BACKGROUND

Chicken skeletal muscle is an important economic product. The late stages of chicken development constitute the main period that affects meat production. LncRNAs play important roles in controlling the epigenetic process of growth and development. However, studies on the role of lncRNAs in the late stages of chicken breast muscle development are still lacking. In this study, to investigate the expression characteristics of lncRNAs during chicken muscle development, 12 cDNA libraries were constructed from Gushi chicken breast muscle samples from 6-, 14-, 22-, and 30-week-old chickens.

RESULTS

A total of 1252 new lncRNAs and 1376 annotated lncRNAs were identified. Furthermore, 53, 61, 50, 153, 117, and 78 DE-lncRNAs were found in the W14 vs. W6, W22 vs. W14, W22 vs. W6, W30 vs. W6, W30 vs. W14, and W30 vs. W22 comparison groups, respectively. After GO enrichment analysis of the DE-lncRNAs, several muscle development-related GO terms were found in the W22 vs. W14 comparison group. Moreover, it was found that the MAPK signaling pathway was one of the most frequently enriched pathways in the different comparison groups. In addition, 12 common target DE-miRNAs of DE-lncRNAs were found in different comparison groups, some of which were muscle-specific miRNAs, such as gga-miR-206, gga-miR-1a-3p, and miR-133a-3p. Interestingly, the precursors of four newly identified miRNAs were found to be homologous to lncRNAs. Additionally, we found some ceRNA networks associated with muscle development-related GO terms. For example, the ceRNA networks contained the DYNLL2 gene with 12 lncRNAs that targeted 2 miRNAs. We also constructed PPI networks, such as IGF-I-EGF and FZD6-WNT11.

CONCLUSIONS

This study revealed, for the first time, the dynamic changes in lncRNA expression in Gushi chicken breast muscle at different periods and revealed that the MAPK signaling pathway plays a vital role in muscle development. Furthermore, MEF2C and its target lncRNA may be involved in muscle regulation through the MAPK signaling pathway. This research provided valuable resources for elucidating posttranscriptional regulatory mechanisms to promote the development of chicken breast muscles after hatching.

The Landscape of DNA Methylation Associated With the Transcriptomic Network of Intramuscular Adipocytes Generates Insight Into Intramuscular Fat Deposition in Chicken

Intramuscular fat (IMF), which regulated by genetics, nutrition and environment is an important factor that influencing meat quality. Up to now, the epigenetic regulation mechanism underlying poultry IMF deposition remains poorly understood. Here, we focused on the DNA methylation, which usually regulate genes in transcription level. To look into the essential role of DNA methylation on the IMF deposition, chicken intramuscular preadipocytes were isolated and cultured in vitro, and a model of intramuscular adipocyte differentiation was constructed. Combined the whole genome bisulfite sequencing (WGBS) and RNA-Seq technologies, we identified several methylated genes, which mainly affecting fatty acid metabolism and muscle development. Furthermore, we reported that DNA methylation regulate intramuscular adipogenesis by regulating the genes, such as collagen, type VI, alpha 1 (COL6A1) thus affecting IMF deposition. Overexpression of COL6A1 increases the lipid droplet and inhibits cell proliferation by regulating CHAD and CAMK2 in intramuscular adipocytes, while knockdown of COL6A1 shows the opposite effect. Taken together, our results reveal that DNA methylation plays an important role in poultry IMF deposition.

The circular RNA aplacirc_13267 upregulates duck granulosa cell apoptosis by the apla-miR-1-13/THBS1 signaling pathway

Follicle development is a key factor that determines the reproductive performance of poultry. The existing evidence suggests that circular RNAs (circRNAs) play an important role in a variety of biological processes, especially in posttranscriptional regulation, but the regulatory mechanism of circRNAs in duck follicle development has rarely been reported. To better explore the molecular mechanism of follicle development in ducks, we sequenced and analyzed the follicular circRNAs; 4,204 circRNAs were predicted in the duck follicles. Fourteen circRNAs were differentially expressed between the white follicles and yellow follicles. The results of our studies showed that aplacirc_013267 promoted cell apoptosis in duck GCs. Moreover, a bioinformatics prediction analysis demonstrated that aplacirc_013267 was involved in a circRNA-miRNA-mRNA coexpression network and was observed to sponge two follicle-related miRNAs by a luciferase activity assay. Furthermore, we found that overexpression of aplacirc_013267 significantly increased thrombospondin-1 (THBS1) expression and downregulated granulosa cell apoptosis. The mechanistic study showed that aplacirc_013267 directly binds to and inhibits apla-mir-1-13; then, aplacirc_013267 increases the expression of THBS1 and upregulates granulosa cell apoptosis. Taken together, our findings demonstrate that circRNAs have potential effects in duck ovarian follicles and that circRNAs may represent a new avenue to understand follicular development.

Small RNA sequencing of pectoral muscle tissue reveals microRNA-mediated gene modulation in chicken muscle growth

Sichuan mountainous black-bone (SMB) chicken is a small-sized black-feathered chicken breed with low amount of meat, while Dahen (DH) chicken has a larger body size and a faster growth rate. MicroRNAs (miRNAs) are involved in various physiological processes, but their role in chicken muscle growth remains unclear. We aimed to investigate the miRNAs and pathways participating in the muscle growth of chicken. MiRNA profiles of four SMB chickens and four DH chickens were detected by small RNA sequencing. A total of 994 known miRNAs were identified, among which gga-miR-1a-3p, gga-miR-148-3p and gga-miR-133a-3p exhibited the highest enrichment in both breeds of chickens. Thirty-two miRNAs were differently expressed between SMB and DH chickens. The differently expressed miRNAs were mainly associated with fatty acid metabolism, immunity and MAPK activation-related processes. Kyoto encyclopaedia of genes and genomes (KEGG) analysis showed that miRNAs were involved in the immunity-related and MAPK signalling pathways. Moreover, miR-204 was downregulated in DH chicken compared with SMB chicken, and significantly inhibited the expression of MAP3K13, which is involved in the MAPK pathway. It was confirmed through luciferase reporter assays that miR-204 specifically inhibited the activity of MAP3K13. Our results helped demonstrate the potential molecular mechanisms of muscle growth in chickens and provide valuable information for chicken breeding.

Krüppel-like factor KLF9 inhibits chicken intramuscular preadipocyte differentiation

1. Poultry meat quality is affected by many factors, among which intramuscular fat (IMF) is predominant. IMF content affects tenderness, juiciness and flavour of meat. Krüppel-like transcriptional factors (KLFs) are important regulators of adipocyte differentiation. However, little is known about the KLF9 gene associated with poultry IMF deposition, especially intramuscular adipocyte differentiation.2. Previous work has shown that chicken KLF9 was differentially expressed during adipogenesis of intramuscular preadipocytes differentiation. In this study, the function of KLF9 in chicken intramuscular preadipocytes differentiation was investigated.3. In the chicken preadipocyte differentiation model, KLF9 expression showed a major increase with adipogenic induction. Overexpression of KLF9 down-regulated the expression of the adipogenic marker gene AP2, and impaired triglyceride accumulation. Knockdown of KLF9 in chicken intramuscular preadipocytes increased the expression of PPARG, CEBPA and AP2. In addition, it was proposed that KLF9 may regulate adipogenesis via lncRNAs NONGGAT002209.2, NONGGAT003346.2, NONGGAT000436.2 and NONGGAT006302.2 in chicken.4. The data supported a novel role of KLF9 in regulating chicken intramuscular preadipocyte differentiation. Such findings may contribute to a more thorough understanding of chicken IMF deposition and the improvement of poultry meat quality.

MicroRNA-15a Regulates the Differentiation of Intramuscular Preadipocytes by Targeting ACAA1, ACOX1 and SCP2 in Chickens

Our previous studies showed that microRNA-15a (miR-15a) was closely related to intramuscular fat (IMF) deposition in chickens; however, its regulatory mechanism remains unclear. Here, we evaluated the expression characteristics of miR-15a and its relationship with the expression of acetyl-CoA acyltransferase 1 (ACAA1), acyl-CoA oxidase 1 (ACOX1) and sterol carrier protein 2 (SCP2) by qPCR analysis in Gushi chicken breast muscle at 6, 14, 22, and 30 weeks old, where we performed transfection tests of miR-15a mimics in intramuscular preadipocytes and verified the target gene of miR-15a in chicken fibroblasts (DF1). The miR-15a expression level at 30 weeks increased 13.5, 4.5, and 2.7-fold compared with the expression levels at 6, 14, and 22 weeks, respectively. After 6 days of induction, miR-15a over-expression significantly promoted intramuscular adipogenic differentiation and increased cholesterol and triglyceride accumulation in adipocytes. Meanwhile, 48 h after transfection with miR-15a mimics, the expression levels of ACAA1, ACOX1 and SCP2 genes decreased by 56.52%, 31.18% and 37.14% at the mRNA level in intramuscular preadipocytes. In addition, the co-transfection of miR-15a mimics and ACAA1, ACOX1 and SCP2 3′UTR (untranslated region) dual-luciferase vector significantly inhibited dual-luciferase activity in DF1 cells. Taken together, our data demonstrate that miR-15a can reduce fatty acid oxidation by targeting ACAA1, ACOX1, and SCP2, which subsequently indirectly promotes the differentiation of chicken intramuscular preadipocytes.

Acetyl-coenzyme A acyltransferase 2 promote the differentiation of sheep precursor adipocytes into adipocytes

The acetyl CoA acyltransferase 2 (ACAA2) is a key enzyme of the fatty acid oxidation pathway, catalyzing the last step of the mitochondrial beta oxidation, thus playing an important role in the fatty acid metabolism. The purpose of this study was to investigate the effect of knocking out ACAA2 on the expression of genes lipoprteinlipase (LPL), peroxisome proliferator-activated receptor-γ (PPAR-γ), fatty acid synthase, fat mass and obesity-associated gene, adipocyte fatty acid-binding protein (AP2) in precursor adipocytes and their differentiation into adipocytes. The knockout vector was constructed using CRISPR-Cas RNA-guided nuclease technology with an efficiency of 23.80%, and the vector was transfected into precursor adipocyte cells, while an overexpression vector of the ACAA2 gene was also transfected in another group of preadipocytes. Quantitative polymerase chain reaction showed that the expression of the PPAR-γ, LPL, and AP2 was significantly lower in the knockout compared with the overexpression group, while there was no difference in cell growth. After induction of adipocyte precursor cells into adipocytes using dexamethasone, insulin, and IBMX, oil red staining showed a significantly different number of lipid droplets in the knockout group. These results provide a preliminary indication for a possible involvement of the ACAA2 gene in adipocyte differentiation in vitro.

Profiling and characterization of miRNAs associated with intramuscular fat content in Yorkshire pigs

miRNAs are short noncoding RNAs that post-transcriptionally regulate gene expression by binding to complementary regions of the target mRNA. The miRNAs associated with the deposition of intramuscular fat (IMF) content in pigs, which is an important meat quality trait, still remain to be investigated. In this study, the longissimus dorsi muscles (LDMs) from 234 individuals were collected from Yorkshire pigs at 90 kg body weight and the miRNA deep sequencing was conducted by using two tailed groups which were taken five individuals each from high (2.94 ± 0.04%) and low (1.62 ± 0.02%) IMF samples. The results showed that total 268 mature miRNAs were identified, of which 70 were previously known, 162 were conserved among species and 36 were identified specifically in pigs. Moreover, 28 miRNAs involved in adipogenesis were differentially expressed in the two groups, and five out of 16 miRNAs were validated by quantitative PCR (qPCR) using stem loop primers. Our results may serve as a fundamental basis for understanding the roles of miRNA in IMF development in pigs. The miRNAs identified in our study can be utilized for research IMF trait in pig population and will provide further clues to the study of meat quality regulatory mechanisms.

Weighted gene coexpression network analysis identifies specific transcriptional modules and hub genes related to intramuscular fat traits in chicken breast muscle

Intramuscular fat (IMF) traits are important factors that influence meat quality. However, the molecular regulatory mechanisms that underlie this trait in chickens are still poorly understood at the gene coexpression level. Here, we performed a weighted gene coexpression network analysis between IMF traits and transcriptome profile in breast muscle in the Chinese domestic Gushi chicken breed at 6, 14, 22, and 30 weeks. A total of 26 coexpressed gene modules were identified. Six modules, which included the dark gray, purple, cyan, pink, light cyan, and blue modules, showed a significant positive correlation (P < 0.05) with IMF traits. The strongest correlation was observed between the dark gray module and IMF content (r = 0.85; P = 4e-04) and between the blue module and different fatty acid content (r = 0.87~0.91; P = 5e-05~2e-04). Enrichment analysis showed that the enrichment of biological processes, such as fatty acid metabolic process, fat cell differentiation, acylglycerol metabolic process, and glycerolipid metabolism were significantly different in the six modules. In addition, the 32, 24, 4, 7, 6, and 25 hub genes were identified from the blue, pink, light cyan, cyan, dark gray, and purple modules, respectively. These hub genes are involved in multiple links to fatty acid metabolism, phospholipid metabolism, cholesterol metabolism, diverse cellular behaviors, and cell events. These results provide novel insights into the molecular regulatory mechanisms for IMF-related traits in chicken and may also help to uncover the formation mechanism for excellent meat quality traits in local breeds of Chinese chicken.

Analyses of MicroRNA and mRNA Expression Profiles Reveal the Crucial Interaction Networks and Pathways for Regulation of Chicken Breast Muscle Development

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.