Genetic diversity, phylogeography, and maternal origin of yak (Bos grunniens).

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.

Age-dependent changes in the expression and localization of LYZL4, LYZL6 and PCNA during testicular development in the Ashidan yak.

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.

Identification and profiling of microRNAs during yak's testicular development.

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.

Transcriptomic analysis reveals that non-forage or forage fiber source promotes rumen development through different metabolic processes in lambs.

Dietary fiber supplementation can stimulate rumen development in lambs during the pre-weaning period. However, it is unclear whether different sources of fiber have varying effects on rumen development. This study aimed to investigate the molecular mechanism of rumen morphological and functional development based on non-forage or forage as a starter dietary fiber source. Twenty-four male Hu lambs with similar body weights (BW, 3.67 ± 0.08 kg) were selected and divided into two groups that received diets supplemented with either alfalfa hay (AH) or soybean hull (SH). At the age of 70 days, six lambs were slaughtered from each treatment group for rumen fermentation and morphological analyses. Three samples of the rumen tissue from the ventral sac were collected for transcriptomic analysis. The results identified 633 differentially expressed genes (DEGs), of which 210 were upregulated and 423 were downregulated in the SH group compared with those in the AH group. The upregulated DEGs were most enriched in the immune function and proteolysis pathways, whereas the downregulated DEGs were mainly involved in cell proliferation, apoptosis, and differentiation pathways. These findings indicated that non-forage as a starter dietary fiber source improved immune function and enhanced nitrogen utilization, whereas forage facilitated rumen morphological development.

Polymorphisms in BMPRIB gene affect litter size in Chinese indigenous sheep breed.

The BMPRIB gene belongs to the TGF-ß superfamily and is considered to be a regulator of sheep reproductive performance. Single nucleotide polymorphisms (SNPs) of BMPRIB gene in the Small Tail Han, Hu, Mongolian, Oula, Gansu Alpine Fine-wool, Dorper and Australian White sheep were detected by Sanger sequencing. Five SNPs (rs427897187 G > A, rs418841713 A > G, rs159952533 T > C, rs429416173 C > A and rs403555643 A > G) of BMPRIB gene were identified. For rs427897187 G > A, further analysis revealed that genotype GG and GA had 0.26 (p < 0.05) and 0.33 (p < 0.05) litter size less than those with genotype AA in Oula sheep. For rs403555643 A > G, further analysis revealed that genotype GG and AG had 0.65 (p < 0.05) and 0.38 (p < 0.05) litter size more than those with genotype AA in Oula sheep, and genotype GG had 0.56 (p < 0.05) litter size more than those with genotype AA in Mongolian sheep. The results showed that rs427897187 G > A and rs403555643 A > G are potential molecular markers wich could improve litter size of Chinese indigenous sheep and be used in Chinese indigenous sheep breeding.

Effects of starter feeding and early weaning on developmental expressions of IGF-I gene in liver and IGF-IR gene in rumen of lambs from birth to eighty-four days of age.

The early weaning and starter feeding have significant effects on lamb growth and digestive tract development. However, it is not clear whether the expression of IGF-1/IGF-1R genes could be affected by feeding starter and weaning. Therefore, a total of 102 Hu male neonatal lambs were randomly divided into 3 groups: fed with starter at 42d + weaned at 56d (group A), fed with starter at 7d + weaned at 28d (group B) and fed with starter at 7d + weaned at 56d (group C), to explore the effects of starter feeding and weaning age on developmental expressions of IGF-I gene in liver and IGF-IR gene in rumen of Hu sheep. The results showed that IGF-I and IGF-IR genes were expressed extensively in various tissues of lambs, the expression of IGF-I was significantly higher in liver (p < 0.01), while the expression of IGF-IR was higher in rumen among gastrointestinal tissues. The mRNA level of IGF-I of group C was strikingly higher than that of group A at 28, 70, and 84d (p < 0.01) in liver, respectively. Rumen mRNA level of IGF-IR of group C was prominently higher than that of group A at 70d and 84d (p < 0.05), respectively. However, IGF-IR expression of group A was higher than that of group C at 14 and 42d (p < 0.05). The mRNA level of IGF-I of group B was significantly higher than that of group C at 42, 56, 70, and 84d (p < 0.01) in liver, while IGF-IR expression of group B was higher than that of group C at 70d (p < 0.01). In conclusion, early weaning and starter feeding affected liver IGF-I mRNA and rumen IGF-IR mRNA expression.

