Skeletal Muscle Transcriptome Analysis of Hanzhong Ma Duck at Different Growth Stages Using RNA-Seq

Full-Length Transcriptome Analysis of Skeletal Muscle of Jiangquan Black Pig at Different Developmental Stages

Skeletal muscle grows in response to a combination of genetic and environmental factors, and its growth and development influence the quality of pork. Elucidating the molecular mechanisms regulating the growth and development of skeletal muscle is of great significance to both animal husbandry and farm management. The Jiangquan black pig is an excellent pig breed based on the original Yimeng black pig, importing the genes of the Duroc pig for meat traits, and cultivated through years of scientific selection and breeding. In this study, full-length transcriptome sequencing was performed on three growth stages of Jiangquan black pigs, aiming to study the developmental changes in Jiangquan black pigs at different developmental stages at the molecular level and to screen the key genes affecting the growth of skeletal muscle in Jiangquan black pigs. We performed an enrichment analysis of genes showing differential expression and constructed a protein-protein interaction network with the aim of identifying core genes involved in the development of Jiangquan black pigs. Notably, genes such as TNNI2, TMOD4, PLDIM3, MYOZ1, and MYH1 may be potential regulators of muscle development in Jiangquan black pigs. Our results contribute to the understanding of the molecular mechanisms of skeletal muscle development in this pig breed, which will facilitate molecular breeding efforts and the development of pig breeds to meet the needs of the livestock industry.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Analysis of circRNA-miRNA-mRNA regulatory network of embryonic gonadal development in Mulard duck

The aim of the study was to explore the regulatory mechanism of differences in embryonic gonadal development between intergeneric distance hybrid offspring Mulard ducks and parent ducks. The morphological differences gonadal tissues of Muscovy ducks, Pekin ducks and Mulard ducks at 12.5-day embryonic age were observed by sectioning and hematoxylin-eosin (HE) staining. Then followed by transcriptome sequencing to screen for gonadal development-related differentially expressed circRNAs and mRNAs to construct a competitive endogenous RNA (ceRNA) regulatory network. Finally, qRT-PCR and luciferase reporter system were used to verify the sequencing data and targeting relationship of ceRNA pairs. The results showed that the seminiferous tubule lumen of Mulard ducks was not obvious, while there were obvious seminiferous tubules and tubular structures in testis of Pekin ducks and Muscovy ducks, with number and shape indicating maturity. There were 18 upregulated circRNAs and 16 downregulated circRNAs in Mulard ducks and Pekin ducks, respectively, and 39 upregulated circRNAs and 1 downregulated circRNA in Mulard ducks and Muscovy ducks, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis found that genes involves in dorso-ventral axis formation, for example, neurogenic locus notch homolog protein 1 (NOTCH1), were significantly enriched (P < 0.05). The novel_circ_0002265-gga-miR-122-5p-PAFAH1B2 regulatory network was constructed. The qRT-PCR results showed that the sequencing results were reliable. The dual-luciferase reporter assay showed that gga-miR-122-5p exists binding site of circ_0002265 and PAFAH1B2, indicating circ_0002265-gga-miR-122-5p-PAFAH1B2 targeting relationship. In summary, the embryonic gonadal development of intergeneric hybrid Mulard ducks may be regulated by differentially expressed circRNAs and genes, such as novel_circ_0000519, novel_circ_0003537, NOTCH1, FGFR2, PAFAH1B1, and PAFAH1B2, among which circ_0002265-gga-miR-122-5p-PAFAH1B2 may participate in the targeted regulation of gonadal development in Mulard ducks. The findings of this study are helpful for analyzing the mechanism of embryonic gonadal development differences in avians.

Characterization and comparative transcriptomic analysis of skeletal muscle in female Pekin duck and Hanzhong Ma duck during different growth stages using RNA-seq

