Search

Yaks are a rare and unique animal species inhabiting the Qinghai-Tibet Plateau; they are renowned for their remarkable ability to thrive in harsh environments. Milk-derived exosomes, tiny vesicles containing various biological molecules, play crucial roles in numerous pathological and physiological processes, including cell growth, development, and immune regulation. This study delved into the microRNA expression profiles of yak milk-derived exosomes collected from both high- and low-altitude populations using small RNA sequencing. These miRNAs were found to be implicated in pathways associated with mammary gland inflammation, virus infection regulation, and heat stress response. Functional enrichment analyses, utilizing GO and KEGG databases, revealed that the target genes of these differentially expressed miRNAs are enriched in signaling pathways crucial for Th17 cell differentiation and the Ras-MAPK signaling pathway. In conclusion, this research illuminates the adaptive mechanisms of yaks through the differential expression of miRNAs in their milk-derived exosomes across varying environmental conditions. These findings provide a valuable foundation for future investigations into yak resilience and the potential of milk-derived exosomes as tools for disease management and immune modulation.

Single-cell sequencing technology provides apparent advantages in cell population heterogeneity, allowing individuals to better comprehend tissues and organs. Sequencing technology is currently moving beyond the standard transcriptome to the single-cell level, which is likely to bring new insights into the function of breast cells. In this study, we examine the primary cell types involved in breast development, as well as achievements in the study of scRNA-seq in the microenvironment, stressing the finding of novel cell subsets using single-cell approaches and analyzing the problems and solutions to scRNA-seq. Furthermore, we are excited about the field's promising future.

Yak ( Bos grunniens ) is a large ruminant endemic to the Tibetan plateau. The addition of enzyme complexes to feed can significantly improve their growth performance. Therefore, studying the effects of ruminant compound enzyme preparations dosage on yak rumen microorganisms and production performance is crucial to promoting the development of the yak industry. This study aimed to determine the effects of feeding yaks with different doses of ruminant enzyme compounds on the performance, meat quality, and rumen microorganisms of yaks. Three kinds of experimental diets with doses of 0.5 g/kg (LE group), 1 g/kg (ME group), and 2 g/kg (HE group) were selected to determine the growth index, meat quality, serum biochemical indexes, rumen fluid pH and other indexes of the three experimental groups. Metagenomics studies were used to investigate the differences in rumen microbial composition and function among yak groups, and transcriptome sequencing of the longest dorsal muscle was performed to reveal the expression of differential genes among different groups. It was determined that the levels of dietary enzyme complexes significantly affected growth performance, rumen fluid pH, and serum biochemical indices. At the phylum level, the dominant phylum in all three treatment groups was Bacteroidota, Bacillota, Kiritimatiellota, and Pseudomonadota. At the genus level, Prevotella, Methanobrevibacter, Oscillibacter. Fibrobacter showed statistically significant differences in abundance ( p  < 0.05). CAZymes family analysis revealed significant differences in GHs, CTs, and CEs among the three groups. Genome-wide differential gene expression in the longest muscle of the yak back was analyzed by RNA-seq between the three experimental groups. Some DEGs were found to be enriched in the ECM, PI3K-Akt, PPAR, and protein digestion and absorption receptor pathways. Combined metagenomics and transcriptomics analyses revealed that some microorganisms were significantly associated with the genes COL11A1, POSTN, and PTHLH, which are involved in growth metabolism. In summary, this study investigated the effects and interrelationships of ruminant complex enzymes on yak performance, meat quality, and rumen environment. The results of this study provide a scientific basis for adding ruminant enzymes to yaks., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Zhang, La, Ma, Zhandui, Wu, Guo, Yan, Dunzhu and Liang.)

Yaks are crucial to local herders' economy and agriculture. However, several diseases pose a significant threat to the health of yaks and cause substantial economic losses for herders. Therefore, studying the immune indicators and breeding of yaks has become an important task. This study aimed to investigate the association between single nucleotide polymorphisms (SNPs) of the G protein-activated inwardly rectifying K+ channel 2 ( KCNJ6 , GIRK2) gene and yak immune indicators, with the goal of identifying potential candidate molecular markers for yak breeding. In this study, we genotyped 192 healthy adult yaks and detected three SNPs (g163684421 C > T, g163688148 C > T, and g163690745 T > C) in the yak KCNJ6 gene. These SNPs were found to be distributed in the yak population. Subsequently, we performed a linkage disequilibrium analysis and found that the linkage disequilibrium levels of g163684421 C > T and g163690745 T > C were relatively high. Through a correlation analysis of yak KCNJ6 gene SNPs and immune indicators, we found that g163684421 C > T and g163690745 T > C were significantly associated with IgA, IgG, IgM, CRP, HP, IL-2, IL-4, IFN-γ, and TNF-α ( p < 0.05), and the mutation of these SNPs leads to a decrease in yak immune indicators. On the other hand, g163688148 C > T was significantly associated with IgG, IL-4, IFN-γ, TNF-α, IgA, CRP, and HP (p < 0.05), and the mutation of this SNP leads to an increase in yak immune indicators. In conclusion, we identified SNPs associated with yak immune indicators and found that KCNJ6 gene polymorphisms can serve as candidate molecular markers for yak immune indicators. This study provides valuable genetic resources for marker-assisted selection in yak breeding. The results of this study are of great importance for the research on yak immune indicators and marker-assisted selection in yak breeding.

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.

Background: The yak (Bos grunniens) is essential to the livelihoods of Tibetan people on the Qinghai-Tibet Plateau; however, its growth and productivity are constrained by the region's harsh climate and high altitude. Yak skeletal muscle myoblasts, which have evolved to thrive under these challenging conditions, offer a valuable model for investigating muscle development. In this study, we performed transcriptome profiling of yak longissimus dorsi muscle at different growth stages, identifying a key long non-coding RNA, LncRNA-XR&#95;314844 (Lnc-MEG8), with a potential role in muscle development., Results: We developed a novel technique to isolate high-quality yak myoblasts, enabling detailed analysis of Lnc-MEG8. Our results indicated that Lnc-MEG8's subcellular localization varies during muscle cell growth: it is found in both the nucleus and cytoplasm during proliferation but shifts mainly to the cytoplasm during differentiation. Functional experiments showed that Lnc-MEG8 promotes cell proliferation and inhibits differentiation, while its silencing had the opposite effect. Further analysis revealed that both Lnc-MEG8 and the gene RTL1 share miR-22-3p as a common target. Dual-luciferase assays confirmed miR-22-3p directly targets both Lnc-MEG8 and RTL1 mRNA. Co-transfection of Lnc-MEG8 and a miR-22-3p mimic restored RTL1 expression, highlighting Lnc-MEG8's regulatory role. Lnc-MEG8 also counteracts miR-22-3p's suppression of key muscle genes such as MyF5 and MyoG, facilitating myotube formation., Conclusion: These findings demonstrate that the Lnc-MEG8-miR-22-3p-RTL1 axis plays a crucial role in yak muscle development, providing insights that could advance muscle tissue engineering and enhance yak meat quality., Competing Interests: Declarations. Ethics approval and consent to participate: All yaks were handled strictly with good animal practices that complied with the Animal Ethics Procedures and Guidelines of the People's Republic of China. The study was approved by the Animal Administration and Ethics Committee of Lanzhou Institute of Husbandry and Pharmaceutical Sciences of the Chinese Academy of Agricultural Sciences (Permit No. 2019–002). The study was also carried out in compliance with the ARRIVE guidelines [46, 47]. The animals used in this study were from the Datong yak Farm in Qinghai Province. During the study, we obtained informed consent from the animal owners. The owners were provided clear information about the purpose, methods, and anticipated outcomes of the study, and they agreed to include their animals in the research. We received the owner's consent in writing and retained it as part of the ethical approval for this study. Consent for publication: Not applicable. Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

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.