Genome-Wide Landscape of mRNAs, lncRNAs, and circRNAs during Testicular Development of 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_00087394 and TCONS_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.

Single-Cell Transcriptomics Analysis Reveals a Cell Atlas and Cell Communication in Yak Ovary.

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.

Single-Cell RNA Sequencing Reveals Atlas of Yak Testis Cells.

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 association of polymorphisms in the ovine PPARGC1B and ZEB2 genes with body weight in Hu sheep.

The aims of this study were to analyze the effects of PPARGC1B and ZEB2 polymorphisms on the body weight of Hu sheep. DNA sequencing and KASPar technologies were used to detect single nucleotide polymorphisms (SNPs) within the PPARGC1B and ZEB2 genes of Hu sheep (n = 207). Two synonymous mutations, g.300 G > A and g.645 C > T, were detected in PPARGC1B and ZEB2, respectively. The body weights of sheep were recorded at 80, 100, 120, 140, 160 and 180-days, and associations between these polymorphisms and body weight changes were analyzed. Association analysis demonstrated that the polymorphisms in PPARGC1B and ZEB2 significantly associated with body weight (p < 0.05). At the g.300 G > A locus, individuals with the GG genotype had significantly higher body weight than those with the AA genotype, and at the g.645 C > T locus, individuals with the TT genotype had significantly higher body weight than those with the TC genotype. Individuals with both polymorphisms exhibited significantly different growth (p < 0.05). These data suggest that polymorphisms in the PPARGC1B and ZEB2 genes can be used as candidate molecular markers for the breeding of desirable growth traits in Hu sheep.

Different expression of LHR, PRLR, GH and IGF1 during testicular development of yak.

Luteinizing hormone receptor (LHR), prolactin receptor (PRLR), growth hormone (GH) and insulin-like growth factor 1 (IGF1) have been shown to be key regulators of germ cell development. However, the role of LHR, PRLR, GH and IGF1 in the development of yak testis remains unclear. In this study, we aimed to describe and compare gene expression and protein localization of LHR, PRLR, GH and IGF1 in the development of yak testes. Testes were collected from 6, 24, 36 and 72 months yak, and the kidney, liver, testicular, lung, skeletal muscle, heart and spleen tissues were collected from 36 months yak. The quantitative real-time PCR (qRT-PCR) results showed that the expression of these four genes was widely expressed in kidney, liver, testicular, lung, skeletal muscle, heart and spleen, while the LHR and PRLR were highly expressed in the kidney, skeletal muscle and testis, and higher levels of GH and IGF were expressed in spleen and testis. Moreover, the mRNA expression of these genes in adults was higher than in pre-pubertal yak. In the testis, the LHR-, PRLR-, GH- and IGF1-positive signals were detected in the Leydig cells of the 6 months, while the intense positive signals were discovered in Leydig cells, spermatogonia and spermatocytes of the 36 and 72 months. Thus, LHR, PRLR, GH and IGF1 may be involved in the development of spermatids and spermatocytes, and in the regulation of spermatogonia proliferation and Leydig cell function.

Identification and Characterization of Piwi-Interacting RNAs for Early Testicular Development in Yak.

