Uncovering the embryonic development-related proteome and metabolome signatures in breast muscle and intramuscular fat of fast-and slow-growing chickens

Unraveling the Genetic Foundations of Broiler Meat Quality: Advancements in Research and Their Impact

As societal progress elevates living standards, the focus on meat consumption has shifted from quantity to quality. In broiler production, optimizing meat quality has become paramount, prompting efforts to refine various meat attributes. Recent advancements in sequencing technologies have revealed the genome's complexity, surpassing previous conceptions. Through experimentation, numerous genetic elements have been linked to crucial meat quality traits in broiler chickens. This review synthesizes the current understanding of genetic determinants associated with meat quality attributes in broilers. Researchers have unveiled the pivotal insights detailed herein by employing diverse genomic methodologies such as QTL-based investigations, candidate gene studies, single-nucleotide polymorphism screening, genome-wide association studies, and RNA sequencing. These studies have identified numerous genes involved in broiler meat quality traits, including meat lightness (COL1A2 and ACAA2), meat yellowness (BCMO1 and GDPD5), fiber diameter (myostatin and LncIRS1), meat pH (PRDX4), tenderness (CAPN1), and intramuscular fat content (miR-24-3p and ANXA6). Consequently, a comprehensive exploration of these genetic elements is imperative to devise novel molecular markers and potential targets, promising to revolutionize strategies for enhancing broiler meat quality.

Multiomics analyses of Jining Grey goat and Boer goat reveal genomic regions associated with fatty acid and amino acid metabolism and muscle development

OBJECTIVE

Jining Grey goat is a local Chinese goat breed that is well known for its high fertility and excellent meat quality but shows low meat production performance. Numerous studies have focused on revealing the genetic mechanism of its high fertility, but its highlighting meat quality and muscle growth mechanism still need to be studied.

METHODS

In this research, an integrative analysis of the genomics and transcriptomics of Jining Grey goats compared with Boer goats was performed to identify candidate genes and pathways related to the mechanisms of meat quality and muscle development.

RESULTS

Our results overlap among five genes (ABHD2, FN1, PGM2L1, PRKAG3, RAVER2) and detected a set of candidate genes associated with fatty acid metabolism (PRKAG3, HADHB, FASN, ACADM), amino acid metabolism (KMT2C, PLOD3, NSD2, SETDB1, STT3B, MAN1A2, BCKDHB, NAT8L, P4HA3) and muscle development (MSTN, PPARGC1A, ANKRD2). Several pathways have also been detected, such as the FoxO signaling pathway and Apelin signaling pathway that play roles in lipid metabolism, lysine degradation, N-glycan biosynthesis, valine, leucine and isoleucine degradation that involving with amino acid metabolism.

CONCLUSION

The comparative genomic and transcriptomic analysis of Jining Grey goat and Boer goat revealed the mechanisms underlying the meat quality and meat productive performance of goats. These results provide valuable information for future breeding of goats.

Metabolomic, lipidomic, and proteomic profiles provide insights on meat quality differences between Shitou and Wuzong geese

A comprehensive comparison of metabolomic, lipidomic, and proteomic profiles was conducted between the breast and leg muscles of Shitou goose (STE) and Wuzhong goose (WZE), which exhibit significant variations in body size and growth rate, to evaluate their impact on meat quality. WZE had higher intramuscular fat content in their breast muscles, which were also chewier and had higher drip and cooking losses than STE. Metabolomic analysis revealed differential regulation of amino acid and purine metabolism between WZE and STE. Lipidomic analysis indicated a higher abundance of PE and PC lipid molecules in WZE. Integration of proteomic and metabolomic data highlighted purine metabolism and amino acid biosynthesis as the major distinguishing pathways between STE and WZE. The primary differential pathways between breast and leg muscles were associated with energy metabolism and fatty acid metabolism. This comprehensive analysis provides valuable insights into the distinct meat quality of STE and WZE.

Large-scale genomic and transcriptomic analyses elucidate the genetic basis of high meat yield in chickens

INTRODUCTION

Investigating the genetic markers and genomic signatures related to chicken meat production by combing multi-omics methods could provide new insights into modern chicken breeding technology systems.

OBJECT

Chicken is one of the most efficient and environmentally friendly livestock, especially the fast-growing white-feathered chicken (broiler), which is well known for high meat yield, but the underlying genetic basis is poorly understood.

METHOD

We generated whole-genome resequencing of three purebred broilers (n = 748) and six local breeds/lines (n = 114), and sequencing data of twelve chicken breeds (n = 199) were obtained from the NCBI database. Additionally, transcriptome sequencing of six tissues from two chicken breeds (n = 129) at two developmental stages was performed. A genome-wide association study combined with cis-eQTL mapping and the Mendelian randomization was applied.

RESULT

We identified > 17 million high-quality SNPs, of which 21.74% were newly identified, based on 21 chicken breeds/lines. A total of 163 protein-coding genes underwent positive selection in purebred broilers, and 83 genes were differentially expressed between purebred broilers and local chickens. Notably, muscle development was proven to be the major difference between purebred broilers and local chickens, or ancestors, based on genomic and transcriptomic evidence from multiple tissues and stages. The MYH1 gene family showed the top selection signatures and muscle-specific expression in purebred broilers. Furthermore, we found that the causal gene SOX6 influenced breast muscle yield and also related to myopathy occurrences. A refined haplotype was provided, which had a significant effect on SOX6 expression and phenotypic changes.

CONCLUSION

