Fatty acid composition of muscle is an important trait contributing to meat quality. Recently, genome-wide association study (GWAS) has been extensively used to explore the molecular mechanism underlying important traits in cattle. In this study, we performed GWAS using high density SNP array to analyze the association between SNPs and fatty acids and evaluated the accuracy of genomic prediction for fatty acids in Chinese Simmental cattle. Using the BayesB method, we identified 35 and 7 regions in Chinese Simmental cattle that displayed significant associations with individual fatty acids and fatty acid groups, respectively. We further obtained several candidate genes which may be involved in fatty acid biosynthesis including elongation of very long chain fatty acids protein 5 (ELOVL5), fatty acid synthase (FASN), caspase 2 (CASP2) and thyroglobulin (TG). Specifically, we obtained strong evidence of association signals for one SNP located at 51.3 Mb for FASN using Genome-wide Rapid Association Mixed Model and Regression-Genomic Control (GRAMMAR-GC) approaches. Also, region-based association test identified multiple SNPs within FASN and ELOVL5 for C14:0. In addition, our result revealed that the effectiveness of genomic prediction for fatty acid composition using BayesB was slightly superior over GBLUP in Chinese Simmental cattle. We identified several significantly associated regions and loci which can be considered as potential candidate markers for genomics-assisted breeding programs. Using multiple methods, our results revealed that FASN and ELOVL5 are associated with fatty acids with strong evidence. Our finding also suggested that it is feasible to perform genomic selection for fatty acids in Chinese Simmental cattle.

China exhibits a great diversity of ecosystems and abundant cattle resources, with nearly 30 million cattle from 53 indigenous breeds reared in specific geographic regions. To explore the genetic diversity and population structure of Chinese indigenous cattle, a population genetic analysis at both the individual and population levels was conducted and the admixture analysis was performed. We genotyped 572 samples from 20 Chinese indigenous cattle breeds using GeneSeek Genomic Profiler Bovine LD (GGP-LD, 30 K) and downloaded the published data of 77 samples from 4 worldwide commercial breeds genotyped with Illumina BovineSNP50 Beadchip (SNP50, 50 K). In principal component analysis (PCA) and neighbour-joining (NJ) tree analysis, samples of the same breeds were grouped together, leading to clear separation from other breeds. And Chinese indigenous cattle were clustered into two groups of southern and northern breeds, originated from Asian Bos indicus lineage and Eurasian Bos taurus lineage, respectively. In STRUCTURE K = 2, a clear transition occurred from the northern breeds to the southern breeds. Additionally, the northern breeds contained a smaller Eurasian taurine (62.5%) descent proportion than that reported previously (more than 90%). In STRUCTURE K = 3, a distinct descent was detected in the southern Chinese breeds, which could reflect a long-term selection history of Chinese indigenous cattle. The results from TreeMix and f3 statistic provided the evidence of an admixture history between southern breeds and northern breeds. Consistent with the observed geographical distributions, Chinese indigenous cattle were divided into two genetic clusters, northern indigenous cattle and southern indigenous cattle. Three improved breeds in the northern area also exhibited northern indigenous ancestry. We found that the breeds distributed in the northern China showed more southern lineage introgression than previously reported. Central-located populations appeared to the admixture between southern and northern lineages, and introgression events from European cattle were observed in Luxi Cattle, Qinchuan Cattle and Jinnan Cattle. The study revealed the population structures and levels of admixture pattern among Chinese indigenous cattle, shedding light on the origin and evolutionary history of these breeds.

Intramuscular fat deposition or marbling is essential for high quality beef. The molecular mechanism of adipogenesis in skeletal muscle remains largely unknown. In this study, we isolated Platelet-derived growth factor receptor α (PDGFRα) positive progenitor cells from fetal bovine skeletal muscle and induced into adipocytes. Using miRNAome sequencing, we revealed that bta-miR-23a was an adipogenic miRNA mediating bovine adipogenesis in skeletal muscle. The expression of bta-miR-23a was down-regulated during differentiation of PDGFRα⁺ progenitor cells. Forced expression of bta-miR-23a mimics reduced lipid accumulation and inhibited the key adipogenic transcription factor peroxisome proliferative activated receptor gamma (PPARγ) and CCAAT/enhancer binding protein alpha (C/EBPα). Whereas down-regulation of bta-miR-23a by its inhibitors increased lipid accumulation and expression of C/EBPα, PPARγ and fatty acid-binding protein 4 (FABP4). Target prediction analysis revealed that ZNF423 was a potential target of bta-miR-23a. Dual-luciferase reporter assay revealed that bta-miR-23a directly targeted the 3′-UTR of ZNF423. Together, our data showed that bta-miR-23a orchestrates early intramuscular adipogeneic commitment as an anti-adipogenic regulator which acts by targeting ZNF423.