Duck is an economically important poultry, and there is currently a major focus on improving its meat quality through breeding. There are wide variations in the growth regulation mechanisms of different duck breeds, that fundamental research on skeletal muscle growth is essential for understanding the regulation of unknown genes. The study aimed to broaden the understanding the duck skeletal muscle development and thereby to improve the performance of domestic ducks. In this study, RNA-seq data from skeletal muscles (breast muscle and leg muscle) of Pekin duck and Hanzhong Ma duck sampled at d 17, 21, and 27 of embryo (E17d, E21d, and E27d), as well as at 6-mo-old following birth (M6), to investigate and compare the mRNA temporal expression profiles and associated pathways that regulate skeletal myogenesis of different duck breeds. There were 331 to 1,440 annotated differentially expressed genes (DEGs) in breast muscle and 380 to 1,790 annotated DEGs in leg muscle from different databases between 2 duck breeds. Gene ontology (GO) enrichment in skeletal muscles indicated that these DEGs were mainly involved in biosynthetic process, developmental process, regulation of protein metabolic process and regulation of gene expression. KEGG analysis in skeletal muscles showed that a total of 41 DEGs were mapped to 7 KEGG pathways, including ECM-receptor interaction, focal adhesion, carbon metabolism, regulation of actin cytoskeleton, calcium signaling pathway, biosynthesis of amino acids and PPAR signaling pathway. The differential expression of 8 selected DEGs was verified by qRT-PCR, and the results were consistent with RNA-seq data. The identified DEGs, such as SDC, SPP1, PAK1, MYL9, PGK1, NOS1, PHGDH, TNNT2, FN1, and AQP4, were specially highlighted, indicating their associations with muscle development in the Pekin duck and Hanzhong Ma duck. This study provides a basis for revealing the differences in skeletal muscle development between Pekin duck and Hanzhong Ma duck.

Whole-Transcriptome RNA Sequencing Uncovers the Global Expression Changes and RNA Regulatory Networks in Duck Embryonic Myogenesis

Duck meat is pivotal in providing high-quality protein for human nutrition, underscoring the importance of studying duck myogenesis. The regulatory mechanisms governing duck myogenesis involve both coding and non-coding RNAs, yet their specific expression patterns and molecular mechanisms remain elusive. To address this knowledge gap, we performed expression profiling analyses of mRNAs, lncRNAs, circRNAs, and miRNAs involved in duck myogenesis using whole-transcriptome RNA-seq. Our analysis identified 1733 differentially expressed (DE)-mRNAs, 1116 DE-lncRNAs, 54 DE-circRNAs, and 174 DE-miRNAs when comparing myoblasts and myotubes. A GO analysis highlighted the enrichment of DE molecules in the extracellular region, protein binding, and exocyst. A KEGG analysis pinpointed pathways related to ferroptosis, PPAR signaling, nitrogen metabolism, cell cycle, cardiac muscle contraction, glycerolipid metabolism, and actin cytoskeleton. A total of 51 trans-acting lncRNAs, including ENSAPLT00020002101 and ENSAPLT00020012069, were predicted to participate in regulating myoblast proliferation and differentiation. Based on the ceRNAs, we constructed lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA ceRNA networks involving five miRNAs (miR-129-5p, miR-133a-5p, miR-22-3p, miR-27b-3p, and let-7b-5p) that are relevant to myogenesis. Furthermore, the GO and KEGG analyses of the DE-mRNAs within the ceRNA network underscored the significant enrichment of the glycerolipid metabolism pathway. We identified five different DE-mRNAs, specifically ENSAPLG00020001677, ENSAPLG00020002183, ENSAPLG00020005019, ENSAPLG00020010497, and ENSAPLG00020017682, as potential target genes that are crucial for myogenesis in the context of glycerolipid metabolism. These five mRNAs are integral to ceRNA networks, with miR-107_R-2 and miR-1260 emerging as key regulators. In summary, this study provides a valuable resource elucidating the intricate interplay of mRNA-lncRNA-circRNA-miRNA in duck myogenesis, shedding light on the molecular mechanisms that govern this critical biological process.

The Integration of Genome-Wide DNA Methylation and Transcriptomics Identifies the Potential Genes That Regulate the Development of Skeletal Muscles in Ducks

