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.

Genome-wide single nucleotide polymorphism (SNP) arrays can be used to explore homozygosity segments, where two haplotypes inherited from the parents are identical. In this study, we identified a total of 27,358 runs of homozygosity (ROH) with an average of 153 ROH events per animal in Chinese local cattle. The sizes of ROH events varied considerably ranging from 0.5 to 66 Mb, with an average length of 1.22 Mb. The highest average proportion of the genome covered by ROH (~11.54% of the cattle genome) was found in Nanda cattle (NDC) from South China, whereas the lowest average proportion (~3.1%) was observed in Yanhuang cattle (YHC). The average estimated FROH ranged from 0.03 in YHC to 0.12 in NDC. For each of three ROH classes with different sizes (Small 0.5–1 Mb, Medium 1–5 Mb and Large >5 Mb), the numbers and total lengths of ROH per individual showed considerable differences across breeds. Moreover, we obtained 993 to 3603 ROH hotspots (which were defined where ROH frequency at a SNP within each breed exceeded the 1% threshold) among eight cattle breeds. Our results also revealed several candidate genes embedded with ROH hotspots which may be related to environmental conditions and local adaptation. In conclusion, we generated baselines for homozygosity patterns in diverse Chinese cattle breeds. Our results suggested that selection has, at least partially, played a role with other factors in shaping the genomic patterns of ROH in Chinese local cattle and might provide valuable insights for understanding the genetic basis of economic and adaptive traits.

Recent advances in high-throughput genotyping technologies have provided the opportunity to map genes using associations between complex traits and markers. Genome-wide association studies (GWAS) based on either a single marker or haplotype have identified genetic variants and underlying genetic mechanisms of quantitative traits. Prompted by the achievements of studies examining economic traits in cattle and to verify the consistency of these two methods using real data, the current study was conducted to construct the haplotype structure in the bovine genome and to detect relevant genes genuinely affecting a carcass trait and a meat quality trait. Using the Illumina BovineHD BeadChip, 942 young bulls with genotyping data were introduced as a reference population to identify the genes in the beef cattle genome significantly associated with foreshank weight and triglyceride levels. In total, 92,553 haplotype blocks were detected in the genome. The regions of high linkage disequilibrium extended up to approximately 200 kb, and the size of haplotype blocks ranged from 22 bp to 199,266 bp. Additionally, the individual SNP analysis and the haplotype-based analysis detected similar regions and common SNPs for these two representative traits. A total of 12 and 7 SNPs in the bovine genome were significantly associated with foreshank weight and triglyceride levels, respectively. By comparison, 4 and 5 haplotype blocks containing the majority of significant SNPs were strongly associated with foreshank weight and triglyceride levels, respectively. In addition, 36 SNPs with high linkage disequilibrium were detected in the GNAQ gene, a potential hotspot that may play a crucial role for regulating carcass trait components.