Understanding the population structure and genetic bases of well-adapted cattle breeds to local environments is one of the most essential tasks to develop appropriate genetic improvement programs. We performed a comprehensive study to investigate the population structure, divergence and selection signatures at genome-wide level in diverse Chinese local cattle using Bovine HD SNPs array, including two breeds from North China, one breed from Northwest China, three breeds from Southwest China and two breeds from South China. Population genetic analyses revealed the genetic structures of these populations were mostly related to the geographic locations. Notably, we detected 294 and 1263 candidate regions under selection using the di and iHS approaches, respectively. A series of group-specific and breed-specific candidate genes were identified, which are involved in immune response, sexual maturation, stature related, birth and bone weight, embryonic development, coat colors and adaptation. Furthermore, haplotype diversity and network pattern for candidate genes, including LPGAT1, LCORL, PPP1R8, RXFP2 and FANCA, suggest that these genes have been under differential selection pressures in various environmental conditions. Our results shed insights into diverse selection during breed formation in Chinese local cattle. These findings may promote the application of genome-assisted breeding for well-adapted local breeds with economic and ecological importance.

The objective of the present study was to perform a genome-wide association study (GWAS) for growth curve parameters using nonlinear models that fit original weight–age records. In this study, data from 808 Chinese Simmental beef cattle that were weighed at 0, 6, 12, and 18 months of age were used to fit the growth curve. The Gompertz model showed the highest coefficient of determination (R² = 0.954). The parameters’ mature body weight (A), time-scale parameter (b), and maturity rate (K) were treated as phenotypes for single-trait GWAS and multi-trait GWAS. In total, 9, 49, and 7 significant SNPs associated with A, b, and K were identified by single-trait GWAS; 22 significant single nucleotide polymorphisms (SNPs) were identified by multi-trait GWAS. Among them, we observed several candidate genes, including PLIN3, KCNS3, TMCO1, PRKAG3, ANGPTL2, IGF-1, SHISA9, and STK3, which were previously reported to associate with growth and development. Further research for these candidate genes may be useful for exploring the full genetic architecture underlying growth and development traits in livestock.

Improving the genetic process of growth traits is one of the major goals in the beef cattle industry, as it can increase meat production and reduce the cost of raising animals. Although several quantitative trait loci affecting growth traits in beef cattle have been identified, the genetic architecture of these economically important traits remains elusive. This study aims to map single nucleotide polymorphisms (SNPs) and genes associated with birth weight (BW), yearling weight (YW), average daily gain from birth to yearling (BYADG), and body weight at the age of 18 months (18MW) in a Chinese Simmental beef cattle population using a weighted, single-step, genome-wide association study (wssGWAS). Phenotypic and pedigree data from 6022 animals and genotypes from 744 animals (596,297 SNPs) were used for an association analysis. The results showed that 66 genomic windows explained 1.01–20.15% of the genetic variance for the four examined traits, together with the genes near the top SNP within each window. Furthermore, the identified genomic windows (>1%) explained 50.56%, 57.71%, 61.78%, and 37.82% of the genetic variances for BW, YW, BYADG, and 18MW, respectively. Genes with potential functions in muscle development and regulation of cell growth were highlighted as candidates for growth traits in Simmental cattle (SQOR and TBCB for BW, MYH10 for YW, RLF for BYADG, and ARHGAP31 for 18MW). Moreover, we found 40 SNPs that had not previously been identified as being associated with growth traits in cattle. These findings will further advance our understanding of the genetic basis for growth traits and will be useful for the molecular breeding of BW, YW, BYADG, and 18MW in the context of genomic selection in beef cattle.