DNA methylation is a pivotal epigenetic regulatory mechanism in the development of skeletal muscles. Nonetheless, the regulators responsible for DNA methylation in the development of embryonic duck skeletal muscles remain unknown. In the present study, whole genome bisulfite sequencing (WGBS) and transcriptome sequencing were conducted on the skeletal muscles of embryonic day 21 (E21) and day 28 (E28) ducks. The DNA methylation pattern was found to fall mainly within the cytosine-guanine (CG) context, with high methylation levels in the intron, exon, and promoter regions. Overall, 7902 differentially methylated regions (DMRs) were identified, which corresponded to 3174 differentially methylated genes (DMGs). By using integrative analysis of both WGBS with transcriptomics, we identified 1072 genes that are DMGs that are negatively associated with differentially expressed genes (DEGs). The gene ontology (GO) analysis revealed significant enrichment in phosphorylation, kinase activity, phosphotransferase activity, alcohol-based receptors, and binding to cytoskeletal proteins. The Kyoto Encyclopedia of Genes and Genomes (KEGGs) analysis showed significant enrichment in MAPK signaling, Wnt signaling, apelin signaling, insulin signaling, and FoxO signaling. The screening of enriched genes showed that hyper-methylation inhibited the expression of Idh3a, Got1, Bcl2, Mylk2, Klf2, Erbin, and Klhl38, and hypo-methylation stimulated the expression of Col22a1, Dnmt3b, Fn1, E2f1, Rprm, and Wfikkn1. Further predictions showed that the CpG islands in the promoters of Klhl38, Klf2, Erbin, Mylk2, and Got1 may play a crucial role in regulating the development of skeletal muscles. This study provides new insights into the epigenetic regulation of the development of duck skeletal muscles.

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.

RNA-Seq-based discovery of genetic variants and allele-specific expression of two layer lines and broiler chicken

Recent advances in the selective breeding of broilers and layers have made poultry production one of the fastest-growing industries. In this study, a transcriptome variant calling approach from RNA-seq data was used to determine population diversity between broilers and layers. In total, 200 individuals were analyzed from three different chicken populations (Lohmann Brown (LB), n = 90), Lohmann Selected Leghorn (LSL, n = 89), and Broiler (BR, n = 21). The raw RNA-sequencing reads were pre-processed, quality control checked, mapped to the reference genome, and made compatible with Genome Analysis ToolKit for variant detection. Subsequently, pairwise fixation index (F ST) analysis was performed between broilers and layers. Numerous candidate genes were identified, that were associated with growth, development, metabolism, immunity, and other economically significant traits. Finally, allele-specific expression (ASE) analysis was performed in the gut mucosa of LB and LSL strains at 10, 16, 24, 30, and 60 weeks of age. At different ages, the two-layer strains showed significantly different allele-specific expressions in the gut mucosa, and changes in allelic imbalance were observed across the entire lifespan. Most ASE genes are involved in energy metabolism, including sirtuin signaling pathways, oxidative phosphorylation, and mitochondrial dysfunction. A high number of ASE genes were found during the peak of laying, which were particularly enriched in cholesterol biosynthesis. These findings indicate that genetic architecture as well as biological processes driving particular demands relate to metabolic and nutritional requirements during the laying period shape allelic heterogeneity. These processes are considerably affected by breeding and management, whereby elucidating allele-specific gene regulation is an essential step towards deciphering the genotype to phenotype map or functional diversity between the chicken populations. Additionally, we observed that several genes showing significant allelic imbalance also colocalized with the top 1% of genes identified by the FST approach, suggesting a fixation of genes in cis-regulatory elements.

High-density single nucleotide polymorphism markers reveal the population structure of 2 local chicken genetic resources