Background: Genomic structural variations (GSVs), notably copy number variations (CNVs), significantly shape genetic diversity and facilitate adaptation in cattle populations. Despite their importance, the genome-wide characterization of CNVs in indigenous Ethiopian cattle breeds-Abigar, Fellata, and Gojjam-Highland remains largely unexplored. In this study, we applied a read-depth approach to whole genome sequencing (WGS) data to conduct the first comprehensive analysis of CNVs in these populations., Results: We identified 3,893 CNV regions (CNVRs) covering 19.15 Mb (0.71% of the cattle genome). These CNVRs ranged from 1.60 kb to 488.0 kb, with an average size of 4.92 kb. These CNVRs included deletions (1713), duplications (1929), and mixed events (251) showing notable differences in distribution among the breeds. Four out of five randomly selected CNVRs were successfully validated using real time polymerase chain reaction (qPCR). Further analyses identified candidate genes associated with high-altitude adaptation (GBE1 and SOD1), heat stress adaptation (HSPA13, DNAJC18, and DNAJC8) and resistance to tick infestations (BoLA and KRT33A). In addition, variance stabilizing transformation (V ST ) statistics highlighted population-specific CNVRs, emphasizing the unique genetic signatures of high-altitude adaptation in the Gojjam-Highland cattle breed. Among the detected CNVRs, 4.93% (192 out of 3,893) overlapped with 520 quantitative traits loci (QTLs) associated with six economically important trait categories suggesting that these CNVRs may significantly contribute to the genetic variation underlying these traits., Conclusions: Our comprehensive analysis reveals significant CNVRs associated with key adaptive traits in Ethiopian cattle breeds highlighting their genetic diversity and resilience. These findings offer valuable insights into the genetic basis of adaptability and can inform sustainable breeding practices and conservation efforts. Future research should prioritize the functional validation of these CNVRs and their integration into breeding programs to enhance traits such as disease resistance and environmental adaptability., Competing Interests: Declarations Ethics approval and consent to participate Blood samples were collected in the study area after obtaining informed consent from each cattle owner. All animals were handled strictly with good animal handling practices, as outlined by the animal ethics procedures of the People’s Republic of China. The study protocol was reviewed and approved by the Animal Administration and Ethics Committee of the Lanzhou Institute of Husbandry and Pharmaceutical Science, CAAS Institutional Review Board (Permit No. SYXK-2014-0002). Consent for publication Not Applicable. Competing interests The authors declare no competing interests., (© 2024. Crown.)

Yak is an important dominant livestock species at high altitude, and the growth performance of yak has obvious differences under different feeding methods. This experiment was conducted to compare the effects of different feeding practices on growth performance and meat quality of yaks through combined transcriptomic and metabolomic analyses. In terms of yak growth performance, compared with traditional grazing, in-house feeding can significantly improve the average daily weight gain, carcass weight and net meat weight of yaks; in terms of yak meat quality, in-house feeding can effectively improve the quality of yak meat. A combined transcriptomic and metabolomic analysis revealed 31 co-enriched pathways, among which arginine metabolism, proline metabolism and glycerophospholipid metabolism may be involved in the development of the longissimus dorsi muscle of yak and the regulation of meat quality-related traits. The experimental results increased our understanding of yak meat quality and provided data materials for subsequent deep excavation of the mechanism of yak meat quality.

Lysozyme like 4 (LYZL4), lysozyme like 6 (LYZL6) and proliferating cell nuclear antigen (PCNA) are implicated in the regulation of testicular function, but there was no research reported available on the expression patterns of LYZL4 , LYZL6 and PCNA genes at different developmental stages of yak testes. In this study, we used the qRT-PCR, western blotting and immunohistochemistry estimated the LYZL4 , LYZL6 and PCNA gene expression and protein lo-calization at different developmental stages of yak testes. The qPCR results showed that the mRNA expression of LYZL4 , LYZL6 and PCNA genes significantly increased with age in the testes of yaks. Western blot results showed that the protein abundance of LYZL4, LYZL6 and PCNA in yak testes was significantly higher after puberty than before puberty. Furthermore, the results of immunohistochemistry indicated that LYZL4, LYZL6 and PCNA may be involved in the regulation of spermatogonia proliferation and Leydig cell function in immature testis. In adult yak testes, LYZL4, LYZL6 and PCNA may involve in the development of round spermatids and primary spermatocytes during testicular development. Our results indicated that LYZL4, LYZL6 and PCNA may be involved in the development of Sertoli cells, Leydig cells and gonocytes in yak testes.

The Meiren yak is an important genetic resource in Gansu Province, China. In this study, we aimed to explore the evolutionary history and population structure of the genetic resource of Meiren yak and to mine the characteristic genes of Meiren yak. We analysed a total of 93 yaks of eight yak breeds based on whole genome resequencing combined with population genomics and used θπ ratio and Fst method to screen the selected sites in the genome region. The results proved that Meiren yak can be used as a potential genetic resource in Gansu Province. The genes in Meiren yak with positive selection in selection signal analysis were subjected to the Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses, which indicated that the genes were related to the adaptability to high altitude and hypoxic environment. By analysing the genetic variation of Meiren yak at the genome-wide level, this study provided a theoretical basis for genetic improvement of Meiren yak and for the development of high-quality yak resources.

This study explored the polymorphism of the leucine aminopeptidase ( LAP3 ) gene and its relationship with milk quality characteristics in Gannan yak. A cohort of 162 Gannan yak was genotyped utilizing the Illumina Yak cGPS 7K BeadChip, and the identified single nucleotide polymorphisms (SNPs) were evaluated for their association with milk protein, casein, lactose, and fat concentrations. The results showed that four SNPs (g.4494G > A, g.5919A > G, g.8033G > C, and g.15,615A > G) in the LAP3 gene exhibited polymorphism with information content values of 0.267, 0.267, 0.293, and 0.114, respectively. All four SNPs were in Hardy-Weinberg equilibrium ( p > 0.05). The g.4494G > A and g.5919A > G SNPs were significantly associated with protein content ( p < 0.05), with homozygous genotypes showing significantly higher protein content than heterozygous genotypes ( p < 0.05). The g.8033G > C SNP was significantly associated with casein content, protein content, non-fat solids, and acidity ( p < 0.05), with the CC genotype having significantly higher casein, protein, and non-fat solids content than the GG and GC genotypes ( p < 0.05). The g.15,615A > G SNP was significantly associated with average fat globule diameter ( p < 0.05). In general, the mutations within the LAP3 gene demonstrated a positive impact on milk quality traits in Gannan yak, with mutated genotypes correlating with enhanced milk quality. These results indicate that the LAP3 gene could be a significant or candidate gene affecting milk quality traits in Gannan yak and offer potential genetic markers for molecular breeding programs in this species.

Yaks are crucial genetic resources in the Tibetan Plateau and surrounding regions. Throughout the long process of domestication, natural and artificial selection pressures have enabled yaks to demonstrate adaptive characteristics to the environment in terms of physiological structure and genetic molecules, but no systematic cell analysis has been carried out on this phenomenon of yaks. Here, the population structure and genetic diversity of yak were studied by WGRS, and the genes related to yak adaptability were excavated. Combined with scRNA-seq method, the transcription map of yak lung tissue and skin tissue was constructed, which provided a new comprehensive insight into yak adaptability. The analysis of yak population structure showed that there was obvious genetic differentiation between TZ &#95; yak and other seven yak populations, while there was significant genetic exchange between PL &#95; yak and SB &#95; yak at high altitude. WGRS and scRNA-seq analysis revealed that the gene HIF1A related to high altitude adaptation was expressed in various cell types, while EPAS1 was predominantly expressed in epithelial and endothelial cells of yak lung tissue. Endothelial cells play a critical role in hypoxia-adapted VEGF signaling, which correlates closely with the high expression of KDR and VEGFA genes in endothelial cells and monocytes. Furthermore, in the selection signal of High &#95; yak vs Low &#95; yak, 19.8 % of the genes overlapped with the genes screened by skin scRNA-seq, including genes related to coat color such as RORA, BNC2, and KIT. Notably, BNC2 is a gene associated with melanin deposition and shows high expression levels in HS cells. Additionally, GRN in melanocytes and SORT1 in IRS play an important role in cell communication between melanocytes and IRS. These findings offer new insights into the natural polymorphism of yaks and provide a valuable reference for future research on high-altitude mammals, and potentially even human genetics., Competing Interests: Declaration of competing interest The authors declare no conflict of interest., (Copyright © 2024. Published by Elsevier Inc.)

