This study aimed to investigate the nutritional properties of yak milk in various areas of Gannan. The milk composition analyzer, automatic amino acid analyzer, and flavor analyzer were used to detect the conventional nutrients, amino acids, and volatile flavor substances of 249 yak milks in Meiren grassland, Xiahe grassland, and Maqu grassland (hereinafter referred to as Meiren yak, Xiahe yak, and Maqu yak) in the Gannan area. The results showed that the fat content of Meiren yak milk was significantly higher than that of Maqu yak and Xiahe yak (p < 0.05). The protein content of Meiren yak milk was significantly higher than that of Xiahe yak (p < 0.05), but not significantly different from that of Maqu yak (p > 0.05). The casein content in the milk of Maqu yak was significantly higher than that of Meiren yak and Xiahe yak (p < 0.05). There was no significant difference in the lactose content of yak milk in the three regions (p > 0.05). The content of glutamic acid in the milk of Meiren yak, Xiahe yak, and Maqu yak was noticeably high, which was 1.03 g/100 g, 1.07 g/100 g, and 1.10 g/100 g, respectively. The total amino acid (TAA) content was 4.78 g/100 g, 4.87 g/100 g, and 5.0 g/100 g, respectively. The ratios of essential amino acids (EAA) and total amino acids (TAA) in the milk of Meiren yak, Xiahe yak, and Maqu yak were 42.26%, 41.27%, and 41.39%, respectively, and the ratios of essential amino acids (EAA) and nonessential amino acids (NEAA) were 73.19%, 70.28%, and 70.61%, respectively. In the yak milk samples collected from three different regions, a total of 34 volatile flavor compounds were detected, including 10 aldehydes, five esters, six ketones, four alcohols, two acids, and seven others. The main flavor substances qualitatively obtained from Meiren yak milk were ethyl acetate, n-valeraldehyde, acetic acid, heptanal, and n-hexanal. Xiahe yak milk mainly contains ethyl acetate, isoamyl alcohol, n-valeraldehyde, heptanal, and ethyl butyrate. Maqu yak milk mainly contains ethyl acetate, n-valeraldehyde, isoamyl alcohol, heptanal, ethyl butyrate, and n-hexanal. Principal component analysis showed that the flavor difference between Xiahe yak and Maqu yak was small, while the flavor difference between Xiahe yak, Maqu yak, and Meiren yak was large. The findings of this research can serve as a foundation for the future advancement and application of yak milk.

The IQ motif containing GTPase activating protein 2 (IQGAP2) gene functions as a tumor suppressor, reducing the malignant properties of breast cancer cells. The circulating cartilage acidic protein 1 (CRTAC1) gene, present in the whey protein fraction of dairy cows throughout lactation, is significantly correlated with fatty acids in milk. In this study, we investigated the correlation between single nucleotide polymorphisms (SNPs) in the IQGAP2 and CRTAC1 genes and milk quality traits in Gannan yaks, aiming to identify potential molecular marker loci for enhancing milk quality. Using the Illumina Yak cGPS 7K liquid chip, we genotyped 162 yaks and identified five SNPs in the IQGAP2 (g.232,769C>G, g.232,922G>C) and CRTAC1 (g.4,203T>C, g.5,348T>G, g.122,451T>C) genes. Genetic polymorphism analysis revealed that these five SNPs were moderately polymorphic and in Hardy–Weinberg equilibrium. An association analysis results showed that, at the g.232,769C>G locus of the IQGAP2 gene, the heterozygous CG genotype had significantly higher lactose content than the CC and GG homozygous genotypes (p < 0.05). Similarly, at the g.232,922G>C locus, the heterozygous GC and mutant CC genotypes significantly increased the contents of milk fat, lactose, and total solids (TS) (p < 0.05). In the CRTAC1 gene (g.4,203T>C, g.5,348T>G, g.122,451T>C), the mutant CC genotype significantly increased milk fat content, while the heterozygous TG genotype significantly increased lactose content (p < 0.05). In summary, mutations at the loci of g.232,769C>G, g.232,922G>C, g.4,203T>C, g.5,348T>G, and g.122,451T>C significantly elevated the lactose, milk fat, and TS content in Gannan yak milk, providing potential molecular marker candidates for improving Gannan yak milk quality.

The tyrosine protein kinase (TXK) gene, as a member of the non-receptor tyrosine kinase Tec family, plays a vital role in signal transduction mediation. Phospholipase C epsilon 1 (PLCE1), a membrane-associated enzyme, is of paramount importance for the differentiation of myoblasts and the normal functioning of muscle tissue. In recent years, both of these genes have been reported to be associated with the economic traits of animals. This study aimed to investigate the relationship between single nucleotide polymorphisms (SNPs) in the TXK and PLCE1 genes and growth traits in Ashidan yaks and to search for potential molecular marker loci that can influence Ashidan yak breeding. A cGPS liquid microarray was utilized to genotype 232 Ashidan yaks and to analyze correlations between two SNP loci in the TXK and PLCE1 genes and yak body weight, body height, body length, and chest circumference at different periods. The results indicated that the g.55,999,531C>T locus of the TXK gene and the g.342,350T>G locus of the PLCE1 gene were significantly correlated with the growth traits of Ashidan yaks. Among these, individuals with the CC genotype at the g.55,999,531C>T locus showed a significantly higher body length at 6 months old compared to TT individuals, and those with the CT genotype at 12 months old had a significantly higher chest circumference than TT individuals. At the g.342,350T>G locus, the body height of GG genotype individuals at 18 months of age was significantly higher than that of TT genotype individuals and TG genotype individuals. The above findings can be used as theoretical support for the subsequent improvement of Ashidan yak breeding.

This study systematically measured gross hair weight and hair length traits across five body regions of 759 Tianzhu White Yak individuals. The BSL trait exhibited moderate heritability, while the BL trait demonstrated high heritability (h² = 0.450). All other traits showed low heritability. GWASs were conducted using whole-genome resequencing data comprising 22,566,255 high-quality SNP loci. The MLM model identified 519 genome-wide significant loci and 767 chromosome-wide significant loci. Chromosome 6 harbored the highest number of significant SNP loci, while the remaining significant loci were distributed across multiple autosomes. Strong long-range linkage disequilibrium (r² > 0.7) was observed between numerous significant SNPs on chromosome 6 associated with Gw and HL traits. A total of 73 candidate genes were annotated, including FGF5, CFAP299, and PRDM8. Functional enrichment analysis based on the GO and KEGG databases revealed significant enrichment in cytoplasm and the MAPK signaling pathway. Sanger sequencing results revealed that mutations in the FGF5, CFAP299, PRDM8, ANTXR2, and GPHB5 genes significantly affected the Gw, HL, and BSL traits of Tianzhu White Yak (p < 0.01). Linkage disequilibrium analysis indicated strong linkage disequilibrium (r² > 0.6) among Sanger-sequenced SNP loci on the same chromosome. From a biological perspective, multiple candidate genes such as FGF5 and CFAP299 are involved in hair follicle cycle regulation, cell proliferation, and metabolic control, suggesting its potential role in hair follicle development and hair shaft growth. This study identifies candidate loci and genes for gross hair weight and hair length traits in Tianzhu White Yak, contributing to elucidating the genetic mechanisms underlying hair production performance.