Our study provides a comprehensive atlas comprising the typical genomic variants and transcriptional characteristics for muscle development and suggests a new regulatory target (SOX6-MYH1s axis) for breast muscle yield and myopathy, which could aid in the development of genome-scale selective breeding aimed at high meat yield in broiler 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.

Identification of key genes related to intramuscular fat deposition in Beijing-You chicken by mRNA and miRNA transcriptome analysis

Intramuscular fat (IMF) is an important factor affecting chicken quality. However, the age-related mechanism of IMF deposition has not yet been elucidated. In this study, the IMF, phospholipids (PL), triglycerides (TG), and fatty acid (FA) content in the breast muscle of Beijing-You chicken (BJY) at 1, 56, 98, and 120 d of age was measured, and mRNA and miRNA sequencing was integrated to explore the regulatory genes of IMF deposition. The results showed that the IMF content of BJY at 1 d of age was significantly higher than that at later stage of birth (P < 0.05). The transcriptome sequencing results showed that 7, 225 differentially expressed genes (DEGs) and 243 differentially expressed miRNAs (DE-miRNAs) were identified. The cluster analysis showed that the expression of DEGs and DE-miRNAs at 1 d of age was significantly different from that at later stages of birth. Furthermore, a potential mRNA-miRNA regulatory network related to IMF deposition was established by weighted gene co-expression network analysis (WGCNA); gga-miR-29c-3p-PIK3R1, gga-miR-6701-3p-PTEN, gga-miR-363-3p-PTEN, gga-miR-1563-WWP1, gga-miR-449c/d-5p-TRAF6, and gga-miR-6701-3p-BMPR1B were identified as key mRNA-miRNA pairs for the regulation of IMF deposition. These results will help elucidate the mechanism of IMF formation mediated by miRNAs in chickens, and provide a theoretical foundation for the genetic improvement of broiler meat quality.

Comparative Analysis of Different Proteins and Metabolites in the Liver and Ovary of Local Breeds of Chicken and Commercial Chickens in the Later Laying Period

The liver and ovary perform a vital role in egg production in hens. In the later laying period, the egg-laying capacity of female hens, particularly that of local breeds, declines significantly. Hence, it is essential to study the features and conditions of the ovary and liver during this period. In this research, we characterized the proteins and metabolites in the liver and ovary of 55-week-old Guangyuan gray chickens (Group G) and Hy-Line gray chickens (Group H) by using liquid chromatography chip/electrospray ionization quadruple time-of-flight/mass spectroscopy (LC-MS/MS). In total, 139 differentially expressed proteins (DEPs) and 186 differential metabolites (DMs) were identified in the liver, and 139 DEPs and 36 DMs were identified in the ovary. The upregulated DEPs and DMs in both the liver and ovary of Group G were primarily enriched in pathways involved in amino acid and carbohydrate metabolism. This suggests that energy metabolism was highly active in the Guangyuan gray chickens. In contrast, the upregulated DEPs and DMs in Group H were mainly enriched in pathways associated with lipid metabolism, which may explain the higher egg production and the higher fatty liver rate in Hy-Line gray hens in the later laying period. Additionally, it was found that the unique protein s-(hydroxymethyl) glutathione dehydrogenase (ADH4) in Group G was implicated in functions such as fatty acid degradation, glycolysis, and pyruvate metabolism, whereas the unique proteins, steroid sulfatase (STS), glucosylceramidase (LOC107050229), and phospholipase A2 Group XV (PLA2G15), in Group H were involved in the metabolism of steroid hormones and glycerol phosphate. In conclusion, variations in how carbohydrates, lipids, and amino acids are processed in the liver and ovary of local breeds of chicken and commercial hens towards the end of their laying period could explain the disparities in their egg production abilities.

Transcriptome-Based Identification of the Muscle Tissue-Specific Expression Gene CKM and Its Regulation of Proliferation, Apoptosis and Differentiation in Chicken Primary Myoblasts

Skeletal muscle is an essential tissue in meat-producing animals, and meat-producing traits have been a hot topic in chicken genetic breeding research. Current research shows that creatine kinase M-type-like (CKM) is one of the most abundant proteins in skeletal muscle and plays an important role in the growth and development of skeletal muscle, but its role in the development of chicken skeletal muscle is still unclear. Via RNA sequencing (RNA-seq), we found that CKM was highly expressed in chicken breast muscle tissue. In this study, the expression profile of CKM was examined by quantitative real-time PCR (qPCR), and overexpression and RNA interference techniques were used to explore the functions of CKM in the proliferation, apoptosis and differentiation of chicken primary myoblasts (CPMs). It was shown that CKM was specifically highly expressed in breast muscle and leg muscle and was highly expressed in stage 16 embryonic muscle, while CKM inhibited proliferation, promoted the apoptosis and differentiation of CPMs and was involved in regulating chicken myogenesis. Transcriptome sequencing was used to identify genes that were differentially expressed in CPMs after CKM disruption, and bioinformatics analysis showed that CKM was involved in regulating chicken myogenesis. In summary, CKM plays an important role in skeletal muscle development during chicken growth and development.

Screening and functional validation of lipid metabolism-related lncRNA-46546 based on the transcriptome analysis of early embryonic muscle tissue in chicken

OBJECTIVE

The study was conducted to screen differentially expressed long noncoding RNA (lncRNA) in chickens by high-throughput sequencing and explore its mechanism of action on intramuscular fat deposition.

METHODS