Energy feed can provide animals with balanced nutrition, thereby enhancing their growth performance. This study aimed to evaluate the effects of dietary energy levels on the growth performance, serum metabolites, and meat quality of Jersey cattle-yaks. A total of 24 male Jersey cattle-yaks were randomly divided into three groups. Each group was fed diets with metabolizable energy levels of 8.21 MJ/kg (LE), 9.50 MJ/kg (ME), and 10.65 MJ/kg (HE), respectively. The HE and ME groups showed significantly higher final body weight, average daily gain (ADG), and feed efficiency compared to the LE group ( p < 0.05). The glucose (GLU) and total cholesterol (TC) concentrations were significantly increased in the serum of the ME and HE groups ( p < 0.05). The low-density lipoprotein cholesterol (LDL-C) and alanine aminotransferase (ALT) levels were significantly higher in the serum of the HE group than in the ME group ( p < 0.05). Blood urea nitrogen (BUN) levels exhibited a significant decrease with increasing metabolizable energy levels in the diet ( p < 0.05). Increasing dietary energy levels enhances the eye muscle area and intramuscular fat content of Jersey cattle-yaks ( p < 0.05), with no effect on pH45 min, pH24 h, and shear force. In the HE group, the levels of heneicosanoic acid (C21:0), palmitoleic acid (C16:1), elaidic acid (C18:1n9t), and eicosadienoic acid (C20:2n6) were notably elevated ( p < 0.05) when compared to the LE group. We concluded that a higher dietary energy level enhanced the growth performance and meat quality traits of male Jersey cattle-yaks.

Male reproductive health is largely determined already in the early development of the testis. Although much work has been carried out to study the mechanisms of testicular development and spermatogenesis, there was previously no information on the differences in the protein composition of yak testicles during early development. In this study, the protein profiles in the testicles of 6- (M6), 18- (M18), and 30-month-old (M30) yaks were comparatively analyzed using TMT proteomics. A total of 5521 proteins were identified, with 13, 1295, and 1397 differentially expressed proteins (DEPs) in 30- vs. 18-, 18- vs. 6-, and 30- vs. 6-month-old testes, respectively. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that DEPs were mainly involved in signaling pathways related to testicular development and spermatogenesis, including the MAPK, PI3K-Akt, Wnt, mTOR, TGF-β, and AMPK signaling pathways. Furthermore, we also identified eight potential proteins (TEX101, PDCL2, SYCP2, SYCP3, COL1A1, COL1A2, ADAM10, and ATF1) that may be related to the testicular development and spermatogenesis of yaks. This study may provide new insights into the molecular mechanisms of the testicular development and spermatogenesis of yaks.

Cattle-yak is a hybrid offspring resulting from the crossbreeding of yak and cattle, and it exhibits substantial heterosis in production performance. However, male sterility in cattle-yak remains a concern. Reports suggest that noncoding RNAs are involved in the regulation of spermatogenesis. Therefore, in this study, we comprehensively compared testicular transcription profiles among cattle, yak, and cattle-yak. Numerous differentially expressed genes (DEGs), differentially expressed circRNAs (DECs), and differentially expressed miRNAs (DEMs) were identified in the intersection of two comparison groups, namely cattle versus cattle-yak and yak versus cattle-yak, with the number of DEGs, DECs, and DEMs being 4968, 360, and 59, respectively. The DEGs in cattle-yaks, cattle, and yaks were mainly associated with spermatogenesis, male gamete generation, and sexual reproduction. Concurrently, GO and KEGG analyses indicated that DEC host genes and DEM source genes were involved in the regulation of spermatogenesis. The construction of a potential competing endogenous RNA network revealed that some differentially expressed noncoding RNAs may be involved in regulating the expression of genes related to testicular spermatogenesis, including miR-423-5p, miR-449b, miR-34b/c, and miR-15b, as well as previously unreported miR-6123 and miR-1306, along with various miRNA-circRNA interaction pairs. This study serves as a valuable reference for further investigations into the mechanisms underlying male sterility in cattle-yaks., Competing Interests: Declaration of competing interest The authors declare that they have no competing interests., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Yak milk, known as the "liquid gold", is a nutritious food with extensive consumption. Compared with cow milk, yak milk contains higher levels of nutrients such as dry matter, milk fat, and milk protein, which demonstrates great potential for exploitation and utilization. Protein kinase cGMP-dependent 1 ( PRKG1 ) is an important functional molecule in the cGMP signaling pathway, and its significant influence on milk fatty acids has been discovered. The aim of this study is to explore the correlation between single nucleotide polymorphisms (SNPs) in the PRKG1 gene and the quality traits of Gannan yak milk in order to identify candidate molecular markers for Gannan yak breeding. In this study, genotyping was performed on 172 healthy, 4-5-year-old lactating Gannan yaks with similar body types, naturally grazed, and two to three parity. Three SNPs (g.404195C>T, g.404213C>T, and g.760138T>C) were detected in the PRKG1 gene of Gannan yaks, which were uniformly distributed in the yak population. Linkage disequilibrium analysis was conducted, revealing complete linkage disequilibrium between g.404195C>T and g.404213C>T. After conducting a correlation analysis between SNPs in the PRKG1 gene and milk quality in Gannan yaks, we found that PRKG1 SNPs significantly increased the content of casein, protein, and SNFs in yak milk. Among them, the TT homozygous genotype at the PRKG1 g.404195C>T loci exhibited higher casein and protein contents compared to the CC and CT genotypes ( p < 0.05). The SNP g.760138T>C locus was associated with casein, protein, SNFs, and TS traits ( p < 0.05). The CC genotype had higher casein and protein contents than the TT and TA genotypes ( p < 0.05). However, there were no significant differences in milk fat, lactose, and acidity among the three genotypes ( p > 0.05). In summary, PRKG1 gene polymorphism can serve as a candidate molecular marker for improving milk quality in Gannan yaks.

Indigenous animal genetic resources play a crucial role in preserving global genetic diversity and supporting the livelihoods of millions of people. In Ethiopia, the majority of the cattle population consists of indigenous breeds. Understanding the genetic architecture of these cattle breeds is essential for effective management and conservation efforts. In this study, we sequenced DNA samples from 70 animals from seven indigenous cattle breeds, generating about two terabytes of pair-end reads with an average coverage of 14X. The sequencing data were pre-processed and mapped to the cattle reference genome (ARS-UCD1.2) with an alignment rate of 99.2%. Finally, the variant calling process produced approximately 35 million high-quality SNPs. These data provide a deeper understanding of the genetic landscape, facilitate the identification of causal mutations, and enable the exploration of evolutionary patterns to assist cattle improvement and sustainable utilization, particularly in the face of unpredictable climate changes., (© 2024. The Author(s).)

In this study, our primary aim was to explore the genomic landscape of Barka cattle, a breed recognized for high milk production in a semi-arid environment, by focusing on genes with known roles in milk production traits. We employed genome-wide analysis and three selective sweep detection methods ( ZF ST , θπ ratio, and ZHp ) to identify candidate genes associated with milk production and composition traits. Notably, ACAA1 , P4HTM , and SLC4A4 were consistently identified by all methods. Functional annotation highlighted their roles in crucial biological processes such as fatty acid metabolism, mammary gland development, and milk protein synthesis. These findings contribute to understanding the genetic basis of milk production in Barka cattle, presenting opportunities for enhancing dairy cattle production in tropical climates. Further validation through genome-wide association studies and transcriptomic analyses is essential to fully exploit these candidate genes for selective breeding and genetic improvement in tropical dairy cattle.

