The Tianzhu white yaks (Bos grunniens) live on the Qinghai–Tibet Plateau. During winter, a lack of resources and low nutritional levels seriously affect their growth performance. In this study, we aimed to explore the effect of supplementation straw pellets on the growth performance and ruminal microbiota of yaks. Overall, at 6 (6M, n = 24), 18 (18M, n = 26), 30 (30M, n = 20), 42 (42M, n = 24), and 54 (54M, n = 22) month old Tianzhu white yaks were selected (total n = 116) and divided into the mixed straw + grazing (MSG), corn straw + grazing (CSG), and the grazing control (G) groups according to age and gender. Their growth performance was measured as per different dietary treatments. The rumen microbial community structure and levels of VFAs were analyzed from the 6M, 30M, and 54M male yaks from each group. The supplementary diets led to an increase in the ADG, which was the highest in the MSG group. The MSG group exhibited the highest level of acetate and total VFAs (TVFAs) among the three groups (p < 0.05). In addition, the 16S rRNA sequencing results proved that the microbial composition was dominated by the members of Firmicutes and Bacteroidetes. Christensenellaceae R-7 group was significantly abundant in the CSG and MSG groups compared to the G group (p < 0.05). Principal coordinate analysis (PCoA) revealed that the bacterial community structure of rumen in the MSG and CSG groups was considerably different from that in the G group; 6M samples exhibited different rumen microbial diversity compared with the other samples. Correlation analysis revealed that Christensenellaceae_R-7_group was positively correlated with the levels of acetate, TVFAs, and ADG. These results demonstrated that mixed straw pellets improved the growth performance of yaks, increased the abundance of Christensenellaceae R-7_group involved in cellulose degradation in the rumen, and produced large amounts of VFAs, which were absorbed by yaks, thus increasing their ADG. This study provides new insights into the effects of straw pellet supplementation on the changes in the rumen microbiota and growth performance of yaks.

Acetyl-CoA carboxylase beta (ACACB) is a functional candidate gene that impacts fat deposition. In the present study, we sequenced exon 37–intron 37, exon 46–intron 46, and intron 47 of yak ACACB using hybrid pool sequencing to search for variants and genotyped the gene in 593 Gannan yaks via Kompetitive allele-specific polymerase chain (KASP) reaction to determine the effect of ACACB variants on carcass and meat quality traits. Seven single nucleotide polymorphisms were detected in three regions. Eight effective haplotypes and ten diplotypes were constructed. Among them, a missense variation g.50421 A > G was identified in exon 37 of ACACB, resulting in an amino acid shift from serine to glycine. Correlation analysis revealed that this variation was associated with the cooking loss rate and yak carcass weight (p = 0.024 and 0.012, respectively). The presence of haplotypes H5 and H6 decreased Warner–Bratzler shear force (p = 0.049 and 0.006, respectively), whereas that of haplotypes H3 and H4 increased cooking loss rate and eye muscle area (p = 0.004 and 0.034, respectively). Moreover, the presence of haplotype H8 decreased the drip loss rate (p = 0.019). The presence of one and two copies of haplotypes H1 and H8 decreased the drip loss rate (p = 0.028 and 0.004, respectively). However, haplotype H1 did not decrease hot carcass weight (p = 0.011), whereas H3 increased the cooking loss rate (p = 0.007). The presence of one and two copies of haplotype H6 decreased Warner–Bratzler shear force (p = 0.014). The findings of the present study suggest that genetic variations in ACACB can be a preferable biomarker for improving yak meat quality.

As an indigenous species on the Tibetan Plateau, the yak is well adapted to the plateau hypoxic environment. The high-altitude hypoxia adaptation of the yak requires the adaptive reshaping of multiple tissues and organs, especially the lungs. To reveal the adaptive development of yak lungs under hypoxic stress at the tissue and molecular levels, we conducted histomorphological observations as well as transcriptomic and metabolomic studies of yak lungs at three ages (0.5, 2.5, and 4.5 years). The results showed that the lung tissue developed significantly with age. The mean alveolar area was higher (p < 0.01) in 4.5 and 2.5-year-old yaks than in 0.5-year-old yaks. The percentage of elastic fibers, micro-arterial wall thickness, and micro-arterial area showed an increasing trend (p < 0.01) from 0.5-year-old yaks to 2.5-year-old yaks and then to 4.5-year-old yaks. In addition, some critical differentially expressed genes related to angiogenesis (MYC, EPHA2, TNF), fiber formation (EREG), smooth muscle proliferation (HBEGF), erythropoiesis (SOCS3), and hypoxia response (ZFP36) were identified. Some metabolites associated with these genes were also found simultaneously. These findings provide a deeper understanding of the molecular strategies underlying this species’ extraordinary ability to survive normally in low-oxygen environments. In conclusion, the lungs of yaks undergo continuous adaptive development under hypoxic stress, and these findings are crucial for understanding the molecular mechanisms by which native species of the Tibetan Plateau survive in harsh environments.