Italy counts a large number of local chicken populations, some without a recognized genetic structure, such as Val Platani (VPL) and Cornuta (COS), which represent noteworthy local genetic resources. In this study, the genotype data of 34 COS and 42 VPL, obtained with the Affymetrix Axiom600KChicken Genotyping Array, were used with the aim to investigate the genetic diversity, the runs of homozygosity (ROH) pattern, as well as the population structure and relationship within the framework of other local Italian and commercial chickens. The genetic diversity indices, estimated using different approaches, displayed moderate levels of genetic diversity in both populations. The identified ROH hotspots harbored genes related to immune response and adaptation to local hot temperatures. The results on genetic relationship and population structure reported a clear clustering of the populations according to their geographic origin. The COS formed a nonoverlapping genomic cluster and clearly separated from the other populations, but showed evident proximity to the Siciliana breed (SIC). The VPL highlighted intermediate relationships between the COS-SIC group and the rest of the sample, but closer to the other Italian local chickens. Moreover, VPL showed a complex genomic structure, highlighting the presence of 2 subpopulations that match with the different source of the samples. The results obtained from the survey on genetic differentiation underline the hypothesis that Cornuta is a population with a defined genetic structure. The substructure that characterizes the Val Platani chicken is probably the consequence of the combined effects of genetic drift, small population size, reproductive isolation, and inbreeding. These findings contribute to the understanding of genetic diversity and population structure, and represent a starting point for designing programs to monitor and safeguard these local genetic resources, in order to define a possible official recognition program as breeds.

Effects of Resveratrol on Muscle Inflammation, Energy Utilisation, and Exercise Performance in an Eccentric Contraction Exercise Mouse Model

Eccentric contraction can easily cause muscle damage and an inflammatory response, which reduces the efficiency of muscle contraction. Resveratrol causes anti-inflammatory effects in muscles, accelerates muscle repair, and promotes exercise performance after contusion recovery. However, whether resveratrol provides the same benefits for sports injuries caused by eccentric contraction is unknown. Thus, we explored the effects of resveratrol on inflammation and energy metabolism. In this study, mice were divided into four groups: a control group, an exercise group (EX), an exercise with low-dose resveratrol group (EX + RES25), and an exercise with high-dose resveratrol group (EX + RES150). The results of an exhaustion test showed that the time before exhaustion of the EX + RES150 group was greater than that of the EX group. Tumour necrosis factor-α (Tnfα) mRNA expression was lower in the EX + RES150 group than in the EX group. The energy utilisation of the EX + RES150 group was greater than that of the EX + RES25 group in different muscles. High-dose resveratrol intervention decreased Tnfα mRNA expression and enhanced the mRNA expressions of sirtuin 1, glucose transporter 4, AMP-activated protein kinase α1, and AMP-activated protein kinase α2 in muscles. These results revealed that high-dose resveratrol supplementation can reduce inflammation and oxidation and improve energy utilisation during short-duration high-intensity exercise.

Analysis of mRNA and lncRNA Expression Profiles of Breast Muscle during Pigeon (Columbalivia) Development

The breast muscle is essential for flight and determines the meat yield and quality of the meat type in pigeons. At present, studies about long non-coding RNA (lncRNA) expression profiles in skeletal muscles across the postnatal development of pigeons have not been reported. Here, we used transcriptome sequencing to examine the White-King pigeon breast muscle at four different ages (1 day, 14 days, 28 days, and 2 years old). We identified 12,918 mRNAs and 9158 lncRNAs (5492 known lncRNAs and 3666 novel lncRNAs) in the breast muscle, and 7352 mRNAs and 4494 lncRNAs were differentially expressed in the process of development. We found that highly expressed mRNAs were mainly related to cell-basic and muscle-specific functions. Differential expression and time-series analysis showed that differentially expressed genes were primarily associated with muscle development and functions, blood vessel development, cell cycle, and energy metabolism. To further predict the possible role of lncRNAs, we also conducted the WGCNA and trans/cis analyses. We found that differentially expressed lncRNAs such as lncRNA-LOC102093252, lncRNA-G12653, lncRNA-LOC110357465, lncRNA-G14790, and lncRNA-LOC110360188 might respectively target UBE2B, Pax7, AGTR2, HDAC1, Sox8 and participate in the development of the muscle. Our study provides a valuable resource for studying the lncRNAs and mRNAs of pigeon muscles and for improving the understanding of molecular mechanisms in muscle development.

Comparative transcriptome analysis of Indian domestic duck reveals candidate genes associated with egg production