Normal testicular development plays a crucial role in male reproduction and is the precondition for spermatogenesis. PIWI-interacting RNAs (piRNAs) are novel noncoding RNAs expressed in animal germ cells that form complexes with PIWI family proteins and are involved in germ cell development, differentiation, and spermatogenesis. However, changes in piRNA expression profiles during early testicular development in yak have not been investigated. In this study, we used small RNA sequencing to evaluate the differences and potential functions of piRNA expression profiles in 6-, 18-, and 30-month-old yak testis tissues. Differential expression analysis found 109, 293, and 336 differentially expressed piRNAs in M30 vs. M18, M18 vs. M6, and M30 vs. M6, respectively, and found 30 common differentially expressed piRNAs in the three groups of M6, M18, and M30. In addition, the functional enrichment analysis of differentially expressed piRNAs target genes indicated that they were related to testicular development and spermatogenesis. Finally, we detected the expression of the PIWI protein family in the yak testis at different developmental stages and found that PIWIL1, PIWIL2, PIWIL3, and PIWIL4 were highly expressed in 18- and 30-month-old yak testis and almost not expressed in 6-month-old yak testis. In conclusion, this study summarizes the changes of piRNA expression patterns during the early development of yak testis and provides new clues for the regulatory role of piRNA in yak testis.

Polymorphisms of the Ovine BMPR-IB, BMP-15 and FSHR and Their Associations with Litter Size in Two Chinese Indigenous Sheep Breeds.

The Small Tailed Han sheep and Hu sheep are two prolific local sheep in China. In this study, the polymorphisms of BMPR-IB (Bone morphogenetic protein receptor IB), BMP-15 (Bone morphogenetic protein 15) and FSHR (follicle stimulating hormone receptor) were investigated to check whether they are associated with litter size in Small Tailed Han sheep and Hu sheep. Consequently, three polymorphisms, FecB mutation in BMPR-IB (c.746A>G), FecG mutation in BMP-15 (c.718C>T) and the mutation (g. 47C>T) in FSHR were found in the above two sheep breeds with a total number of 1630 individuals. The single marker association analysis showed that the three mutations were significantly associated with litter size. The ewes with genotype FecBB/FecBB and FecBB/FecB+ had 0.78 and 0.58 more lambs (p < 0.01) than those with genotype FecB+/FecB+, respectively. The heterozygous Han and Hu ewes with FecXG/FecX+ genotype showed 0.30 (p = 0.05) more lambs than those with the FecX+/FecX+ genotype. For FSHR gene, the ewes with genotype CC had 0.52 (p < 0.01) and 0.75 (p < 0.01) more lambs than those with genotypes TC and TT, respectively. Combined effect analyses indicated an extremely significant interaction (p < 0.01) between the random combinations of BMPR-IB, BMP-15 and FSHR genes on litter size. In addition, the Han and Hu ewes with BB/G+/CC genotype harbor the highest litter size among ewes analyzed in current study. In conclusion, BMPR-IB, BMP-15 and FSHR polymorphisms could be used as genetic markers in multi-gene pyramiding for improving litter size in sheep husbandry.

Polymorphisms within the PRKG1 Gene of Gannan Yaks and Their Association with Milk Quality Characteristics.

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.

Polymorphisms of ITGA9 Gene and Their Correlation with Milk Quality Traits in Yak (Bos grunniens).

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.

Serum Metabolomic Analysis of Synchronous Estrus in Yaks Based on UPLC-Q-TOF MS Technology.

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.

Proteomic Analysis Reveals the Effects of Different Dietary Protein Levels on Growth and Development of Jersey-Yak.

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.

Genome-Wide Transcriptome Profiling Reveals the Mechanisms Underlying Hepatic Metabolism under Different Raising Systems in Yak.

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.

Association between Single Nucleotide Polymorphisms of PRKD1 and KCNQ3 Gene and Milk Quality Traits in Gannan Yak (Bos grunniens).

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 cGPS 7K liquid-phase microarray technology. Milk composition was analyzed using a MilkoScanTM 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.

Whole-transcriptome analysis of longissimus dorsi muscle in cattle-yaks reveals the regulatory functions of ADAMTS6 gene in myoblasts.

Cattle-yak, which is the hybrid F1 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.