During the postnatal stages, skeletal muscle development undergoes a series of meticulously regulated alterations in gene expression. However, limited studies have employed chromatin accessibility to unravel the underlying molecular mechanisms governing muscle development in yak species. Therefore, we conducted an analysis of both gene expression levels and chromatin accessibility to comprehensively characterize the dynamic genome-wide chromatin accessibility during muscle growth and development in the Tianzhu white yak, thereby elucidating the features of accessible chromatin regions throughout this process. Initially, we compared the differences in chromatin accessibility between two groups and observed that calves exhibited higher levels of chromatin accessibility compared to adult cattle, particularly within ±2 kb of the transcription start site (TSS). In order to investigate the correlation between alterations in chromatin accessible regions and variations in gene expression levels, we employed a combination of ATAC-seq and RNA-seq techniques, leading to the identification of 18 central transcriptional factors (TFs) and 110 key genes with significant effects. Through further analysis, we successfully identified several TFs, including Sp1, YY1, MyoG, MEF2A and MEF2C, as well as a number of candidate genes (ANKRD2, ANKRD1, BTG2 and LMOD3) which may be closely associated with muscle growth and development. Moreover, we constructed an interactive network program encompassing hub TFs and key genes related to muscle growth and development. This innovative approach provided valuable insights into the molecular mechanism underlying skeletal muscle development in the postnatal stages of Tianzhu white yaks while also establishing a solid theoretical foundation for future research on yak muscle development.

The production performance of Jeryak, resulting from the F1 generation of the cross between Gannan yak and Jersey cattle, exhibits a significantly superior outcome compared with that of Gannan yak. Therefore, we used an RNA-seq approach to identify differentially expressed mRNAs (DEMs) and differentially expressed lncRNAs (DELs) influencing muscle growth and development in Gannan yaks and Jeryaks. A total of 304 differentially expressed lncRNAs and 1819 differentially expressed mRNAs were identified based on the screening criteria of |log 2 FC| > 1 and FDR < 0.05. Among these, 132 lncRNAs and 1081 mRNAs were found to be down-regulated, while 172 lncRNAs and 738 mRNAs were up-regulated. GO and KEGG analyses showed that the identified DELs and DEMs were enriched in the entries of pathways associated with muscle growth and development. On this basis, we constructed an lncRNA–mRNA interaction network. Interestingly, two candidate DELs (MSTRG.16260.9 and MSTRG.22127.1) had targeting relationships with 16 (MYC, IGFBP5, IGFBP2, MYH4, FGF6, etc.) genes related to muscle growth and development. These results could provide a basis for further studies on the roles of lncRNAs and mRNAs in muscle growth in Gannan yaks and Jeryak breeds.

Jeryak is the F1 generation of the cross between Gannan yak and Jersey cattle, which has the advantages of fast growth and high adaptability. The growth and development of skeletal muscle is closely linked to meat production and the quality of meat. However, the molecular regulatory mechanisms of muscle growth differences between Gannan yak and Jeryak analyzed from the perspective of chromatin opening have not been reported. In this study, ATAC-seq was used to analyze the difference of chromatin openness in longissimus muscle of Gannan yak and Jeryak. It was found that chromatin accessibility was more enriched in Jeryak compared to Gannan yak, especially in the range of the transcription start site (TSS) ± 2 kb. GO and KEGG enrichment analysis indicate that differential peak-associated genes are involved in the negative regulation of muscle adaptation and the Hippo signaling pathway. Integration analysis of ATAC-seq and RNA-seq revealed overlapping genes were significantly enriched during skeletal muscle cell differentiation and muscle organ morphogenesis. At the same time, we screened FOXO1, ZBED6, CRY2 and CFL2 for possible involvement in skeletal muscle development, constructed a genes and transcription factors network map, and found that some transcription factors (TFs), including YY1, KLF4, KLF5 and Bach1, were involved in skeletal muscle development. Overall, we have gained a comprehensive understanding of the key factors that impact skeletal muscle development in various breeds of cattle, providing new insights for future analysis of the molecular regulatory mechanisms involved in muscle growth and development.