Egg production is an important economic trait and a key indicator of reproductive performance in ducks. Egg production is regulated by several factors including genes. However the genes involved in egg production in duck remain unclear. In this study, we compared the ovarian transcriptome of high egg laying (HEL) and low egg laying (LEL) ducks using RNA-Seq to identify the genes involved in egg production. The HEL ducks laid on average 433 eggs while the LEL ducks laid 221 eggs over 93 weeks. A total of 489 genes were found to be significantly differentially expressed out of which 310 and 179 genes were up and downregulated, respectively, in the HEL group. Thirty-eight differentially expressed genes (DEGs), including LHX9, GRIA1, DBH, SYCP2L, HSD17B2, PAR6, CAPRIN2, STC2, and RAB27B were found to be potentially related to egg production and folliculogenesis. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis suggested that DEGs were enriched for functions related to glutamate receptor activity, serine-type endopeptidase activity, immune function, progesterone mediated oocyte maturation and MAPK signaling. Protein-protein interaction network analysis (PPI) showed strong interaction between 32 DEGs in two distinct clusters. Together, these findings suggest a mix of genetic and immunological factors affect egg production, and highlights candidate genes and pathways, that provides an understanding of the molecular mechanisms regulating egg production in ducks and in birds more broadly.

Identification of differentially expressed long non-coding RNAs and messenger RNAs involved with muscle development in Dazu black goats through RNA sequencing

This study aimed to explore the genetic basis of muscle development in goats. The transcriptome dataset for differentially expressed lncRNAs (DELs) and differentially expressed genes (DEGs) of goat muscle at different developmental stages were obtained using RNA-Seq. A total of 447,806,481 and 587,559,465 clean reads in the longissimus dorsi muscle of Dazu black goats between 75d embryonic stage and 1d after birth were generated through Illumina paired-end sequencing, and their mapping rates were 89.82 and 90.99%, respectively. Moreover, 4517 DEGs and 648 DELs were identified, and 4784 lncRNA-mRNA targeting relationships were predicted. Gene function annotation results showed that 4101 DEGs were significantly enriched to 1098 GO terms, and 2014 DEGs were significantly enriched to 40 KEGG pathways, including many GO terms and pathways related to muscle development, such as cell differentiation and Wnt signaling pathway. Then, 10 DELs and 20 DEGs were randomly selected for RT-qPCR verification, and the agreement rate between the verification and RNA-Seq results was 90%, indicating the high reliability of the RNA-Seq data analysis. In conclusion, this study obtained several mRNAs and lncRNAs related to the muscle development of Dazu black goats and identified several targeted regulatory pairs of lncRNA-mRNA. This study may serve as a reference to understand the genetic basis and molecular mechanism of muscle development in goats.

Transcriptomic analysis of thigh muscle of Lueyang black-bone chicken in free-range and caged feeding

Lueyang black-bone chicken is free-range in hilly areas and has unique genetic characteristics and excellent muscle quality. However, the molecular mechanisms of breeding mode influence growth and meat quality in Lueyang black-bone chicken are still unclear. Here we analyzed the meat quality and transcriptome data of thigh muscle by comparing free-range and caged modes at the age of 60 and 120 days in Lueyang black-boned chicken. The results demonstrated that the free-range mode could improve the pH value, tenderness, and reducing the hardness of the thigh muscle. Intramuscular fat (IMF) content of the thigh muscle was markedly higher in the caged chickens compared with free-range animals at the age of 60 days. Functional pathway analysis illustrated that tight junction signaling was associated with the formation of slow-twitch fibers in free-range chickens at age of 120 days. All research data proved that the free-range mode could improve muscle quality by promoting the formation of slow-twitch fibers and IMF in thigh muscle in Lueyang black-bone chicken. Based on the animal benefit and healthy, the free-range feeding should be considered during the breeding process of broiler chicken. The results provide good knowledge of the functional molecular mechanisms associated with muscle quality in Lueyang black-bone chicken.

Genome-Wide Association Study Based on Random Regression Model Reveals Candidate Genes Associated with Longitudinal Data in Chinese Simmental Beef Cattle

Simple Summary