Alternative splicing is a prevalent phenomenon in testicular tissues. Due to the low assembly accuracy of short-read RNA sequencing technology in analyzing post-transcriptional regulatory events, full-length (FL) transcript sequencing is highly demanded to accurately determine FL splicing variants. In this study, we performed FL transcriptome sequencing of testicular tissues from 0.5, 1.5, 2.5, and 4-year-old yaks and 4-year-old cattle-yaks using Oxford Nanopore Technologies. The obtained sequencing data were predicted to have 47,185 open reading frames (ORFs), including 26,630 complete ORFs, detected 7645 fusion transcripts, 15,355 alternative splicing events, 25,798 simple sequence repeats, 7628 transcription factors, and 35,503 long non-coding RNAs. A total of 40,038 novel transcripts were obtained from the sequencing data, and the proportion was almost close to the number of known transcripts identified. Structural analysis and functional annotation of these novel transcripts resulted in the successful annotation of 9568 transcripts, with the highest and lowest annotation numbers in the Nr and KOG databases, respectively. Weighted gene co-expression network analysis revealed the key regulatory pathways and hub genes at various stages of yak testicular development. Our findings enhance our comprehension of transcriptome complexity, contribute to genome annotation refinement, and provide foundational data for further investigations into male sterility in cattle-yaks., Competing Interests: Declaration of competing interest The authors declare that they have no competing interests., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

A single-nucleotide polymorphism (SNP) is a genome-level trait that arises from a variation in a single nucleotide, leading to diversity in DNA sequences. SNP screening is commonly used to provide candidate genes for yak breeding efforts. Integrin Subunit Alpha 9 ( ITGA9 ) is an integrin protein. It plays an important role in cell adhesion, signalling, and other processes. The aim of this study was to discuss the association between genetic polymorphisms in the ITGA9 gene and milk quality traits and to identify potential molecular marker loci for yak breeding quality. We genotyped 162 yaks using an Illumina Yak cGPS 7K liquid chip and identified the presence of polymorphisms at nine SNP loci in the ITGA9 gene of yaks. The results showed that the mutant genotypes in the loci g.285,808T>A, g.306,600T>C, and g.315,413C>T were positively correlated with the contents of casein, protein, total solids (TS), and solid nonfat (SNF) in yak milk. In other loci, heterozygous genotypes had a positive correlation with nutrient content in yak milk. Then, two ITGA9 haplotype blocks were constructed based on linkage disequilibrium, which facilitated a more accurate screening of ITGA9 as a candidate gene for yak milk quality improvement. In conclusion, we identified SNPs and haplotype blocks related to yak milk quality traits and provided genetic resources for marker-assisted selection in yak breeding.

Background: The hair follicle development process is regulated by sophisticated genes and signaling networks, and the hair grows from the hair follicle. The Tianzhu white yak population exhibits differences in hair length, especially on the forehead and shoulder region. However, the genetic mechanism is still unclear. Isoform sequencing (Iso-seq) technology with advantages in long reads sequencing. Hence, we combined the Iso-seq and RNA-seq methods to investigate the transcript complexity and difference between long-haired yak (LHY) and normal-haired yak (NHY)., Results: The hair length measurement result showed a significant difference between LHY and NHY on the forehead and the shoulder (P-value < 0.001). The skin samples from the forehead and the shoulder of LHY and NHY were pooled for isoform sequencing (Iso-seq). We obtained numerous long transcripts, including novel isoforms, long non-coding RNA, alternative splicing events, and alternative polyadenylation events. Combined with RNA-seq data, we performed differential isoforms (DEIs) analysis between LHY and NHY. We found that some hair follicle and skin development-related DEIs, like BMP4, KRT2, IGF2R, and COL1A2 in the forehead skin; BMP1, KRT1, FGF5, COL2A1, and IGFBP5 in the shoulder skin. Enrichment analysis revealed that DEIs in both two comparable groups significantly participated in skin and hair follicle development-related pathways, like ECM-receptor interaction, focal adhesion, and PI3K-Akt signaling pathways. The results indicated that the hair follicle development of Tianzhu white yak may influence the hair length difference. Besides, the protein-protein interaction (PPI) network of DEIs showed COL2A1 and COL3A1 exhibited a high degree of centrality, and these two genes were suggested as potential candidates for the hair length growth of Tianzhu white yak., Conclusions: The results provided a comprehensive analysis of the transcriptome complexity and identified differential transcripts that enhance our understanding of the molecular mechanisms underlying the variation in hair length growth in Tianzhu white yak., (© 2024. The Author(s).)

Background: There is no consensus as to the origin of the domestic yak (Bos grunniens). Previous studies on yak mitochondria mainly focused on mitochondrial displacement loop (D-loop), a region with low phylogenetic resolution. Here, we analyzed the entire mitochondrial genomes of 509 yaks to obtain greater phylogenetic resolution and a comprehensive picture of geographical diversity., Results: A total of 278 haplotypes were defined in 509 yaks from 21 yak breeds. Among them, 28 haplotypes were shared by different varieties, and 250 haplotypes were unique to specific varieties. The overall haplotype diversity and nucleotide diversity of yak were 0.979 ± 0.0039 and 0.00237 ± 0.00076, respectively. Phylogenetic tree and network analysis showed that yak had three highly differentiated genetic branches with high support rate. The differentiation time of clades I and II were about 0.4328 Ma, and the differentiation time of clades (I and II) and III were 0.5654 Ma. Yushu yak is shared by all haplogroups. Most (94.70%) of the genetic variation occurred within populations, and only 5.30% of the genetic variation occurred between populations. The classification showed that yaks and wild yaks were first clustered together, and yaks were clustered with American bison as a whole. Altitude had the highest impact on the distribution of yaks., Conclusions: Yaks have high genetic diversity and yak populations have experienced population expansion and lack obvious phylogeographic structure. During the glacial period, yaks had at least three or more glacial refugia., (© 2024. The Author(s).)

The yak is a unique species of livestock found in the Qinghai-Tibet Plateau and its surrounding areas. Due to factors such as late sexual maturity and a low rate of estrus, its reproductive efficiency is relatively low. The process of estrus synchronization in yaks plays a crucial role in enhancing their reproductive success and ensuring the continuation of their species. In order to clarify the characteristics of the serum metabolites of yak estrus synchronization, the yaks with inactive ovaries were compared with the estrus synchronization yaks. In this study, yaks were divided into the inactive ovaries group (IO), gonarelin-induced yak estrus group (GnRH), and chloprostenol sodium-induced yak estrus group (PGF). After the completion of the estrus synchronization treatment, blood samples were collected from the jugular veins of the non-estrus yaks in the control group and the yaks with obvious estrus characteristics in the GnRH and PGF groups. Metabolites were detected by ultra-high performance liquid chromatography-mass spectrometry, and differential metabolites were screened by multivariate statistical analysis. The results showed that a total of 70 significant differential metabolites were screened and identified in the GnRH vs. IO group, and 77 significant differential metabolites were screened and identified in the PGF vs. IO group. Compared with non-estrus yaks, 36 common significant differential metabolites were screened out after the induction of yak estrus by gonarelin (GnRH) and cloprostenol sodium (PGF), which were significantly enriched in signaling pathways such as the beta oxidation of very long chain fatty acids, bile acid biosynthesis, oxidation of branched chain fatty acids, steroidogenesis, steroid biosynthesis, and arginine and proline metabolism. This study analyzed the effects of gonadotropin releasing hormone (GnRH) and prostaglandin F (PGF) on the reproductive performance of yaks treated with estrus synchronization, which provides a theoretical basis for the optimization and application of yak estrus synchronization technology and promotes the healthy development of the yak industry.

Protein kinase D1 (PRKD1) functions primarily in normal mammary cells, and the potassium voltage-gated channel subfamily Q member 3 ( KCNQ3 ) gene plays an important role in controlling membrane potential and neuronal excitability, it has been found that this particular gene is linked to the percentage of milk fat in dairy cows. The purpose of this study was to investigate the relationship between nucleotide polymorphisms (SNPs) of PRKD1 and KCNQ3 genes and the milk quality of Gannan yak and to find molecular marker sites that may be used for milk quality breeding of Gannan yak. Three new SNPs were detected in the PRKD1 (g.283,619T>C, g.283,659C>A) and KCNQ3 gene (g.133,741T>C) of 172 Gannan lactating female yaks by Illumina yak c GPS 7K liquid-phase microarray technology. Milk composition was analyzed using a MilkoScan TM milk composition analyzer. We found that the mutations of these three loci significantly improved the lactose, milk fat, casein, protein, non-fat milk solid (SNF) content and acidity of Gannan yaks. The lactose content of the TC heterozygous genotype population at g.283,619T>C locus was significantly higher than that of the TT wild-type population ( p < 0.05); the milk fat content of the CA heterozygous genotype population at g.283,659C>A locus was significantly higher than that of the CC wild-type and AA mutant populations ( p < 0.05); the casein, protein and acidity of the CC mutant and TC heterozygous groups at the g.133,741T>C locus were significantly higher than those of the wild type ( p < 0.05), and the SNF of the TC heterozygous group was significantly higher than that of the mutant group ( p < 0.05). The results showed that PRKD1 and KCNQ3 genes could be used as candidate genes affecting the milk traits of Gannan yak.