Herein, Rose crown and Cbb broiler chicken embryo breast and leg muscle lncRNA and mRNA expression profiles were constructed by RNA sequencing. A total of 96 and 42 differentially expressed lncRNAs were obtained in Rose crown vs Cobb broiler chicken breast and leg muscle, respectively. lncRNA-ENSGALT00000046546, with high interspecific variability and a potential regulatory role in lipid metabolism, and its predicted downstream target gene 1-acylglycerol-3-phosphate-O-acyltransferase 2 (AGPAT2), were selected for further study on the preadipocytes.

RESULTS

lncRNA-46546 overexpression in chicken preadipocyte 2 cells significantly increased (p<0.01) the expression levels of AGPAT2 and its downstream genes diacylglycerol acyltransferase 1 and diacylglycerol acyltransferase 2 and those of the fat metabolism-related genes peroxisome proliferator-activated receptor γ, CCAAT/enhancer binding protein α, fatty acid synthase, sterol regulatory element-binding transcription factor 1, and fatty acid binding protein 4. The lipid droplet concentration was higher in the overexpression group than in the control cells, and the triglyceride content in cells and medium was also significantly increased (p<0.01).

CONCLUSION

This study preliminarily concludes that lncRNA-46546 may promote intramuscular fat deposition in chickens, laying a foundation for the study of lncRNAs in chicken early embryonic development and fat deposition.

CircUBE3A promotes myoblasts proliferation and differentiation by sponging miR-28-5p to enhance expression

circRNAs have been found to be involved in the regulatory network of skeletal muscle development in studies. However, their precise functions and regulatory mechanisms remain unknown. The expression patterns and alterations of circRNAs in the longissimus dorsi muscle of two major developmental stages of goats (D75 fetus and D1 kid) were studied using high-throughput sequencing technology and bioinformatics tools in this study. In kid skeletal muscles, 831 differently expressed circRNAs were found, comprising 486 up-regulated circRNAs and 345 down-regulated circRNAs. In skeletal muscle, we focused on the highly expressed and variably expressed circUBE3A. CircUBE3A levels were discovered to be much higher in kid skeletal muscle and differentiated myoblasts. Knocking down circUBE3A resulted in a significant increase in cell proliferation and differentiation in goat myoblasts. CircUBE3A specifically binds to and inhibits miR-28-5p, boosting the expression of Hydroxyacyl Coenzyme A Dehydrogenase Beta (HADHB) and contributing to goat myoblast proliferation and differentiation, according to the mechanistic investigation. The above results indicated that circUBE3A could regulate HADHB expression by targeting miR-28-5p, consequently increasing goat myoblast proliferation and differentiation. Our findings offer fresh perspectives on goat breeding and growth regulation, as well as substantial theoretical basis.

Assessing Population Structure and Signatures of Selection in Wanbei Pigs Using Whole Genome Resequencing Data

Wanbei pig (WBP) is one of the indigenous pig resources in China and has many germplasm characteristics. However, research on its genome is lacking. To assess the genomic variation, population structure, and selection signatures, we resequenced 18 WBP for the first time and performed a comprehensive analysis with resequenced data of 10 Asian wild boars. In total, 590.03 Gb of data and approximately 41 million variants were obtained. Polymorphism level (θπ) ratio and genetic differentiation (fixation index)-based cross approaches were applied, and 539 regions, which harbored 176 genes, were selected. Functional analysis of the selected genes revealed that they were associated with lipid metabolism (SCP2, APOA1, APOA4, APOC3, CD36, BCL6, ADCY8), backfat thickness (PLAG1, CACNA2D1), muscle (MYOG), and reproduction (CABS1). Overall, our results provide a valuable resource for characterizing the uniqueness of WBP and a basis for future breeding.

Proteomic analysis of the liver regulating lipid metabolism in Chaohu ducks using two-dimensional electrophoresis

In this study, we aimed to characterize the liver protein profile of Chaohu ducks using two-dimensional electrophoresis and proteomics. The livers were quickly collected from 120 healthy, 84-day-old Chaohu ducks. The intramuscular fat (IMF) content of the left pectoralis muscle was determined using the Soxhlet extraction method. The total protein of liver tissues from the high and low IMF groups was extracted for proteomics. Functional enrichment analysis of the differentially expressed proteins (DEPs) was conducted using gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG). In total, 43 DEPs were identified. Functional enrichment analysis indicated that these DEPs were significantly related to four lipid metabolic processes: carboxylic acid metabolic process, ATP metabolic process, oxoacid metabolic process, and organic acid metabolic process. Three pathways correlated with lipid metabolism were identified using KEGG analysis: glycolysis/gluconeogenesis, pentose phosphate pathway, fructose, and mannose metabolism. Eight key proteins associated with lipid metabolism were identified: ALDOB, GAPDH, ENO1, RGN, TPI1, HSPA9, PRDX1, and GPX1. Protein–protein interaction analysis revealed that the glycolysis/gluconeogenesis pathway mediated the interaction relationship. Key proteins and metabolic pathways were closely related to lipid metabolism and showed a strong interaction in Chaohu ducks.

Breast Meat Fatty Acid Profiling and Proteomic Analysis of Beijing-You Chicken During the Laying Period

