Skeletal muscle is a complex heterogeneous tissue and characterizing its cellular heterogeneity and transcriptional and epigenetic signatures are important for understanding the details of its ontogeny. In our study, we applied scRNA‐seq and scATAC‐seq to investigate the cell types, molecular features, transcriptional and epigenetic regulation, and patterns of developing bovine skeletal muscle from gestational, lactational and adult stages. Detailed molecular analyses were used to dissect cellular heterogeneity, and we deduced the differentiation trajectory of myogenic cells and uncovered their dynamic gene expression profiles. SCENIC analysis was performed to demonstrate key regulons during cell fate decisions. We explored the future expression states of these heterogeneous cells by RNA velocity analysis and found extensive networks of intercellular communication using the toolkit CellChat. Moreover, the transcriptomic and chromatin accessibility modalities were confirmed to be highly concordant, and integrative analysis of chromatin accessibility and gene expression revealed key transcriptional regulators acting during myogenesis. In bovine skeletal muscle, by scRNA‐seq and scATAC‐seq analysis, different cell types such as adipocytes, endothelial cells, fibroblasts, lymphocytes, monocytes, pericyte cells and eight skeletal myogenic subpopulations were identified at the three developmental stages. The pseudotime trajectory exhibited a distinct sequential ordering for these myogenic subpopulations and eight distinct gene clusters were observed according to their expression pattern. Moreover, specifically expressed TFs (such as MSC, MYF5, MYOD1, FOXP3, ESRRA, BACH1, SIX2 and ATF4) associated with muscle development were predicted, and likely future transcriptional states of individual cells and the developmental dynamics of differentiation among neighbouring cells were predicted. CellChat analysis on the scRNA‐seq data set then classified many ligand–receptor pairs among these cell clusters, which were further categorized into significant signalling pathways, including BMP, IGF, WNT, MSTN, ANGPTL, TGFB, TNF, VEGF and FGF. Finally, scRNA‐seq and scATAC‐seq results were successfully integrated to reveal a series of specifically expressed TFs that are likely to be candidates for the promotion of cell fate transition during bovine skeletal muscle development. Overall, our results outline a single‐cell dynamic chromatin/transcriptional landscape for normal bovine skeletal muscle development; these provide an important resource for understanding the structure and function of mammalian skeletal muscle, which will promote research into its biology.

Based on the previous GWAS research related to bovine heat tolerance trait, this study aimed to explore the effect of myosin-1a (MYO1A) gene on bovine heat tolerance trait, and find the molecular markers related to the heat tolerance of Chinese cattle. In our study, four novel candidate SNPs highly conserved in B. indicus breeds but barely existed in B. taurus were identified in MYO1A gene according to Bovine Genome Variation Database and Selective Signatures (BGVD). PCR and DNA sequencing were used to genotype 1072 individuals including 34 Chinese indigenous cattle breeds as well as Angus and Indian zebu. Two synonymous mutations (rs208210464 and rs110123931), one missense mutation (rs209999142; Phe172Ser), and one intron mutation (rs135771836) were detected. The frequencies of mutant alleles of the four SNPs gradually increased from northern groups to southern groups of Chinese cattle, which was consistent with the distribution of various climatic conditions of China. Additionally, four SNPs were significantly associated with four climatic conditions including annual mean temperature (T), relative humidity (H), temperature-humidity index (THI), and average annual sunshine hours (100-cloudiness) (SR). Among these, rs209999142 and Hap 1/1 had better performance than others. Our results suggested that rs209999142 was associated with heat-tolerance trait and rs208210464, rs110123931, and rs135771836 showed high phenotypic effect on heat-tolerance trait because of the strong linkage with rs209999142. These SNPs could be used as candidates for marker-assisted selection (MAS) in cattle breeding.

Copy number variation is a part of genomic structural variation and has caused widespread concern. According to the results of high-throughput screening of the MLLT10 gene, we found that the copy number variation region of the MLLT10 gene was correlated with bovine growth traits. We aimed to detect the MLLT10 gene copy number variation and provide materials for the Chinese yellow cattle breed. In this study, the SPSS software was used to analyze the correlation among the copy number type of six different cattle breeds (i.e., Qinchuan, Xianan, Jiaxian, Yanbian, Sinan, Yunling) and the corresponding growth traits. The results showed the following: In Qinchuan cattle, the copy number duplication type was greater than the deletion and normal types; in Xianan cattle, the copy number duplication and normal types were less as compared with the deletion type; and in Yunling cattle, the frequency of the duplication type was dominant among the three types of copy number variants. The correlation analysis result showed that there is a significant correlation between the copy number variation (CNV) of the MLLT10 gene and the growth traits of three cattle breeds. Furthermore, correlation analysis showed that MLLT10 CNV had positive effects on growth traits such as hip width, rump length, hucklebone width, and cannon bone circumference (p < 0.05). This study provides a basis for the molecular-assisted marker breeding of cattle and contributes to the breeding of cattle.