The FOXO3 gene, a prominent member of the FOXO family, has been identified as a potential quantitative trait locus for muscle atrophy and lipid metabolism in livestock. It is also considered a promising candidate gene for meat quality traits such as Warner–Bratzler shear force (WBSF) and water holding capacity (WHC). The aim of this study was to identify sequence mutations in the FOXO3 gene of yaks and to analyze the association of genotypes and haplotypes with meat traits such as WBSF and WHC. Quantitative reverse-transcriptase PCR (RT-qPCR) was applied to determine the expression levels of FOXO3 in yak tissues, with the results revealing a high expression in the yak longissimus dorsi muscle. Exons of the FOXO3 gene were then sequenced in 572 yaks using hybrid pool sequencing. Five single nucleotide polymorphisms were identified. Additionally, four effective haplotypes and four combined haplotypes were constructed. Two mutations of the FOXO3 gene, namely C>G at exon g.636 and A>G at exon g.1296, were associated with cooked meat percentage (CMP) (p < 0.05) and WBSF (p < 0.05), respectively. Furthermore, the WBSF of the H2H3 haplotype combination was significantly lower than that of other combinations (p < 0.05). The findings of this study suggest that genetic variations in FOXO3 could be a promising biomarker for improving yak meat traits.

Hormone-sensitive lipase (HSL) is involved in the breakdown of triacylglycerols in adipose tissue, which influences muscle tenderness and juiciness by affecting the intramuscular fat content (IMF). This study analyzed the association between different genotypes and haplotypes of the yak HSL gene and carcass and meat quality traits. We used hybridization pool sequencing to detect exon 2, exon 8, and intron 3 variants of the yak HSL gene and genotyped 525 Gannan yaks via KASP to analyze the effects of the HSL gene variants on the carcass and meat quality traits in yaks. According to the results, the HSL gene is highly expressed in yak adipose tissue. Three single nucleotide polymorphisms (SNPs) were identified, with 2 of them located in the coding region and one in the intron region. Variants in the 2 coding regions resulted in amino acid changes. The population had 3 genotypes of GG, AG, and AA, and individuals with the AA genotype had lower WBSF values (p < 0.05). The H3H3 haplotype combinations could improve meat tenderness by reducing the WBSF values and the cooking loss rate (CLR) (p < 0.05). H1H1 haplotype combinations were associated with the increased drip loss rate (DLR) (p < 0.05). The presence of the H1 haplotype was associated the increased CLR in yaks, while that of the H2 haplotype was associated with the decreased DLR in yaks (p < 0.05). These results demonstrated that the HSL gene may influence the meat quality traits in yaks by affecting the IMF content in muscle tissues. Consequently, the HSL gene can possibly be used as a biomarker for improving the meat quality traits in yaks in the future.

As a class of regulatory factors, microRNAs (miRNAs) play an important role in regulating normal muscle development and fat deposition. Muscle and adipose tissues, as major components of the animal organism, are also economically important traits in livestock production. However, the effect of miRNA expression profiles on the development of muscle and adipose tissues in yak is currently unknown. In this study, we performed RNA sequencing (RNA-Seq) on Tianzhu white yak longissimus dorsi muscle tissue obtained from calves (6 months of age, M6, n = 6) and young (30 months of age, M30, n = 6) and adult yak (54 months of age, M54, n = 6) to identify which miRNAs are differentially expressed and to investigate their temporal expression profiles, establishing a regulatory network of miRNAs associated with the development of muscle and adipose. The results showed that 1191 miRNAs and 22061 mRNAs were screened across the three stages, of which the numbers of differentially expressed miRNAs (DE miRNAs) and differentially expressed mRNAs (DE mRNAs) were 225 and 450, respectively. The expression levels of the nine DE miRNAs were confirmed using a reverse transcription quantitative PCR (RT-qPCR) assay, and the trend of the assay results was generally consistent with the trend of the transcriptome profiles. Based on the expression trend, DE miRNAs were categorized into eight different expression patterns. Regarding the expression of DE miRNAs in sub-trends Profile 1 and Profile 2 (p < 0.05), the gene expression patterns were upregulated (87 DE miRNAs). Gene ontology (GO) and Kyoto Encyclopedia of Genes Genomes (KEGG) analyses showed that the identified DE miRNAs and DE mRNAs were enriched in pathway entries associated with muscle and intramuscular fat (IMF) growth and development. On this basis, we constructed a DE miRNA–mRNA interaction network. We found that some DE mRNAs of interest overlapped with miRNA target genes, such as ACSL3, FOXO3, FBXO30, FGFBP4, TSKU, MYH10 (muscle development), ACOX1, FADS2, EIF4E2, SCD1, EL0VL5, and ACACB (intramuscular fat deposition). These results provide a valuable resource for further studies on the molecular mechanisms of muscle tissue development in yak and also lay a foundation for investigating the interactions between genes and miRNAs.