The disparity in fatty acids (FA) composition exhibits a significant impact on meat quality, however, the molecular regulatory mechanisms underlying this trait in chicken are far from clear. In this study, a total of 45 female Beijing-You chicken (BYC) hens, fed on the same diet, were collected at the slaughter age of 150, 300, or 450 days (D150, D300, and D450) from sexual maturation stage to culling stage (15 birds per age). Gas chromatography-mass spectrometry (GC-MS) and tandem mass tag labeling technology based on liquid chromatography mass spectrometry (TMT-LC-MS/MS) analysis strategies were applied to profile FA compositions and to compare differential expressed proteins (DEPs) between these different slaughter ages, respectively. The FA profiling showed that increasing hen ages resulted in increased contents of both saturated and unsaturated fatty acids. Proteomic analyses showed a total of 4,935 proteins in chicken breast muscle with the false discovery rate (FDR) < 1% and 664 of them were differentially expressed (fold change > 1.50 or < 0.67 and P < 0.01). There were 410 up- and 116 down-regulated proteins in D150 vs. D300 group, 32 up- and 20 down-regulated in D150 vs. D450 group, and 72 up- and 241 down-regulated in D300 vs. D450 group. A total of 57 DEPs related to FA/lipid-related metabolisms were obtained according to the enrichment analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). These DEPs were involved in 21 significantly enriched (P < 0.05) pathways, including well-known pathways for FA synthesis (metabolism, desaturation, and elongation) and the signaling pathways for lipid metabolism (PPAR, adipocytokine, calcium, VEGF, MAPK, and Wnt). In addition, there existed several representative DEPs (FABP, FABP3, apoA-I, apoA-IV, apoC-III, apoB, VTG1, and VTG2) involved in the regulation of FA/lipid transportation. The construction of the interaction networks indicated that HADH, ACAA2, HADHA, ACSL1, CD36, CPT1A, PPP3R1, and SPHK1 were the key core nodes. Finally, eight DEPs were quantified using parallel reaction monitoring (PRM) to validate the results from TMT analysis. These results expanded our understanding of how the laying age affects the FA compositions and metabolism in hen breast meat.

Integrative analysis of circRNA, miRNA, and mRNA profiles to reveal ceRNA regulation in chicken muscle development from the embryonic to post-hatching periods

Background

The growth and development of skeletal muscle are regulated by protein-coding genes and non-coding RNA. Circular RNA (circRNA) is a type of non-coding RNA involved in a variety of biological processes, especially in post-transcriptional regulation. To better understand the regulatory mechanism of circRNAs during the development of muscle in chicken, we performed RNA-seq with linear RNA depletion for chicken breast muscle in 12 (E 12) and17 (E 17) day embryos, and 1 (D 1), 14 (D 14), 56 (D 56), and 98 (D 98) days post-hatch.

Results

We identified 5755 differentially expressed (DE)-circRNAs during muscle development. We profiled the expression of DE-circRNAs and mRNAs (identified in our previous study) at up to six time points during chicken muscle development and uncovered a significant profile (profile 16) for circRNA upregulation during aging in muscle tissues. To investigate competing endogenous RNA (ceRNA) regulation in muscle and identify muscle-related circRNAs, we constructed a circRNA-miRNA-mRNA regulatory network using the circRNAs and mRNAs from profile 16 and miRNAs identified in our previous study, which included 361 miRNAs, 68 circRNAs, 599 mRNAs, and 31,063 interacting pairs. Functional annotation showed that upregulated circRNAs might contribute to glycolysis/gluconeogenesis, biosynthesis of amino acids, pyruvate metabolism, carbon metabolism, glycogen and sucrose metabolism through the ceRNA network, and thus affected postnatal muscle development by regulating muscle protein deposition. Of them, circRNA225 and circRNA226 from the same host gene might be key circRNAs that could regulate muscle development by interacting with seven common miRNAs and 207 mRNAs. Our experiments also demonstrated that there were interactions among circRNA225, gga-miR-1306-5p, and heat shock protein alpha 8 (HSPA8).

Conclusions

Our results suggest that adequate supply of nutrients such as energy and protein after hatching may be a key factor in ensuring chicken yield, and provide several candidate circRNAs for future studies concerning ceRNA regulation during chicken muscle development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08525-5.

Chromatin Interaction Responds to Breast Muscle Development and Intramuscular Fat Deposition Between Chinese Indigenous Chicken and Fast-Growing Broiler

Skeletal muscle development and intramuscular fat (IMF) content, which positively contribute to meat production and quality, are regulated by precisely orchestrated processes. However, changes in three-dimensional chromatin structure and interaction, a newly emerged mediator of gene expression, during the skeletal muscle development and IMF deposition have remained unclear. In the present study, we analyzed the differences in muscle development and IMF content between one-day-old commercial Arbor Acres broiler (AA) and Chinese indigenous Lushi blue-shelled-egg chicken (LS) and performed Hi-C analysis on their breast muscles. Our results indicated that significantly higher IMF content, however remarkably lower muscle fiber diameter was detected in breast muscle of LS chicken compared to that of AA broiler. The chromatin intra-interaction was prior to inter-interaction in both AA and LS chicken, and chromatin inter-interaction was heavily focused on the small and gene-rich chromosomes. For genomic compartmentalization, no significant difference in the number of B type compartments was found, but AA had more A type compartments versus LS. The A/B compartment switching of AA versus LS showed more A to B switching than B to A switching. There were no significant differences in the average sizes and distributions of topologically associating domains (TAD). Additionally, approximately 50% of TAD boundaries were overlapping. The reforming and disappearing events of TAD boundaries were identified between AA and LS chicken breast muscles. Among these, the HMGCR gene was located in the TAD-boundary regions in AA broilers, but in TAD-interior regions in LS chickens, and the IGF2BP3 gene was located in the AA-unique TAD boundaries. Both HMGCR and IGF2BP3 genes exhibited increased mRNA expression in one-day-old AA broiler breast muscles. It was demonstrated that the IGF2BP3 and HMGCR genes regulated by TAD boundary sliding were potential biomarkers for chicken breast muscle development and IMF deposition. Our data not only provide a valuable understanding of higher-order chromatin dynamics during muscle development and lipid accumulation but also reveal new insights into the regulatory mechanisms of muscle development and IMF deposition in chicken.