Cattle-yak, which is the hybrid F 1 generation of cattle and yak, demonstrates better production performance compared to yak. However, there is limited research on the molecular mechanisms responsible for the muscle development of cattle-yak. To address this knowledge gap, a comprehensive transcriptomic survey of the longissimus dorsi muscle in cattle-yak was conducted. Three transcript types, namely lncRNAs, miRNAs, and circRNAs, along with protein-coding genes were characterized at two developmental stages (6 m, 18 m) of cattle-yak. The results revealed significant enrichment of these transcripts into pathways related to myoblast differentiation and muscle development signaling. Additionally, the study identified the TCONS00024465/circHIPK3-bta-miR-499-ADAMTS6 regulatory network, which may play a crucial role in the muscle development of cattle-yak by combining the transcriptome data of yak and constructing the ceRNA co-expression network. HEK 293 T cells were used to validate that TCONS00024465 and circHIPK3 are located upstream of bta-miR-499, and can competitively bind to bta-miR-499 as ceRNA. The study also verified that ADAMTS6 regulates skeletal muscle development by inhibiting myoblast proliferation, promoting myoblast differentiation, and positively regulating the apoptosis of myoblasts. Taken together, this study provides new insights into the advantages of cattle-yak production performance and offers a molecular basis for further research on muscle development., Competing Interests: Declaration of competing interest The authors declare no conflict of interest., (Copyright © 2024. Published by Elsevier B.V.)

In mammals, epigenetic modifications involving DNA methylation are necessary for the completion of the cell differentiation process. However, the global DNA methylation landscape and its dynamics during yak adipocyte differentiation remain unexplored. Here, we performed whole-genome bisulfite sequencing (WGBS) to asses DNA methylation in yak preadipocytes and adipocytes, combining these results with those of our previous studies on changes in chromatin accessibility and gene expression during yak adipogenesis. The results showed that CG methylation levels were lower in promoter than in exon and intron, and gradually decreasing from the distal regions to transcription start site (TSS). There was a significant negative correlation between CG methylation levels located in promoter and gene expression levels. The 46 genes shared by CG differentially methylated regions (DMRs) and differential chromatin accessibility were significantly enriched in Hedgehog and PI3K-Akt signaling pathways. ATAC-seq peaks with high chromatin accessibility located in both promoter (≤ 2 kb from TSS) and distal (> 2 kb from TSS) regions corresponded to low methylation levels. Taken together, these findings demonstrated that DNA methylation and its interactions with chromatin accessibility regulate gene expression during yak adipocyte differentiation, contributing to the understanding of mechanisms of various epigenetic modifications and their interactions in adipogenesis., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier B.V.)

Background: Yaks (Bos grunniens), prized for their ability to thrive in high-altitude environments, are indispensable livestock in the plateau region. Modifying their feeding systems holds significant promise for improving their growth and meat quality. Tenderness, a key determinant of yak meat quality and consumer appeal, is demonstrably influenced by dietary regimen. Indoor feeding regimes have been shown to enhance tenderness by lowering shear stress and optimizing pH values. CircRNAs, well-known modulators of circulatory function, also play a crucial role in skeletal muscle development across various animal species. However, their functional significance in yak skeletal muscle remains largely unexplored., Results: In this study, we identified a total of 5,534 circRNAs within the longissimus dorsi muscle, and we found 51 differentially expressed circRNAs (20 up-regulated and 31 down-regulated) between the two feeding groups. Constructing a comprehensive ceRNA network illuminated intricate regulatory mechanisms, with PGP and circRNA&#95;0617 converging on bta-miR-2285q, mirrored by KLF15/circRNA&#95;0345/bta-miR-20b and CTSF/circRNA&#95;0348/bta-miR-146a. These findings shed light on the potential of circRNAs to influence yak muscle development and meat quality, offering valuable insights for future research., Conclusions: This investigation unraveled a complex interaction network between circRNAs、mRNAs and miRNAs in yak skeletal muscle. We further elucidated the target genes regulated by these target genes within the network, offering valuable insights into the potential regulatory mechanisms governing muscle development and meat quality-related traits in yaks., (© 2024. The Author(s).)

Yak meat is nutritionally superior to beef cattle but has a low fat content and is slow-growing. The liver plays a crucial role in lipid metabolism, and in order to determine whether different feeding modes affect lipid metabolism in yaks and how it is regulated, we employed RNA sequencing (RNA-seq) technology to analyze the genome-wide differential gene expression in the liver of yaks maintained under different raising systems. A total of 1663 differentially expressed genes (DEGs) were identified (|log2FC| ≥ 0 and p -value ≤ 0.05), including 698 down-regulated and 965 up-regulated genes. According to gene ontology (GO) and KEGG enrichment analyses, these DEGs were significantly enriched in 13 GO terms and 26 pathways ( p < 0.05). Some DEGs were enriched in fatty acid degradation, PPAR, PI3K-Akt, and ECM receptor pathways, which are associated with lipid metabolism. A total of 16 genes are well known to be related to lipid metabolism (e.g., APOA1 , FABP1 , EHHADH , FADS2 , SLC27A5 , ACADM , CPT1B , ACOX2 , HMGCS2 , PLIN5 , ACAA1 , IGF1 , FGFR4 , ALDH9A1 , ECHS1 , LAMA2 ). A total of 11 of the above genes were significantly enriched in the PPAR signaling pathway. The reliability of the transcriptomic data was verified using qRT-PCR. Our findings provide new insights into the mechanisms regulating yak meat quality. It shows that fattening improves the expression of genes that regulate lipid deposition in yaks and enhances meat quality. This finding will contribute to a better understanding of the various factors that determine yak meat quality and help develop strategies to improve yield and quality.

Jersey-yak is a hybrid offspring of Jersey cattle and yak ( Bos grunniens ). Changing the feeding system of Jersey-yak can significantly improve its growth performance. In this study, tandem mass tag (TMT) proteomics technology was used to determine the differentially expressed proteins (DEPs) of the longissimus lumborum (LL) muscle of Jersey-yak fed different protein levels of diet. The results showed that compared with the traditional grazing feeding, the growth performance of Jersey-yaks was significantly improved by crude protein supplementation after grazing. A total of 3368 proteins were detected in these muscle samples, of which 3365 were quantified. A total of 434 DEPs were identified. Through analyses, it was found that some pathways related to muscle growth and development were significantly enriched, such as Rap1 signaling pathway, mTOR signaling pathway, and TGF-beta signaling pathway. A number of DEPs enriched in these pathways are related to muscle cell development, differentiation, and muscle development, including integrin subunit alpha 7 (ITGA7), myosin heavy chain 8 (MYH8), and collagen type XII alpha 1 chain (COL12A1). In conclusion, the results of this study provide insights into the proteomics of different feeding patterns of Jersey-yak, providing a stronger basis for further understanding the biological mechanism of hybrid varieties.