Copy number variant (CNV), a common genetic polymorphism, is closely related to the phenotypic variation traits of organisms. Vesicle-associated membrane protein 7 gene (VAMP7) codes a protein, which is a member of the SNARE proteins family and plays an important role in the process of intracellular vesicle transport. In this study, a total of four cattle breeds (Yunling cattle, Xianan cattle, Pinan cattle, Jiaxian red cattle) were used to investigate the copy numbers, and we found an association relationship between CNV of VAMP7 gene and growth traits of cattle by SPSS 20.0 software. The results showed that the CNV type of VAMP7 gene in four cattle breeds had the same distribution, Duplication type occupies a dominant position among the four varieties. In Yunling cattle, the Duplication type of VAMP7 is significantly related to the height at the hip cross (p < 0.05), Individuals with Duplication type commonly have less performance on growth and development, which indicates that the Duplication type of the VAMP7 gene may have a negative effect on cattle growth. Individuals with the other two CNV types may become the breeding direction of the VAMP7 gene. This study provided a new perspective and basic material for the molecular genetics of the CNV of the VAMP7 gene, and also promoted the breeding progress of Chinese local cattle.

Serpin peptidase inhibitor, clade A (alpha-1 antiproteinase, antitrypsin), member 3 (SERPINA3) belongs to the serine protease inhibitor family A subtype, and contains 8 genes from SERPINA3-1 to SERPINA3-8. Although the regulatory effects of these 8 genes have been revealed one by one in recent years, the related effects of SERPINA3-1 gene on cattle growth is still unclear. This study used quantitative Real time PCR (qPCR) to detect the type of copy number variation (CNV) of SERPINA3-1 gene in a total of 542 Chinese cattle, and expression of SERPINA3-1 gene in different tissues of Qinchuan cattles (adult) on mRNA level. Then association analysis was conducted between the detection results and cattle growth traits. The results showed that the Duplication type in SERPINA3-1 gene performed better on the growth traits and the CNV was significantly correlated with multiple growth traits (p < 0.05). In addition, SERPINA3-1 gene has different expression conditions in different tissues, results showed that SERPINA3-1 gene has a low expression in muscle. In conclusion, we speculate that the SERPINA3-1 gene can be used as a molecular marker and the result of this study could be a basic material for candidate functional genes for beef cattle growth and development. In order to detect the gene expression diversification of the SERPINA3-1 gene, blood samples were collected from five Chinese cattle breeds, we detected related signal and made an associated analyze with cattle growth traits. We determined the copy number variation distribution of the SERPINA3-1 gene in cattle populations and found that the SERPINA3-1 gene has a certain promoting effect on the growth and development of Chinese cattle. For example, Pinan cattle with Duplication type copy number have a better performance on growth traits. This study has enriched the candidate genes of Chinese cattle molecular breeding and provided basic data for Chinese cattle breeding.

Mutations are heritable changes at the base level of genomic DNA. Furthermore, mutations lead to genetic polymorphisms and may alter animal growth phenotypes. Our previous study found that mutations in the bovine Actin-like protein 8 (ACTL8) gene may be involved in muscle growth and development. This study explored several mutations of the ACTL8 gene and their influence on body size in Chinese beef cattle, as well as tested the tissue expression profile of the ACTL8 gene in Qinchuan cattle at different ages. Five single nucleotide polymorphisms (SNPs) (including one synonymous mutation (c.2135552895G > A)) and two insertion/deletion polymorphisms (indels) were identified in the ACTL8 gene from 1138 cattle by DNA-seq, RFLP and other methods. Then, the expression profile of the ACTL8 gene in Qinchuan cattle showed that it was expressed in heart, spleen, lung, liver, muscle, and fat tissues. Moreover, the expression level of ACTL8 was increased with cattle growth (p < 0.01). The ACTL8 mRNA expression level in kidney and muscle tissues was the highest in the calves, while lowest in the fetal stage. Overall, we showed that the mutations could act as markers in beef molecular breeding and selection of the growth traits of cattle.