Transcriptome Analysis Reveals the Profile of Long Non-coding RNAs During Chicken Muscle Development

The developmental complexity of muscle arises from elaborate gene regulation. Long non-coding RNAs (lncRNAs) play critical roles in muscle development through the regulation of transcription and post-transcriptional gene expression. In chickens, previous studies have focused on the lncRNA profile during the embryonic periods, but there are no studies that explore the profile from the embryonic to post-hatching period. Here, we reconstructed 14,793 lncRNA transcripts and identified 2,858 differentially expressed lncRNA transcripts and 4,282 mRNAs from 12-day embryos (E12), 17-day embryos (E17), 1-day post-hatch chicks (D1), 14-day post-hatch chicks (D14), 56-day post-hatch chicks (D56), and 98-day post-hatch chicks (D98), based on our published RNA-seq datasets. We performed co-expression analysis for the differentially expressed lncRNAs and mRNAs, using STEM, and identified two profiles with opposite expression trends: profile 4 with a downregulated pattern and profile 21 with an upregulated pattern. The cis- and trans-regulatory interactions between the lncRNAs and mRNAs were predicted within each profile. Functional analysis of the lncRNA targets showed that lncRNAs in profile 4 contributed to the cell proliferation process, while lncRNAs in profile 21 were mainly involved in metabolism. Our work highlights the lncRNA profiles involved in the development of chicken breast muscle and provides a foundation for further experiments on the role of lncRNAs in the regulation of muscle development.

Deciphering the miRNA transcriptome of breast muscle from the embryonic to post-hatching periods in chickens

Background

miRNAs play critical roles in growth and development. Various studies of chicken muscle development have focused on identifying miRNAs that are important for embryo or adult muscle development. However, little is known about the role of miRNAs in the whole muscle development process from embryonic to post-hatching periods. Here, we present a comprehensive investigation of miRNA transcriptomes at 12-day embryo (E12), E17, and day 1 (D1), D14, D56 and D98 post-hatching stages.

Results

We identified 337 differentially expressed miRNAs (DE-miRNAs) during muscle development. A Short Time-Series Expression Miner analysis identified two significantly different expression profiles. Profile 4 with downregulated pattern contained 106 DE-miRNAs, while profile 21 with upregulated pattern contained 44 DE-miRNAs. The DE-miRNAs with the upregulated pattern mainly played regulatory roles in cellular turnover, such as pyrimidine metabolism, DNA replication, and cell cycle, whereas DE-miRNAs with the downregulated pattern directly or indirectly contributed to protein turnover metabolism such as glycolysis/gluconeogenesis, pyruvate metabolism and biosynthesis of amino acids.

Conclusions

The main functional miRNAs during chicken muscle development differ between embryonic and post-hatching stages. miRNAs with an upregulated pattern were mainly involved in cellular turnover, while miRNAs with a downregulated pattern mainly played a regulatory role in protein turnover metabolism. These findings enrich information about the regulatory mechanisms involved in muscle development at the miRNA expression level, and provide several candidates for future studies concerning miRNA-target function in regulation of chicken muscle development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07374-y.

Identification of diverse cell populations in skeletal muscles and biomarkers for intramuscular fat of chicken by single-cell RNA sequencing

BACKGROUND

The development of skeletal muscle is closely related to the efficiency of meat production and meat quality. Chicken skeletal muscle development depends on myogenesis and adipogenesis and occurs in two phases-hyperplasia and hypertrophy. However, cell profiles corresponding to the two-phase muscle development have yet to be determined. Single-cell RNA-sequencing (scRNA-seq) can elucidate the cell subpopulations in tissue and capture the gene expression of individual cells, which can provide new insights into the myogenesis and intramuscular adipogenesis.

RESULTS

Ten cell clusters at the post-hatching developmental stage at Day 5 and seven cell clusters at the late developmental stage at Day 100 were identified in chicken breast muscles by scRNA-seq. Five myocyte-related clusters and two adipocyte clusters were identified at Day 5, and one myocyte cluster and one adipocyte cluster were identified at Day 100. The pattern of cell clustering varied between the two stages. The cell clusters showed clear boundaries at the terminal differentiation stage at Day 100; by contrast, cell differentiation was not complete at Day 5. APOA1 and COL1A1 were selected from up-regulated genes in the adipocyte cluster and found to be co-expressed with the ADIPOQ adipocyte marker gene in breast muscles by RNA in situ hybridization.

CONCLUSIONS

This study is the first to describe the heterogeneity of chicken skeletal muscle at two developmental stages. The genes APOA1 and COL1A1 were identified as biomarkers for chicken intramuscular fat cells.

The Role of Incubation Conditions in the Onset of Avian Myopathies