Genome-wide association study (GWAS) has become the main approach for detecting functional genes that affects complex traits. For growth traits, the conventional GWAS method can only deal with the single-record traits observed at specific time points, rather than the longitudinal traits measured at multiple time points. Previous studies have reported the random regression model (RRM) for longitudinal data could overcome the limitation of the traditional GWAS model. Here, we present an association analysis based on RRM (GWAS-RRM) for 808 Chinese Simmental beef cattle at four stages of age. Ultimately, 37 significant single-nucleotide polymorphisms (SNPs) and several important candidate genes were screened to be associated with the body weight. Enrichment analysis showed these genes were significantly enriched in the signaling transduction pathway and lipid metabolism. This study not only offers a further understanding of the genetic basis for growth traits in beef cattle, but also provides a robust analytics tool for longitudinal traits in various species.

Abstract

Body weight (BW) is an important longitudinal trait that directly described the growth gain of bovine in production. However, previous genome-wide association study (GWAS) mainly focused on the single-record traits, with less attention paid to longitudinal traits. Compared with traditional GWAS models, the association studies based on the random regression model (GWAS-RRM) have better performance in the control of the false positive rate through considering time-stage effects. In this study, the BW trait data were collected from 808 Chinese Simmental beef cattle aged 0, 6, 12, and 18 months, then we performed a GWAS-RRM to fit the time-varied SNP effect. The results showed a total of 37 significant SNPs were associated with BW. Gene functional annotation and enrichment analysis indicated FGF4, ANGPT4, PLA2G4A, and ITGA5 were promising candidate genes for BW. Moreover, these genes were significantly enriched in the signaling transduction pathway and lipid metabolism. These findings will provide prior molecular information for bovine gene-based selection, as well as facilitate the extensive application of GWAS-RRM in domestic animals.

Transcriptome analysis of embryonic muscle development in Chengkou Mountain Chicken

Background

Muscle is the predominant portion of any meat product, and growth performance and product quality are the core of modern breeding. The embryonic period is highly critical for muscle development, the number, shape and structure of muscle fibers are determined at the embryonic stage. Herein, we performed transcriptome analysis to reveal the law of muscle development in the embryonic stage of Chengkou Mountain Chicken at embryonic days (E) 12, 16, 19, 21.

Results

Diameter and area of muscle fibers exhibited significant difference at different embryonic times(P < 0.01). A total of 16,330 mRNAs transcripts were detected, including 109 novel mRNAs transcripts. By comparing different embryonic muscle development time points, 2,262 in E12vsE16, 5,058 in E12vsE19, 6139 in E12vsE21, 1,282 in E16vsE19, 2,920 in E16vsE21, and 646 in E19vsE21differentially expressed mRNAs were identified. It is worth noting that 7,572 mRNAs were differentially expressed. The time-series expression profile of differentially expressed genes (DEGs) showed that the rising and falling expression trends were significantly enriched. The significant enrichment trends included 3,150 DEGs. GO enrichment analysis provided three significantly enriched categories of significantly enriched differential genes, including 65 cellular components, 88 molecular functions, and 453 biological processes. Through KEGG analysis, we explored the biological metabolic pathways involved in differentially expressed genes. A total of 177 KEGG pathways were enriched, including 19 significant pathways, such as extracellular matrix-receptor interactions. Similarly, numerous pathways related to muscle development were found, including the Wnt signaling pathway (P < 0.05), MAPK signalingpathway, TGF-beta signaling pathway, PI3K-Akt signaling pathway and mTOR signaling pathway. Among the differentially expressed genes, we selected those involved in developing 4-time points; notably, up-regulated genes included MYH1F, SLC25A12, and HADHB, whereas the down-regulated genes included STMN1, VASH2, and TUBAL3.

Conclusions

Our study explored the embryonic muscle development of the Chengkou Mountain Chicken. A large number of DEGs related to muscle development have been identified ,and validation of key genes for embryonic development and preliminary explanation of their role in muscle development. Overall, this study broadened our current understanding of the phenotypic mechanism for myofiber formation and provides valuable information for improving chicken quality.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07740-w.