Yak ( Bos grunniens ) is a unique large ruminant species in the Qinghai-Tibetan Plateau (QTP). Changing the energy levels of their rations can significantly improve their growth performance. Therefore, studying the effects of dietary energy levels on the rumen microflora and metabolites of yak is crucial for enhancing the development of the yak industry. Currently, there is a lack of understanding regarding the impact of feeding energy diets on rumen fermentation parameters, microbial functions, and metabolites. This study was designed to determine the appropriate energy level for feeding yak. Three test diets with metabolizable energy levels of 7.57 MJ/kg, 9.44 MJ/kg, and 11.9 MJ/kg were used and the concentration of volatile fatty acids (VFA) in rumen fluid was measured. The microbial communities, functions, and metabolites in yaks were studied by 16S rRNA sequencing, metagenome, and LC-MS non-targeted metabolomics to investigate the relationships among rumen fermentation parameters, microbial diversity, and metabolites. Ration energy levels significantly affect total VFA, acetate, propionate, butyrate, iso -valerate, valerate, and acetate/propionate ( p < 0.05). At the phylum level, the dominant phyla in all three treatment groups were Bacteroidota, Firmicutes, and Actinobacteriota. At the genus level, the abundance of the unclassified&#95;o&#95;&#95;Bacteroidales , norank&#95;f&#95;Muribaculaceae , Lachnospiraceae&#95;NK4A136&#95;group , and Family &#95;XIII&#95;AD3011&#95;group showed significant differences ( p < 0.05) and were significantly correlated with differential metabolites screened for phosphatidylcholine [PC(16:0/0:0), PC(18:3/0:0)], uridine 3'-monophosphate, and adenosine monophosphate, etc. CAZymes family analysis showed that GHs and CEs differed significantly among the three groups. In addition, differential metabolites were mainly enriched in the pathways of lipid metabolism, nucleotide metabolism, and biosynthesis of other secondary metabolites, and the concentrations of differential metabolites were correlated with microbial abundance. In summary, this study analyzed the effects of ration energy levels on rumen microorganisms and metabolites of yaks and their relationships. The results provided a scientific basis for the selection of dietary energy for yaks in the house feeding period in the future., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Ma, La, Yang, Dai, Zhang, Zhang, Jin, Ma, Guo, Chu, Yan, Zhang and Liang.)