White striping, wooden breast, and spaghetti muscle have become common myopathies in broilers worldwide. Several research reports have indicated that the origin of these lesions is metabolic disorders. These failures in normal metabolism can start very early in life, and suboptimal incubation conditions may trigger some of the key alterations on muscle metabolism. Incubation conditions affect the development of muscle and can be associated with the onset of myopathies. A series of experiments conducted with broilers, turkeys, and ducks are discussed to overview primary information showing the main changes in breast muscle histomorphology, metabolism, and physiology caused by suboptimal incubation conditions. These modifications may be associated with current myopathies. Those effects of incubation on myopathy occurrence and severity have also been confirmed at slaughter age. The impact of egg storage, temperature profiles, oxygen concentrations, and time of hatch have been evaluated. The effects have been observed in diverse species, genetic lines, and both genders. Histological and muscle evaluations have detected that myopathies could be induced by extended hypoxia and high temperatures, and those effects depend on the genetic line. Thus, these modifications in muscle metabolic responses may make hatchlings more susceptible to develop myopathies during grow out due to thermal stress, high-density diets, and fast growth rates.

RNA-Seq Analysis Reveals Hub Genes Involved in Chicken Intramuscular Fat and Abdominal Fat Deposition During Development

Fat traits are important in the chicken industry where there is a desire for high intramuscular fat (IMF) and low abdominal fat. However, there is limited knowledge on the relationship between the dynamic status of gene expression and the body fat deposition in chicken. Transcriptome data were obtained from breast muscle and abdominal fat of female chickens from nine developmental stages (from embryonic day 12 to hatched day 180). In total, 8,545 genes in breast muscle and 6,824 genes in abdominal fat were identified as developmentally dynamic genes. Weighted correlation network analysis was used to identify gene modules and the hub genes. Twenty-one hub genes were identified, e.g., ENSGALG00000041996, which represents a candidate for high IMF, and CREB3L1, which relates to low abdominal fat weight. The transcript factor L3MBTL1 and the transcript factor cofactors TNIP1, HAT1, and BEND6 related to both high breast muscle IMF and low abdominal fat weight. Our results provide a resource of developmental transcriptome profiles in chicken breast muscle and abdominal fat. The candidate genes can be used in the selection for increased IMF content and/or a decrease in abdominal fat weight which would contribute to the improvement of these traits.

The whole-transcriptome landscape of muscle and adipose tissues reveals the ceRNA regulation network related to intramuscular fat deposition in yak

Background

The Intramuscular fat (IMF) content in meat products, which is positively correlated with meat quality, is an important trait considered by consumers. The regulation of IMF deposition is species specific. However, the IMF-deposition-related mRNA and non-coding RNA and their regulatory network in yak (Bos grunniens) remain unknown. High-throughput sequencing technology provides a powerful approach for analyzing the association between transcriptome-related differences and specific traits in animals. Thus, the whole transcriptomes of yak muscle and adipose tissues were screened and analyzed to elucidate the IMF deposition-related genes. The muscle tissues were used for IMF content measurements.

Results

Significant differences were observed between the 0.5- and 2.5-year-old yaks. Several mRNAs, miRNAs, lncRNAs and circRNAs were generally expressed in both muscle and adipose tissues. Between the 0.5- and 2.5-year-old yaks, 149 mRNAs, 62 miRNAs, 4 lncRNAs, and 223 circRNAs were differentially expressed in muscle tissue, and 72 mRNAs, 15 miRNAs, 9 lncRNAs, and 211 circRNAs were differentially expressed in adipose tissue. KEGG annotation revelved that these differentially expressed genes were related to pathways that maintain normal biological functions of muscle and adipose tissues. Moreover, 16 mRNAs, 5 miRNAs, 3 lncRNAs, and 5 circRNAs were co-differentially expressed in both types of tissue. We suspected that these co-differentially expressed genes were involved in IMF-deposition in the yak. Additionally, LPL, ACADL, SCD, and FASN, which were previously shown to be associated with the IMF content, were identified in the competing endogenous RNA (ceRNA) regulatory network that was constructed on the basis of the IMF deposition-related genes. Three ceRNA subnetworks also revealed that TCONS-00016416 and its target SIRT1 “talk” to each other through the same miR-381-y and miR-208 response elements, whereas TCONS-00061798 and its target PRKCA, and TCONS-00084092 and its target LPL “talk” to each other through miR-122-x and miR-499-y response elements, respectively.

Conclusion

Taken together, our results reveal the potential mRNA and noncoding RNAs involved in IMF deposition in the yak, providing a useful resource for further research on IMF deposition in this animal species.

Association Analysis of Single-Cell RNA Sequencing and Proteomics Reveals a Vital Role of Ca2+ Signaling in the Determination of Skeletal Muscle Development Potential

This study is aimed at exploring the mechanism underlying the homeostasis between myogenesis and adipogenesis in skeletal muscle using a special porcine model with a distinct phenotype on muscle growth rate and intramuscular fat deposition. Differentiation potential of muscle-derived Myo-lineage cells of lean-type pigs was significantly enhanced relative to obese-type pigs, while that of their Adi-lineage cells was similar. Single-cell RNA sequencing revealed that lean-type pigs reserved a higher proportion of Myo-lineage cells in skeletal muscle relative to obese-type pigs. Besides, Myo-lineage cells of the lean-type pig settled closer to the original stage of muscle-derived progenitor cells. Proteomics analysis found that differentially expressed proteins between two sources of Myo-lineage cells are mainly involved in muscle development, cell proliferation and differentiation, ion homeostasis, apoptosis, and the MAPK signaling pathway. The regulation of intracellular ion homeostasis, Ca²⁺ in particular, significantly differed between two sources of Myo-lineage cells. Ca²⁺ concentration in both cytoplasm and endoplasmic reticulum was lower in Myo-lineage cells of lean-type pigs relative to obese-type pigs. In conclusion, a higher proportion and stronger differentiation capacity of Myo-lineage cells are the main causes for the higher capability of myogenic differentiation and lower intramuscular fat deposition. Relative low concentration of cellular Ca²⁺ is advantageous for Myo-lineage cells to keep a potent differentiation potential.

