The yak is the largest meat-producing mammal around the Tibetan Plateau, and it plays an important role in the economic development and maintenance of the ecological environment throughout much of the Asian highlands. Understanding the genetic components of body weight is key for future improvement in yak breeding; therefore, genome-wide association studies (GWAS) were performed, and the results were used to mine plant and animal genetic resources. We conducted whole genome sequencing on 406 Maiwa yaks at 10 × coverage. Using a multiple loci mixed linear model (MLMM), fixed and random model circulating probability unification (FarmCPU), and Bayesian-information and linkage-disequilibrium iteratively nested keyway (BLINK), we found that a total of 25,000 single-nucleotide polymorphisms (SNPs) were distributed across chromosomes, and seven markers were identified as significantly (p-values < 3.91 × 10⁻⁷) associated with the body weight trait,. Several candidate genes, including MFSD4, LRRC37B, and NCAM2, were identified. This research will help us achieve a better understanding of the genotype–phenotype relationship for body weight.

Sirtuin 1 (SIRT1) overexpression significantly inhibits lipid deposition during yak intramuscular preadipocyte (YIMA) differentiation; however, the regulatory mechanism remains unknown. We elucidated the role of SIRT1 in YIMA differentiation using lentivirus-mediated downregulation technology and conducted mRNA-seq and ChIP-seq assays using H3K9ac antibodies after SIRT1 overexpression in order to reveal SIRT1 targets during YIMA adipogenesis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed in order to identify the functional annotation of common genes. In addition, a potential target of SIRT1 was selected to verify its effects on the differentiation and proliferation of YIMAs. SIRT1 interfered with lipid deposition and promoted YIMA differentiation. In total, 143,518 specific peaks were identified after SIRT1 overexpression, where genes associated with downregulation peaks were enriched in transcription, gene expression, lipid-related processes, and classical lipid-related pathways. The H3K9ac signal in the whole genome promoter region (2 kb upstream and downstream of the transcription start site (TSS)) was weakened, and the peaks were distributed across all gene functional regions. Genes that lost signals in their TSS region or gene body region were enriched in both biological processes and pathways associated with lipogenesis. The ChIP-seq results revealed 714 common differential genes in mRNA-seq, which were enriched in “MAPK signaling”, “lipid and atherosclerosis”, “mTOR signaling”, and “FoxO signaling” pathways. A total of 445 genes were downregulated in both their H3K9ac signals and mRNA expression, and one of their most significantly enriched pathways was FoxO signaling. Nine genes (FBP2, FPGT, HSD17B11, KCNJ15, MAP3K20, SLC5A3, TRIM23, ZCCHC10, and ZMYM1) lost the H3K9ac signal in their TSS regions and had low mRNA expression, and three genes (KCNJ15, TGM3, and TRIM54) had low expression but lost their H3K9ac signal in the gene body region. The interference of TRIM23 significantly inhibited fat deposition during preadipocyte differentiation and promoted cell proliferation by increasing S-phase cell numbers. The present study provides new insights into the molecular mechanism of intramuscular fat content deposition and the epigenetic role of SIRT1 in adipocyte differentiation.

Due to its prominent secretory activity, adipose tissue (AT) is now considered a major player in the crosstalk between organs, especially with skeletal muscle. In which, exosomes are effective carriers for the intercellular material transfer of a wide range of molecules that can influence a series of physiological and pathological processes in recipient cells. Considering their underlying roles, the regulatory mechanisms of adipose-secreted exosomes and their cellular crosstalk with skeletal muscle have received great attention in the field. In this review, we describe what is currently known of adipose-secreted exosomes, as well as their applications in skeletal muscle pathophysiology.

The concentration of intramuscular fat (IMF) is a crucial determinant of yak meat quality. However, the molecular mechanisms that regulate IMF in yak remain largely elusive. In our study, we conducted transcriptome sequencing on the longissimus dorsi muscle tissues of yaks with varying IMF contents. We then filtered differentially expressed genes (DEGs), microRNAs (DEMs), and long non-coding RNAs (DELs) to elucidate potential regulatory pathways of adipogenesis in yaks. Overall, our research sheds light on an array of potential mRNAs and noncoding RNAs implicated in IMF deposition and elaborates on the role of HIF1α in yaks. These findings contribute valuable insights that can serve as a guide for further research into the molecular mechanisms governing IMF deposition.

Yaks have evolved several breeds or genetic resources owing to their geographical and ecological environment, and investigating the genetic construction of body size among breeds is key for breeding. Here, a genome-wide association study (GWAS) was performed for five body size traits in 31 yak breeds and genetic resources. The information from clustering individuals according to their habitats was used for kinship grouping in the compressed mixed linear model (CMLM). We named this approach the pCMLM method. A total of 3,584,464 high-quality single nucleotide polymorphisms (SNPs) were obtained, and six markers were found to be significantly associated with height by pCMLM. Four candidate genes, including FXYD6, SOHLH2, ADGRB2, and OSBPL6, were identified. Our results show that when CMLM cannot identify optimal clustering groups, pCMLM can provide sufficient associated results based on population information. Moreover, this study provides basic information on the gene localization of quantitative traits of body size among yak breeds.

Meiotic recombination is initiated by the formation of programmed DNA double-strand breaks during spermatogenesis. PRDM9 determines the localization of recombination hotspots by interacting with several protein complexes in mammals. The function of PRDM9 is not well understood during spermatogenesis in mice or yaks. In this study, we applied yeast two-hybrid assays combined with next-generation sequencing techniques to screen the complete set of PRDM9-interacting proteins and explore its novel functions in yak spermatogenesis. Our results showed that 267 PRDM9-interacting proteins were identified. The gene ontology (GO) analysis of the interacting proteins revealed that the GO terms were primarily associated with spermatogenesis, positive regulation of double-strand break repair via homologous recombination, RNA splicing, the ubiquitin-dependent ERAD pathway, and other biological processes. MKX and PDCD5 were verified to be strongly interacting with PRDM9 and expressed in prophase I of meiosis in both mouse and yak testes. The localizations of RNA splicing genes including THOC5, DDX5, and XRCC6 were expressed in spermatocytes. Cattleyak is the hybrid offspring of a yak and a domestic cow, and the male offspring are sterile. The gene expression of the interacting proteins was also examined in the sterile male hybrid of yak and cattle. Among the 58 detected genes, 55 were downregulated in cattleyak. In conclusion, we established a complete PRDM9 interaction network, and a novel function of PRDM9 was identified, which will further promote our understanding of spermatogenesis. It also provides new insights for the study of hybrid male sterility.