The variety of species can be efficiently increased by interspecific hybridization. However, because the males in the hybrid progeny are usually sterile, this heterosis cannot be employed when other cattle and yaks are hybridized. While some system-level studies have sought to explore the etiological basis for male cattle-yak sterility, no systematic cellular analyses of this phenomenon have yet been performed. Here, single-cell RNA sequencing and UPHLC-MS/MS targeted metabolomics methods were used to study the differences in testicular tissue between 4-year-old male yak and 4-year-old male cattle-yak, providing new and comprehensive insights into the causes of male cattle-yak sterility. Cattle-yak testes samples detected 6 somatic cell types and one mixed germ cell type. Comparisons of these cell types revealed the more significant differences in Sertoli cells (SCs) and [Leydig cells and myoid cells (LCs&#95;MCs)] between yak and cattle-yak samples compared to other somatic cell clusters. Even though the LCs and MCs from yaks and cattle-yaks were derived from the differentiation of the same progenitor cells, a high degree of overlap between LCs and MCs was observed in yak samples. Still, only a small overlap between LCs and MCs was observed in cattle-yak samples. Functional enrichment analyses revealed that genes down-regulated in cattle-yak SCs were primarily enriched in biological activity, whereas up-regulated genes in these cells were enriched for apoptotic activity. Furthermore, the genes of up-regulated in LCs&#95;MCs of cattle-yak were significantly enriched in enzyme inhibitor and molecular function inhibitor activity. On the other hand, the genes of down-regulated in these cells were enriched for signal receptor binding, molecular function regulation, positive regulation of biological processes, and regulation of cell communication activity. The most significant annotated differences between yak and cattle-yak LCs&#95;MCs were associated with cell-to-cell communication. While yak LCs&#95;MCs regulated spermatogenic cells at spermatogonia, spermatocyte, and spermatid levels, no such relationships were found between cattle-yak LCs&#95;MCs and germ cells. This may suggest that the somatic niche in male cattle-yak testes is a microenvironment that is ultimately not favorable for spermatogenesis., Competing Interests: Declaration of competing interest The authors declare that they have no competing interests., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

The morphogenesis of hair follicle structure is accompanied by the differentiation of skin tissue. Mammalian coats are produced by hair follicles. The formation of hair follicles requires signal transmission between the epidermis and dermis. However, knowledge of the transcriptional regulatory mechanism is still lacking. We used single-cell RNA sequencing to obtain 26,573 single cells from the scapular skin of yaks at hair follicle telogen and anagen stages. With the help of known reference marker genes, 11 main cell types were identified. In addition, we further analyzed the DP cell and dermal fibroblast lineages, drew a single-cell map of the DP cell and dermal fibroblast lineages, and elaborated the key genes, signals, and functions involved in cell fate decision making. The results of this study provide a very valuable resource for the analysis of the heterogeneity of DP cells and dermal fibroblasts in the skin and provide a powerful theoretical reference for further exploring the diversity of hair follicle cell types and hair follicle morphogenesis.

Coiled-coil serine-rich protein 1 ( CCSER 1 ) gene is a regulatory protein gene. This gene has been reported to be associated with various economic traits in large mammals in recent years. The aim of this study was to investigate the association between CCSER1 gene single nucleotide polymorphisms (SNPs) and Gannan yaks and to identify potential molecular marker loci for breeding milk quality in Gannan yaks. We genotyped 172 Gannan yaks using Illumina Yak c GPS 7K liquid microarrays and analyzed the correlation between the three SNPs loci of the CCSER1 gene and the milk qualities of Gannan yaks, including milk fat, protein and casein. It was found that mutations at the g.183,843A>G, g.222,717C>G and g.388,723G>T loci all affected the fat, protein, casein and lactose traits of Gannan yak milk to varying extents, and that the milk quality of individuals with mutant phenotypes was significantly improved. Among them, the milk fat content of AG heterozygous genotype population at g.183,843A>G locus was significantly higher than that of AA and GG genotype populations ( p < 0.05); the casein and protein content of mutant GG and CG genotype populations at g.222,717C>G locus was significantly higher than that of wild-type CC genotype population ( p < 0.05); and the g.388,723G>T locus of the casein and protein contents of the mutant TT genotype population were significantly higher ( p < 0.05) than those of the wild-type GG genotype population. These results provide potential molecular marker sites for Gannan yak breeding.

Sperm cryodamage caused by cryopreservation limits the use of frozen yak spermatozoa in artificial insemination (AI). However, the proteomic changes involved in the cryodamage of yak spermatozoa have not been investigated to date. Therefore, this study aimed to identify proteins related to freezing tolerance. Tandem mass tag (TMT) were used in combination with liquid chromatography-tandem mass spectrometry (LC-MS/MS) for identifying differentially expressed proteins (DEPs) between high-motility (HM) and low-motility (LM) frozen-thawed yak spermatozoa. A total of 116 DEPs were identified (>1.5-fold, P < 0.05); of which, 104 proteins were upregulated in HM spermatozoa and 12 proteins were upregulated in LM spermatozoa. The results of functional annotation analysis revealed that the DEPs were mainly involved in metabolic processes. A total of 20 DEPs that were abundantly expressed in HM spermatozoa were strongly associated with carbohydrate metabolism. The results of KEGG analysis revealed that the DEPs were enriched in glycolysis/gluconeogenesis, PPAR signaling pathway, and Ras signaling pathway. In addition, many antioxidant enzymes such as superoxide dismutase (SOD1), peroxiredoxin-6 (PRDX6), and Parkinson disease protein 7 (PARK7) were upregulated in HM spermatozoa, suggesting that these enzymes affect the motility of spermatozoa by regulating the levels of reactive oxygen species (ROS) in frozen-thawed spermatozoa. Altogether, the findings of this study elucidate the mechanisms through which cryopreservation affects the movement of yak spermatozoa and offer a novel basis for refining freezing techniques and modifying cryopreservation extender components., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Yak has been subject to natural selection, human domestication and interspecific introgression during its evolution. However, genetic variants favored by each of these processes have not been distinguished previously. We constructed a graph-genome for 47 genomes of 7 cross-fertile bovine species. This allowed detection of 57,432 high-resolution structural variants (SVs) within and across the species, which were genotyped in 386 individuals. We distinguished the evolutionary origins of diverse SVs in domestic yaks by phylogenetic analyses. We further identified 334 genes overlapping with SVs in domestic yaks that bore potential signals of selection from wild yaks, plus an additional 686 genes introgressed from cattle. Nearly 90% of the domestic yaks were introgressed by cattle. Introgression of an SV spanning the KIT gene triggered the breeding of white domestic yaks. We validated a significant association of the selected stratified SVs with gene expression, which contributes to phenotypic variations. Our results highlight that SVs of different origins contribute to the phenotypic diversity of domestic yaks., (© 2023. Springer Nature Limited.)

The changes in semen and cryodamage after the cryopreservation process negatively affect sperm function and motility. However, possible proteomic alterations of yak semen during cryopreservation have not yet been achieved. In this study, we compared proteomes of fresh and frozen thawed yak sperm using iTRAQ combined with LC-MS/MS proteome approach. Totally, 2064 proteins were quantitatively identified, including 161 in fresh sperm that showed significant differences compared to frozen thawed sperm. According to the Gene ontology (GO) enrichment analysis, differentially expressed proteins (DEPs) are predominantly associated with spermatogenesis, tricarboxylic acid cycle, ATP synthesis, and differentiation biological process. Furthermore, the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that DEPs were mainly involved in metabolic pathways related to pyruvate metabolism, carbon metabolism, glycolysis/gluconeogenesis, together with the citrate (TCA) cycle. In the analysis of the protein-protein interaction (PPI) network, 15 potential proteins (PDHB, DLAT, PDHA2, PGK1, TP5C1, etc.) that could be related to the sperm quality of the yaks were obtained. Furthermore, 6 DEPs were validated by parallel reaction monitoring (PRM), confirming that the iTRAQ data were reliable. These results indicate that cryopreservation alters the proteome of yak sperm, which is possibly related to cryodamage and fertilization ability., Competing Interests: Declaration of competing interest The authors declare that there is no conflict of interest regarding the publication of this paper., (Copyright © 2023. Published by Elsevier B.V.)

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.

Understanding animal growth plays an important role in improving animal genetics and breeding. In order to explore the early growth and development law of Ashidan yak, the body weight (BW), wither height (WH), body oblique length (BL) and chest girth (CG) of 260 female Ashidan yaks were measured. These individuals grew under grazing conditions, and growth traits were measured at 6, 12, 18 and 30 months of age. Then the absolute growth and relative growth of Ashidan yak were calculated, and five nonlinear models (Logistic model, Gompertz model, Brody model, von Bertalanffy model and Richards model) were used to fit the growth curve of Ashidan yak. The fitting effect of the model was evaluated according to MSE, AIC and BIC. The results showed that the growth rate of Ashidan yak was the fastest from 12 to 18 months old, and the growth was slow or even stagnant from 6 to 12 months old. The AIC and BIC values of the Richards model were the lowest among the five models, with an AIC value of 4543.98 and a BIC value of 4563.19. The Richards model estimated body weight at 155.642 kg. In summary, the growth rate of female Ashidan yak changes with the seasons, growing faster in warm seasons and slower in cold seasons. Richards model is the best model to describe the growth curve of female Ashidan yak in five nonlinear models.

Introduction: The crude protein level in the diet will affect the fermentation parameters, microflora, and metabolites in the rumen of ruminants. It is of great significance to study the effect of crude protein levels in supplementary diet on microbial community and metabolites for improving animal growth performance. At present, the effects of crude protein level in supplementary diet on rumen fermentation parameters, microbial community, and metabolites of Jersey-Yak (JY) are still unclear., Methods: The purpose of this experiment was to study the appropriate crude protein level in the diet of JY. The rumen fermentation indexes (volatile fatty acids and pH) were determined by supplementary diets with crude protein levels of 15.16 and 17.90%, respectively, and the microbial community and metabolites of JYs were analyzed by non-target metabonomics and metagenome sequencing technology, and the changes of rumen fermentation parameters, microbial flora, and metabolites in the three groups and their interactions were studied., Results and Discussion: The crude protein level in the supplementary diet had significant effects on pH, valeric acid, and the ratio of acetic acid to propionic acid ( p < 0.05). The protein level had no significant effect on the dominant microflora at the phylum level ( p > 0.05), and all three groups were Bacteroides and Firmicutes. The results of metabolite analysis showed that the crude protein level of supplementary diet significantly affected the metabolic pathways such as Bile secretion and styrene degradation ( p < 0.05), and there were different metabolites between the LP group and HP group, and these different metabolites were related to the dominant microbial to some extent. To sum up, in this experiment, the effects of crude protein level in supplementary diet on rumen microorganisms and metabolites of JY and their relationship were studied, which provided the theoretical basis for formulating a more scientific and reasonable supplementary diet in the future., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Dai, Ma, Dingkao, Huang, La, Li, Ma, Wu, Chu, Guo, Pei, Yan and Liang.)

Spermatogenesis is a complex process that involves proliferation and differentiation of diploid male germ cells into haploid flagellated sperm and requires intricate interactions between testicular somatic cells and germ cells. The cellular heterogeneity of this process presents a challenge in analyzing the different cell types at various developmental stages. Single-cell RNA sequencing (scRNA-seq) provides a useful tool for exploring cellular heterogeneity. In this study, we performed a comprehensive and unbiased single-cell transcriptomic study of spermatogenesis in sexually mature 4-year-old yak using 10× Genomics scRNA-seq. Our scRNA-seq analysis identified six somatic cell types and various germ cells, including spermatogonial stem cells, spermatogonia, early-spermatocytes, late-spermatocytes, and spermatids in yak testis. Pseudo-timing analysis showed that Leydig and myoid cells originated from common progenitor cells in yaks. Moreover, functional enrichment analysis demonstrated that the top expressed genes in yak testicular somatic cells were significantly enriched in the cAMP signaling pathway, PI3K-Akt signaling pathway, MAPK signaling pathway, and ECM receptor interactions. Throughout the spermatogenesis process, genes related to spermatogenesis, cell differentiation, DNA binding, and ATP binding were expressed. Using immunohistochemical techniques, we identified candidate marker genes for spermatogonial stem cells and Sertoli cells. Our research provides new insights into yak spermatogenesis and the development of various types of cells in the testis, and presents more reliable marker proteins for in vitro culture and identification of yak spermatogonial stem cells in the later stage.

The growth hormone receptor ( GHR ) is a member of the cytokine/hematopoietic factor receptor superfamily, which plays an important role in the growth and development, immunity, and metabolism of animals. This study identified a 246 bp deletion variant in the intronic region of the GHR gene, and three genotypes, including type II, type ID, and type DD, were observed. Genotype analysis of structural variation (SV) was performed on 585 individuals from 14 yak breeds, and it was found that 246 bp deletion was present in each breed. The II genotype was dominant in all yak breeds except for SB yak. The association analysis of gene polymorphisms and growth traits in the ASD yak population showed that the 246 bp SV was significantly associated with body length at 6 months ( p < 0.05). GHR messenger RNA (mRNA) was expressed in all the tested tissues, with significantly higher levels in the liver, muscle, and fat than in other organs. The results of transcription activity showed that the luciferase activity of the pGL4.10-DD vector was significantly higher than that of the pGL4.10-II vector ( p < 0.05). Additionally, the transcription-factor binding prediction results showed that the SV in the runt-related transcription factor 1 (Runx1) transcription-factor binding site may affect the transcriptional activity of the GHR gene, regulating yak growth and development. This study showed that the novel SV of the GHR gene could be used as a candidate molecular marker for the selection of the early growth trait in ASD yak.

Testicular development is a tightly regulated process in mammals. Understanding the molecular mechanisms of yak testicular development will benefit the yak breeding industry. However, the roles of different RNAs, such as mRNA, lncRNA, and circRNA in the testicular development of yak, are still largely unclear. In this study, transcriptome analyses were performed on the expression profiles of mRNAs, lncRNAs, and circRNAs in testis tissues of Ashidan yak at different developmental stages, including 6-months-old (M6), 18-months-old (M18), and 30-months-old (M30). A total of 30, 23, and 277 common differentially expressed (DE) mRNAs, lncRNAs, and circRNAs were identified in M6, M18, and M30, respectively. Furthermore, functional enrichment analysis showed that the common DE mRNAs during the entire developmental process were mainly involved in gonadal mesoderm development, cell differentiation, and spermatogenesis processes. Additionally, co-expression network analysis identified the potential lncRNAs related to spermatogenesis, e.g., TCONS&#95;00087394 and TCONS&#95;00012202. Our study provides new information about changes in RNA expression during yak testicular development, which greatly improves our understanding of the molecular mechanisms regulating testicular development in yaks.

Background: Normal testicular development is highly crucial for male reproduction and is a precondition for spermatogenesis that is the production of spermatozoa in the testes. MiRNAs have been implicated in several testicular biological processes, including cell proliferation, spermatogenesis, hormone secretion, metabolism and reproductive regulation. In the present study, we used deep sequencing data to study the functions of miRNAs during testicular development and spermatogenesis, by analyzing the expression patterns of small RNAs in 6-, 18- and 30-month-old yak testis tissues., Results: A total of 737 known and 359 novel miRNAs were obtained from 6-, 18- and 30-month-old yak testes. In all, we obtained 12, 142 and 139 differentially expressed (DE) miRNAs in 30- vs. 18-, 18- vs. 6-, and 30- vs. 6-month-old testes, respectively. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of all DE miRNA target genes revealed BMP2, TGFB2, GDF6, SMAD6, TGFBR2 and other target genes as participants in different biological processes, including TGF-β, GnRH, Wnt, PI3K-Akt, MAPK signaling pathways and several other reproductive pathways. In addition, quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) was used to detect the expression of seven randomly selected miRNAs in 6-, 18- and 30-month-old testes, and the results were consistent with the sequencing data., Conclusions: The differential expression of miRNAs in yak testes at different development stages was characterized and investigated using deep sequencing technology. We believe that the results will contribute to further understanding the functions of miRNAs in regulating the development of yak testes and improving the reproductive performance of male yaks., (© 2023. The Author(s).)