Dynamic accumulation of fatty acids in duck (Anas platyrhynchos) breast muscle and its correlations with gene expression

Background

Fatty acid composition contributes greatly to the quality and nutritional value of meat. However, the molecular regulatory mechanisms underlying fatty acid accumulation in poultry have not yet been cleared. The aims of this study were to characterize the dynamics of fatty acid accumulation in duck breast muscle and investigate its correlations with gene expression.

Results

Here, we analyzed the fatty acid profile and transcriptome of breast muscle derived from Pekin ducks and mallards at the ages of 2 weeks, 4 weeks, 6 weeks and 8 weeks. Twenty fatty acids were detected in duck breast muscle, with palmitic acid (C16:0, 16.6%~ 21.1%), stearic acid (C18:0, 9.8%~ 17.7%), oleic acid (C18:1n-9, 15.7%~ 33.8%), linoleic acid (C18:2n-6, 10.8%~ 18.9%) and arachidonic acid (C20:4n-6, 11.7%~ 28.9%) as the major fatty acids. Our results showed that fatty acid composition was similar between the two breeds before 6 weeks, but the compositions diverged greatly after this point, mainly due to the stronger capacity for C16:0 and C18:1n-9 deposition in Pekin ducks. By comparing the multistage transcriptomes of Pekin ducks and mallards, we identified 2025 differentially expressed genes (DEGs). Cluster analysis of these DEGs revealed that the genes involved in oxidative phosphorylation, fatty acid degradation and the PPAR signaling pathway were upregulated in mallard at 8 weeks. Moreover, correlation analysis of the DEGs and fatty acid composition traits suggested that the DEGs involved in lipogenesis, lipolysis and fatty acid β-oxidation may interact to influence the deposition of fatty acids in duck breast muscle.

Conclusions

We reported the temporal progression of fatty acid accumulation and the dynamics of the transcriptome in breast muscle of Pekin ducks and mallards. Our results provide insights into the transcriptome regulation of fatty acid accumulation in duck breast muscle, and will facilitate improvements of fatty acid composition in duck breeding.

Dynamic Transcriptomic Analysis of Breast Muscle Development From the Embryonic to Post-hatching Periods in Chickens

Skeletal muscle development and growth are closely associated with efficiency of poultry meat production and its quality. We performed whole transcriptome profiling based on RNA sequencing of breast muscle tissue obtained from Shouguang chickens at embryonic days (E) 12 and 17 to post-hatching days (D) 1, 14, 56, and 98. A total of 9,447 differentially expressed genes (DEGs) were filtered (Q < 0.01, fold change > 2). Time series expression profile clustering analysis identified five significantly different expression profiles that were divided into three clusters. DEGs from cluster I with downregulated pattern were significantly enriched in cell proliferation processes such as cell cycle, mitotic nuclear division, and DNA replication. DEGs from cluster II with upregulated pattern were significantly enriched in metabolic processes such as glycolysis/gluconeogenesis, insulin signaling pathway, calcium signaling pathway, and biosynthesis of amino acids. DEGs from cluster III, with a pattern that increased from E17 to D1 and then decreased from D1 to D14, mainly contributed to lipid metabolism. Therefore, this study may help us explain the mechanisms underlying the phenotype that myofiber hyperplasia occurs predominantly during embryogenesis and hypertrophy occurs mainly after birth at the transcriptional level. Moreover, lipid metabolism may contribute to the early muscle development and growth. These findings add to our knowledge of muscle development in chickens.

Transcriptional insights into key genes and pathways controlling muscle lipid metabolism in broiler chickens

Background

Intramuscular fat (IMF) is one of the most important factors positively associated with meat quality. Triglycerides (TGs), as the main component of IMF, play an essential role in muscle lipid metabolism. This transcriptome analysis of pectoralis muscle tissue aimed to identify functional genes and biological pathways likely contributing to the extreme differences in the TG content of broiler chickens.

Results

The study included Jingxing-Huang broilers that were significantly different in TG content (5.81 mg/g and 2.26 mg/g, p < 0.01) and deposition of cholesterol also showed the same trend. This RNA sequencing analysis was performed on pectoralis muscle samples from the higher TG content group (HTG) and the lower TG content group (LTG) chickens. A total of 1200 differentially expressed genes (DEGs) were identified between two groups, of which 59 DEGs were related to TG and steroid metabolism. The HTG chickens overexpressed numerous genes related to adipogenesis and lipogenesis in pectoralis muscle tissue, including the key genes ADIPOQ, CD36, FABP4, FABP5, LPL, SCD, PLIN1, CIDEC and PPARG, as well as genes related to steroid biosynthesis (DHCR24, LSS, MSMO1, NSDHL and CH25H). Additionally, key pathways related to lipid storage and metabolism (the steroid biosynthesis and peroxisome proliferator activated receptor (PPAR) signaling pathway) may be the key pathways regulating differential lipid deposition between HTG group and LTG group.

Conclusions

This study showed that increased TG deposition accompanying an increase in steroid synthesis in pectoralis muscle tissue. Our findings of changes in gene expression of steroid biosynthesis and PPAR signaling pathway in HTG and LTG chickens provide insight into genetic mechanisms involved in different lipid deposition patterns in pectoralis muscle tissue.