Muscle formation directly determines meat production and quality. The non-SMC condensin I complex subunit G (NCAPG) is strongly linked to the growth features of domestic animals because it is essential in controlling muscle growth and development. This study aims to elucidate the tissue expression level of the bovine NCAPG gene, and determine the key transcription factors for regulating the bovine NCAPG gene. In this study, we observed that the bovine NCAPG gene exhibited high expression levels in longissimus dorsi and spleen tissues. Subsequently, we cloned and characterized the promoter region of the bovine NCAPG gene, consisting of a 2039 bp sequence, through constructing the deletion fragment double-luciferase reporter vector and site-directed mutation-identifying core promoter region with its key transcription factor binding site. In addition, the key transcription factors of the core promoter sequence of the bovine NCAPG gene were analyzed and predicted using online software. Furthermore, by integrating overexpression experiments and the electrophoretic mobility shift assay (EMSA), we have shown that cAMP response element binding protein 1 (CREB1) and myogenic differentiation 1 (MYOD1) bind to the core promoter region (−598/+87), activating transcription activity in the bovine NCAPG gene. In conclusion, these findings shed important light on the regulatory network mechanism that underlies the expression of the NCAPG gene throughout the development of the muscles in beef cattle.

Unsaturated fatty acids (UFAs) in beef play a vital role in promoting human health. Long-chain fatty acyl-CoA synthase 1 (ACSL1) is a crucial gene for UFA synthesis in bovine adipocytes. To investigate the protein expression profile during UFA synthesis, we performed a proteomic analysis of bovine adipocytes by RNA interference and non-interference with ACSL1 using label-free techniques. A total of 3558 proteins were identified in both the NC and si-treated groups, of which 1428 were differentially expressed proteins (DEPs; fold change ≥ 1.2 or ≤ 0.83 and p-value < 0.05). The enrichment analysis of the DEPs revealed signaling pathways related to UFA synthesis or metabolism, including cAMP, oxytocin, fatty acid degradation, glycerol metabolism, insulin, and the regulation of lipolysis in adipocytes (p-value < 0.05). Furthermore, based on the enrichment analysis of the DEPs, we screened 50 DEPs that potentially influence the synthesis of UFAs and constructed an interaction network. Moreover, by integrating our previously published transcriptome data, this study established a regulatory network involving differentially expressed long non-coding RNAs (DELs), highlighting 21 DEPs and 13 DELs as key genes involved in UFA synthesis. These findings present potential candidate genes for further investigation into the molecular mechanisms underlying UFA synthesis in bovines, thereby offering insights to enhance the quality of beef and contribute to consumer health in future studies.

The Gannan yak, a superior livestock breed found on the Tibetan Plateau, exhibits significantly enhanced body size, weight, and growth performance in comparison to the Tianzhu white yak. MiRNAs play a pivotal role in regulating muscle growth by negatively modulating target genes. In this study, we found the average diameter, area, and length of myofibers in Gannan yaks were significantly higher than those of Tianzhu white yaks. Further, we focused on analyzing the longissimus dorsi muscle from both Gannan yaks and Tianzhu white yaks through transcriptome sequencing to identify differentially expressed (DE)miRNAs that influence skeletal muscle development. A total of 254 DE miRNAs were identified, of which 126 miRNAs were up-regulated and 128 miRNAs were down-regulated. GO and KEGG enrichment analysis showed that the target genes of these DE miRNAs were significantly enriched in signaling pathways associated with muscle growth and development. By constructing a DE miRNA- DE mRNA interaction network, we screened 18 key miRNAs, and notably, four of the candidates (novel-m0143-3p, novel-m0024-3p, novel-m0128-5p, and novel-m0026-3p) targeted six genes associated with muscle growth and development (DDIT4, ADAMTS1, CRY2, AKIRIN2, SIX1, and FOXO1). These findings may provide theoretical references for further studies on the role of miRNAs in muscle growth and development in Gannan yaks.