Yaks ( Bos grunniens ) are the only bovine species that adapt well to the harsh high-altitude environment in the Qinghai-Tibetan plateau. However, the reproductive adaptation to the climate of the high elevation remains to be elucidated. Cell composition and molecular characteristics are the foundation of normal ovary function which determines reproductive performance. So, delineating ovarian characteristics at a cellular molecular level is conducive to elucidating the mechanism underlying the reproductive adaption of yaks. Here, the single-cell RNA-sequencing (scRNA-seq) was employed to depict an atlas containing different cell types with specific molecular signatures in the yak ovary. The cell types were identified on the basis of their specifically expressed genes and biological functions. As a result, a cellular atlas of yak ovary was established successfully containing theca cells, stromal cells, endothelial cells, smooth muscle cells, natural killer cells, macrophages, and proliferating cells. A cell-to-cell communication network between the distinct cell types was constructed. The theca cells were clustered into five subtypes based on their biological functions. Further, CYP11A1 was confirmed as a marker gene for the theca cells by immunofluorescence staining. Our work reveals an ovarian atlas at the cellular molecular level and contributes to providing insights into reproductive adaption in yaks.

Yaks ( Bos grunniens ) are a critical livestock breed in the plateau region, and changing the feeding system of yaks can significantly improve their growth performance. The effects of different feeding regimes on the growth performance and meat quality of yaks were comprehensively compared here. The transcriptome and proteome of the Longissimus dorsi muscle were determined using RNA-seq and Tandem Mass Tag (TMT) techniques. Indoor feeding significantly improved the growth performance (such as the average daily gain and carcass weight) and meat quality characteristics compared with traditional grazing feeding. In the grazing (Group G) vs. in-house fed group (Group HF) comparison, 40 differentially expressed genes/differentially abundant proteins exhibited the same mRNA and protein expression trends. These genes were associated with collagen binding, the lipoxygenase pathway, and the arachidonic acid metabolic process. Parallel reaction monitoring verified whether the TMT results were reliable. Moreover, some pathways, such as the AMPK signaling pathway, FoxO signaling pathway, PPAR signaling pathway, and fatty acid metabolism, were significantly enriched. These results expand our knowledge about meat quality in yaks and provide practical information and more evidence for further insight into the biological mechanisms underlying meat quality traits.

Long-haired individuals in the Tianzhu white yak population are a unique genetic resource, and have important landscape value. Copy number variation (CNV) is an important source of phenotypic variation in mammals. In this study, we used resequencing technology to detect the whole genome of 10 long-haired Tianzhu white yaks (LTWY) and 10 normal-haired Tianzhu white yaks (NTWY), and analyzed the differences of CNV in the genome of LTWYs and NTWYs. A total of 110268 CNVs were identified, 2006 CNVRs were defined, and the distribution map of these CNVRs on chromosomes was constructed. The comparison of LTWYs and NTWYs identified 80 differential CNVR-harbored genes, which were enriched in lipid metabolism, cell migration and other functions. Notably, some differential genes were identified as associated with hair growth and hair-follicle development (e.g., ASTN2 , ATM , COL22A1 , GK5 , SLIT3 , PM20D1 , and SGCZ ). In general, we present the first genome-wide analysis of CNV in LTWYs and NTWYs. Our results can provide new insights into the phenotypic variation of different hair lengths in Tianzhu white yaks.

The interaction between the dermal papilla cells (DPCs) and epidermal hair matrix cells (HMCs) of hair follicles (HFs) is crucial for the growth and development of HFs, but the molecular mechanism is complex and remains unclear. MicroRNAs (miRNAs) are the key signaling molecules for cellular communication. In this study, the DPCs and HMCs of yak were isolated and cultured, and the differentially expressed mRNA and miRNA were characterized to analyze the molecular basis of the interaction between DPCs and HMCs during hair follicle (HF) development in yak. The mRNA differential expression and functional enrichment analysis revealed that there were significant differences between DPCs and HMCs, and they showed the molecular functional characteristics of dermal cells and epidermal cells, respectively. Multiple KEGG pathways related to HF development were enriched in the highly expressed genes in DPCs, while the pathways associated with microbiota and immunity were significantly enriched in the highly expressed genes in HMCs. By combining analysis with our previous 10× genomics single-cell transcriptome data, 39 marker genes of DPCs of yak were identified. A total of 123 relatively specifically expressed miRNAs were screened; among these, the miRNAs associated with HF development such as miR-143, miR-214, miR-125b, miR-31, and miR-200 were presented. In conclusion, the large changes in yak DPCs and HMCs for both mRNA and miRNA expression were revealed, and numerous specifically expressed mRNAs and miRNAs in DPCs or HMCs were identified, which may contribute to the interaction and cellular communication between DPCs and HMCs during HF development in yak.

Background: As the direct organ of villus, hair follicles have obvious seasonal cycles. The hair follicle cycle is orchestrated by multiple cell types that together direct cell renewal and differentiation. But the regulation property of hair follicle cells from anagen to catagen in yak is still unknown., Results: In this study, single-cell RNA sequencing was performed on 24,124 single cells of the scapular skin from white yak. Based on tSNE cluster analysis, the cell types of IFE-DC, epidermal cell lines, fibroblasts, keratinocytes, IRS, DS, INFU, and other cells in yak hair follicles during anagen and catagen were successfully identified, and the gene expression profiles were described. The GO enrichment analysis indicated the different cells characteristic genes to be mainly enriched in the epidermal development, epithelial cell differentiation and wound healing pathways. The pseudotime trajectory analysis described the differentiation trajectory of the epidermal lineage and dermal lineage of the hair follicle during anagen and catagen. Moreover, the dynamic changes of the genes like LHX2, KRT25, and KRT71 were found to be highly expressed in HS and IRS, but not in the IFE-DC, INFU, and keratinocyte during differentiation., Conclusions: Our results analyzed the time-varying process of gene expression in the dermal cell lineage and epidermal cell lineage of hair follicles during anagen and catagen during fate differentiation was expounded at the single cell level, revealing the law of fate specialization of different types of cells. In addition, based on the enrichment analysis, the transcriptional regulatory factors involved in the different cell fates were also revealed. These results will help to enhance our understanding of yak hair follicle cycle and promote the development and utilization of yak villus., (© 2022. The Author(s).)