A new chicken 55K SNP genotyping array

BACKGROUND

China has the richest local chicken breeding resources in the world and is the world's second largest producer of meat-type chickens. Development of a moderate-density SNP array for genetic analysis of chickens and breeding of meat-type chickens taking utility of those resources is urgently needed for conventional farms, breeding industry, and research areas.

RESULTS

Eight representative local breeds or commercial broiler lines with 3 pools of 48 individuals within each breed/line were sequenced and supplied the major SNPs resource. There were 7.09 million - 9.41 million SNPs detected in each breed/line. After filtering using multiple criteria such as preferred incorporation of trait-related SNPs and uniformity of distribution across the genome, 52.18 K SNPs were selected in the final array. It consists of: (i) 19.22 K SNPs from the genomes of yellow-feathered, cyan-shank partridge and white-feathered chickens; (ii) 5.98 K SNPs related to economic traits from the Illumina 60 K SNP Bead Chip, which were found as significant associated SNPs with 15 traits in a Beijing-You crossed Cobb F2 resource population by genome-wide association study analysis; (iii) 7.63 K SNPs from 861 candidate genes of economic traits; (iv) the 0.94 K SNPs related to residual feed intake; and (v) 18.41 K from chicken SNPdb. The polymorphisms of 9 extra local breeds and 3 commercial lines were examined with this array, and 40 K - 47 K SNPs were polymorphic (with minor allele frequency > 0.05) in those breeds. The MDS result showed that those breeds can be clearly distinguished by this newly developed genotyping array.

CONCLUSIONS

We successfully developed a 55K genotyping array by using SNPs segregated from typical local breeds and commercial lines. Compared to the existing Affy 600 K and Illumina 60 K arrays, there were 21,41 K new SNPs included on our Affy 55K array. The results of the 55K genotyping data can therefore be imputed to high-density SNPs genotyping data. The array offers a wide range of potential applications such as genomic selection breeding, GWAS of interested traits, and investigation of diversity of different chicken breeds.

Intramuscular preadipocytes impede differentiation and promote lipid deposition of muscle satellite cells in chickens

Background

Skeletal muscle satellite cells (MSC) are crucial for postnatal growth and regeneration of skeletal muscle. An interaction exists between MSC and intramuscular preadipocytes (IMPA). This study is the first to investigate the effects of IMPA on MSC in chickens and unveil the molecular mechanisms by transcriptome analysis.

Results

Primary MSC and IMPA were isolated from the pectoralis major muscle of 7-day-old chickens. After both cell types reached confluence, MSC were cultured alone or co-cultured with IMPA for 2 or 4 d. MSC treated for 2 d were subjected to RNA-seq. A total of 1653 known differentially expressed genes (DEG) were identified between co-cultured and mono-cultured MSC (|log2 FC| ≥ 1, FDR < 0.01). Based on Gene Ontology analysis, 48 DEG related to muscle development were screened, including the key genes MYOD1, MYOG, PAX7, and TMEM8C. The 44 DEG related to lipid deposition included the key genes CD36, FABP4, ACSBG2, CYP7A1 and PLIN2. Most of the DEG related to muscle development were downregulated in co-cultured MSC, and DEG related to lipid deposition were upregulated. Immunofluorescence of MHC supported IMPA impeding differentiation of MSC, and Oil Red O staining showed concurrent promotion of lipid deposition. Pathway analysis found that several key genes were enriched in JNK/MAPK and PPAR signaling, which may be the key pathways regulating differentiation and lipid deposition in MSC. Additionally, pathways related to cell junctions may also contribute to the effect of IMPA on MSC.

Conclusions

The present study showed that IMPA impeded differentiation of MSC while promoting their lipid deposition. Pathway analysis indicated that IMPA might inhibit differentiation via the JNK/MAPK pathway, and promote lipid deposition via the PPAR pathway. This study supplies insights into the effect of IMPA on MSC, providing new clues on exposing the molecular mechanisms underlying the interplay between skeletal muscle and intramuscular fat in chickens.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5209-5) contains supplementary material, which is available to authorized users.

miRNA-223 targets the GPAM gene and regulates the differentiation of intramuscular adipocytes

Intramuscular fat (IMF) has significant effects on the tenderness, juiciness, and flavor of chicken, which are important determinants of poultry meat quality. Although many studies have focused on microRNAs (miRNAs) involved in adipogenesis, little is known about miRNAs associated with poultry IMF deposition or intramuscular adipocyte differentiation. Bioinformatic analysis identified mitochondrial glycerol‑3‑phosphate acyltransferase (GPAM) as a putative target of miR-223. To explore the role of miR-223 in the process of chicken IMF deposition, loss and gain of function experiments were performed in primary intramuscular preadipocytes using miR-223 mimics, miR-223 inhibitor, and si-GPAM. Our results showed that miR-223 is significantly down-regulated in the breast muscle tissues of Gushi hens at the later-laying period in comparison with hens at the pre-laying period. Using qRT-PCR, we found that miR-223 expression in chicken tissues and intramuscular adipocytes correlates negatively with GPAM expression. Cell transfection experiments suggest that miR-223 inhibits intramuscular adipocyte differentiation via targeting GPAM. Experiments using a dual luciferase reporter system show that GPAM is a direct target of miR-223. Taken together, our results support the hypothesis that miR-223 regulates intramuscular fat deposition in chickens.