Search

MicroRNAs (miRNAs) were identified to be involved in various biological functions by regulating the degradation or suppressing the translation of their downstream target genes. Recent studies have identified miR-29a play a key role in functions of mammal cell differentiation, proliferation, apoptosis, and signal transduction. However, the underlying functions for miR-29a in jejunal epithelial cells biological function still to be investigated. In order to explore the yak jejunal epithelial cells proliferation and barrier dysfunction with over expression of miR-29a gene, three 0-day-old Pamir male yaks were randomly selected and slaughtered in present study, and the jejunal epithelial cells were isolated and cultured to determine yak jejunal epithelial cells proliferation and protein composition on differential expression of miR-29a gene in Pamir plateau. Here, we demonstrated that the overexpression of miR-29a gene could inhibit the proliferation of Pamir yaks jejunum epithelial cells, and contribute to the apoptosis of Pamir yaks jejunal epithelial cells with some extent. A total of 133 differentially expressed proteins were identified in different expression of miR-29a groups by label-free Mass Spectrometry (MS), which could be concluded to two predominant themes: cell proliferation and inflammatory response. Interestingly, GPR41, as a bridge protein, was contacted two predominant themes to involved in Pamir Yaks jejunal mechanical barrier PPI network, and the target proteins displayed strong mutual interactions in the complex PPI network. Overall, our study suggested that the over-expression miR-29a inhibited the jejunal epithelial cells proliferation and the expressions of specific proteins, which damaged jejunal barrier function to slow down the intestine structure and function advanced mature development during young livestock period for influence the enhanced performance of production efficiency.

Milk lipids greatly affect the volatile flavor of milk, and the relationship between lipids and volatile flavor in yak milk was explored in this study. The volatile flavor compounds (VFCs), lipids profile, fatty acids in yak ordinary milk and colostrum were detected with HP/SPME-GC–MS, the semiquantitative lipidomics based on LC-MS/MS, GC–MS, respectively. The VFCs differences in yak milk were closely related to 1-((1 s,3ar,4r,7 s,7as)-4-hydroxy-7-isopropyl-4-methyloctahloctahydro-1h-inden-1-Yl)-ethanone,2,6,6-trimethyl-2,4-cycloheptadien-1-one, pentanal, 2-phenylethyl propionate, octanoic acid methyl ester, diphosphoric acid diisooctyl ester, (Z)-3,4,4-trimethyl-5-oxo-2-hexenoic acid and acetic acid. The volatile flavor in yak milk was well correlated with milk lipids, and TG(4:0_12:3_18:1), TG(6:0_8:0_18:1), TG(4:0_12:3_18:1), TG(12:0_18:2_18:3) and TG(16:0e_18:1_22:5) were the crucial lipid molecules affecting volatile flavor. The degeneration of above lipids by hydrolysis produced some fatty acids and alcohol, then these compounds were further derived into other VFCs especially above crucial 8 molecules. This study provided a theoretical basis for improving the volatile flavor by controlling lipids in yak milk.

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.

Abstract 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

Abstract 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.

Yaks (Bos grunniens) exhibit exceptional adaptation to the challenging high-altitude environment of the Qinghai-Tibetan plateau, making them the sole bovine species capable of thriving in such exreme conditions. Investigating the cellular and molecular characteristics of yak ovaries across different reproductive states is crucial for gaining insight into their ovarian functions. Herein, the cellular atlases of yak ovaries in different reproductive states were depicted by single-cell RNA-sequencing (scRNA-seq). The cellular atlases of the ovaries were established by identifying specific gene expression patterns of various cell types, including granulosa cells, theca cells, stromal cells, smooth muscle cells, endothelial cells, glial cell, macrophages, natural killer cells, and proliferating cells. The cellular compositions of the ovaries vary among different reproductive states. Furthermore, the granulosa cells comprise six cell subtypes, while theca cells consist of eight cell subtypes. The granulosa cells and theca cells exhibit distinct biological functions throughout different reproductive states. The two cell types were aligned along their respective pseudotime trajectories. Moreover, a cell-to-cell communication network was constructed among distinct cell types within the ovary, spanning the three reproductive states. Notably, during the estrus period, the granulosa cells demonstrated more prominent interactions with other cell types compared to the remaining reproductive states.

Abstract 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_0617 converging on bta-miR-2285q, mirrored by KLF15/circRNA_0345/bta-miR-20b and CTSF/circRNA_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.

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.

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.

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.

Abstract 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.

Abstract 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.

Abstract Background The hair coat is available for the yak to live in the harsh environment of the plateau. Besides, improving the hair production of yak is necessary for its textile industry development. Hair grows from hair follicles (HFs). The HFs undergo periodic growth after birth and are regulated by the complex gene regulatory network. However, the molecular mechanism of HFs regeneration in the Tianzhu white yak remains unclear. RNA editing is a post-transcriptional mechanism that regulates gene expression and produces new transcripts. Hence, we investigated the influence of the A-to-I RNA editing events on the HFs cycle of the Tianzhu white yak. Results We finally identified 54,707 adenosine-to-inosine (A-to-I) RNA editing sites (RESs) from RNA sequencing data of the HFs cycle in the Tianzhu white yak. Annotation results showed RESs caused missense amino acid changes in 7 known genes. And 202 A-to-I editing sites altered 23 target genes of 140 microRNAs. A total of 1,722 differential RESs were identified during the HFs cycle of Tianzhu white yak. GO and KEGG enrichment analysis revealed several signaling pathways and GO terms involved skin development, hair growth, and HFs cycle. Such as genes with differential RNA editing levels were significantly enriched in the peroxisome, metabolic pathways, Notch signaling pathway, and PPAR signaling pathway. Besides, the editing sites in HFs development-related genes FAS, APCDD1, WWOX, MPZL3, RUNX1, KANK2, DCN, DSC2, LEPR, HEPHL1, and PTK2B were suggested as the potential RESs involving HFs development. Conclusion This study investigated the global A-to-I RNA editing events during the HFs cycle of yak skin tissue and expanded the knowledge of A-to-I RNA editing on the HFs cycle. Furthermore, this study revealed that RNA editing-influenced genes may regulate the HFs cycle by participating in the HFs development-related pathways. The findings might provide new insight into the regulation of RNA editing in hair growth.

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.

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.

Abstract Background High spontaneous miscarriage rate in yak, especially during late pregnancy, have caused a great economic loss to herdsmen living in the Qinghai-Tibet plateau. However, the mechanism underlying spontaneous miscarriage is still poorly understood. In the present study, placenta protein markers were identified to elucidate the pathological reasons for yak spontaneous miscarriage through isobaric tags for relative and absolute quantification (iTRAQ) proteomic technology and bioinformatic approaches. Results Subsequently, a total of 415 differentially expressed proteins (DEPs) were identified between aborted and normal placentas. The up-regulated DEPs in the aborted placentas were significantly associated with “spinocerebellar ataxia”, “sphingolipid signalling”, “relaxin signalling”, “protein export”, “protein digestion and absorption” and “aldosterone synthesis and secretion” pathway. While the down-regulated DEPs in the aborted placentas mainly participated in “valine, leucine and isoleucine degradation”, “PPAR signalling”, “peroxisome”, “oxidative phosphorylation”, “galactose metabolism”, “fatty acid degradation”, “cysteine and methionine metabolism” and “citrate cycle” pathway. Conclusions The results implied that the identified DEPs could be considered as placental protein markers for yak miscarriage during late pregnancy, and biomacromolecule metabolic abnormality and oxidative damage might be responsible for the high spontaneous miscarriage rate in yak. These findings provide an important theoretical basis for deciphering the pathologic mechanism of late spontaneous miscarriage in yak.

Abstract Background and aim Yak estrus is a seasonal phenomenon, probably involving epigenetic regulation of synthesis and secretion of sex hormones as well as growth and development of follicles. N6-methyladenosine (m6A) is the most common internal modification of the eukaryotic mRNA. However, there are no detailed reports on the m6A transcriptome map of yak ovary. Therefore, this study aimed to collected the yak ovarian tissues at three different states of anestrus (YO-A), estrus (YO-F), and pregnancy (YO-P), and obtained the full transcriptome m6A map in yak by MeRIP-seq. Results The HE staining revealed that the number of growing follicles and mature follicles in the ovary during the estrus period was relatively higher than those in the anestrus period and the pregnancy period. The RT-qPCR showed that the expression of METTL3, METTL14, FTO, YTHDC1 were significantly different across different periods in the ovaries, which suggests that m6A may play a regulatory role in ovarian activity. Next, we identified 20,174, 19,747 and 13,523 m6A peaks in the three ovarian samples of YO-A, YO-F and YO-P using the methylated RNA immunoprecipitation sequencing (MeRIP-seq). The m6A peaks are highly enriched in the coding sequence (CDS) region and 3′untranslated region (3′UTR) as well as the conserved sequence of “RRACH.” The GO, KEGG and GSEA analysis revealed the involvement of m6A in many physiological activities of the yak’s ovary during reproductive cycle. The association analysis found that some genes such as BNC1, HOMER1, BMP15, BMP6, GPX3, and WNT11 were related to ovarian functions. Conclusions The comparison of the distribution patterns of methylation peaks in the ovarian tissues across different periods further explored the m6A markers related to the regulation of ovarian ovulation and follicular development in the yak ovary. This comprehensive map provides a solid foundation for revealing the potential function of the mRNA m6A modification in the yak ovary.

Studying the mechanism of spermatogenesis is key to exploring the reproductive characteristics of male yaks. Although N6-methyladenosine (m6A) RNA modification has been reported to regulate spermatogenesis and reproductive function in mammals, the molecular mechanism of m6A in yak testis development and spermatogenesis remains largely unknown. Therefore, we collected testicular tissue from juvenile and adult yaks and found that the m6A level significantly increased after sexual maturity in yaks. In MeRIP-seq, 1702 hypermethylated peaks and 724 hypomethylated peaks were identified. The hypermethylated differentially methylated RNAs (DMRs) (CIB2, AK1, FOXJ2, PKDREJ, SLC9A3, and TOPAZ1) mainly regulated spermatogenesis. Functional enrichment analysis showed that DMRs were significantly enriched in the adherens junction, gap junction, and Wnt, PI3K, and mTOR signaling pathways, regulating cell development, spermatogenesis, and testicular endocrine function. The functional analysis of differentially expressed genes showed that they were involved in the biological processes of mitosis, meiosis, and flagellated sperm motility during the sexual maturity of yak testis. We also screened the key regulatory factors of testis development and spermatogenesis by combined analysis, which included BRCA1, CREBBP, STAT3, and SMAD4. This study indexed the m6A characteristics of yak testicles at different developmental stages, providing basic data for further research of m6A modification regulating yak testicular development.

Abstract Background Ovarian follicle fluid (FF) as a microenvironment surrounding oocyte plays critical roles in physio-biochemical processes of follicle development and oocyte maturation. It is hypothesized that proteins in yak FF participate in the physio-biochemical pathways. The primary aims of this study were to find differentially expressed proteins (DEPs) between mature and immature FF, and to elucidating functions of the mature and immature FF in yak. Results The mature and immature FF samples were obtained from three healthy yaks that were nonpregnant, aged from four to five years, and free from any anatomical reproductive disorders. The FF samples were subjected to mass spectrometry with the isobaric tags for relative and absolute quantification (iTRAQ). The FF samples went through correlation analysis, principle component analysis, and expression pattern analysis based on quantification of the identified proteins. Four hundred sixty-three DEPs between mature and immature FF were identified. The DEPs between the mature and immature FF samples underwent gene ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG), and protein-protein interaction (PPI) analysis. The DEPs highly expressed in the mature FF mainly took parts in the complement and coagulation cascades, defense response, acute-phase response, response to other organism pathways to avoid invasion of exogenous microorganisms. The complement activation pathway contains eight DEPs, namely C2, C5, C6, C7, C9, C4BPA, CFH, and MBL2. The three DEPs, CATHL4, CHGA, and PGLYRP1, take parts in defense response pathway to prevent invasion of exogenetic microorganism. The coagulation cascades pathway involves many coagulation factors, such as F7, F13A1, FGA, FGB, FGG, KLKB1, KNG1, MASP1, SERPINA1, and SERPIND1. While the DEPs highly expressed in the immature FF participated in protein translation, peptide biosynthetic process, DNA conformation change, and DNA geometric change pathways to facilitate follicle development. The translation pathway contains many ribosomal proteins, such as RPL3, RPL5, RPS3, RPS6, and other translation factors, such as EIF3J, EIF4G2, ETF1, MOV10, and NARS. The DNA conformation change and DNA geometric change involve nine DEPs, DDX1, G3BP1, HMGB1, HMGB2, HMGB3, MCM3, MCM5, MCM6, and RUVBL2. Furthermore, the expressed levels of the main DEPs, C2 and SERPIND1, were confirmed by western blot. Conclusions The differential proteomics revealed the up-regulated DEPs in mature FF take parts in immunoreaction to prevent invasion of microorganisms and the up-regulated DEPs in immature FF participate in protein synthesis, which may improve our knowledge of the follicular microenvironment and its biological roles for reproductive processes in yak. The DEPs, C2 and SERPIND1, can be considered as protein markers for mature yak follicle.

Background: Qaidam cattle are local breeds that habitats in northwest China. It has many excellent characteristics, such as high cold and roughage tolerance, low oxygen adaptability, and tender meat quality. Copy number variation (CNV) can induce phenotypic changes in animals by a variety of effects, and thus affects the biological functions of the animals. To explore the molecular mechanism of its adaptation to extreme cold weather and muscle fat development, the CNV variations in the genome of three Qaidam cattle were detected by whole-genome sequencing, in this study. Results: A total of 16,743 CNVs and 9498 copy number variable regions (CNVRs) were obtained after the screening, which accounts for 2.18% of the bovine genome. The CNVR length detected ranged from 0.3 KB to 10.77 KB, with a total length of 58.17 MB and an average length of 6.12 KB/ CNVR. Through functional enrichment of CNVR related genes, LDHB, and ME1 genes were screened as the key genes for Qaidam cattle to adapt to the cold and low oxygen environments, whereas KIT and FGF18 genes might be related to the coat color and growth. In the CNVR overlapped with QTLs, variation in CAPN1 and CAST genes might be closely related to the tender meat quality of Qaidam cattle. Conclusions: Therefore, this study provides new genetic insights on the environmental adaptability and important economic traits of Qaidam cattle.How to cite: Guo S, Wu X, Pei J, et al. Genome-wide CNV analysis reveals variants associated with high-altitude adaptation and meat traits in Qaidam cattle. Electron J Biotechnol 2021;54. https://doi.org/10.1016/j.ejbt.2021.07.006

Mammalian spermatozoa are not mature after ejaculation and must undergo additional functional and structural changes within female reproductive tracts to achieve subsequent fertilization, including both capacitation and acrosome reaction (AR), which are dominated by post-translational modifications (PTMs), especially phosphorylation. However, the mechanism of protein phosphorylation during frozen-thawed sperm capacitation and AR has not been well studied. In this study, the phosphoproteomics approach was employed based on tandem mass tag (TMT) labeling combined with liquid chromatography-tandem mass spectrometry (LC-MS/MS) strategy to analyze frozen-thawed sperm in Ashidan yak under three sequential conditions (density gradient centrifugation-based purification, incubation in the capacitation medium and induction of AR processes by the calcium ionophore A23187 treatment). The identification of 1,377 proteins with 5,509 phosphorylation sites revealed changes in phosphorylation levels of sperm-specific proteins involved in regulation of spermatogenesis, sperm motility, energy metabolism, cilium movement, capacitation and AR. Some phosphorylated proteins, such as AKAP3, AKAP4, SPA17, PDMD11, CABYR, PRKAR1A, and PRKAR2A were found to regulate yak sperm capacitation and AR though the cAMP/PKA signaling pathway cascades. Notably, the phosphorylation level of SPA17 at Y156 increased in capacitated sperm, suggesting that it is also a novel functional protein besides AKAPs during sperm capacitation. Furthermore, the results of this study suggested that the phosphorylation of PRKAR1A and PRKAR2A, and the dephosphorylation of CABYR both play key regulatory role in yak sperm AR process. Protein-protein interaction analysis revealed that differentially phosphorylated proteins (AKAP3, AKAP4, FSIP2, PSMD11, CABYR, and TPPP2) related to capacitation and AR process played a key role in protein kinase A binding, sperm motility, reproductive process, cytoskeleton and sperm flagella function. Taken together, these data provide not only a solid foundation for further exploring phosphoproteome of sperm in yak, but an efficient way to identify sperm fertility-related marker phosphorylated proteins.

N6-methyladenosine (m6A) is the most common form of eukaryotic mRNA modification, and it has been shown to exhibit broad regulatory activity in yeast, plants, and mammals. The specific role of m6A methylation as a regulator of spermatogenesis, however, has yet to be established. In this experiment, through a series of preliminary studies and methylated RNA immunoprecipitation sequencing, the m6A map of cattle-yak testicular tissue was established as a means of exploring how m6A modification affects cattle-yak male infertility. Cattle-yak testis tissues used in this study were found to contain sertoli cells and spermatogonia. Relative to sexually mature yak samples, those isolated from cattle-yak testis exhibited slightly reduced levels of overall methylation, although these levels were significantly higher than those in samples from pre-sexually mature yaks. Annotation analyses revealed that differentially methylated peaks were most concentrated in exonic regions, with progressively lower levels of concentration in the 3'-untranslated region (UTR) and 5'-UTR regions. To further explore the role of such m6A modification, enrichment analyses were performed on differentially methylated and differentially expressed genes in these samples. For the cattle-yaks vs. 18-months-old yaks group comparisons, differentially methylated genes were found to be associated with spermatogenesis-related GO terms related to the cytoskeleton and actin-binding, as well as with KEGG terms related to the regulation of the actin cytoskeleton and the MAPK signaling pathway. Similarly, enrichment analyses performed for the cattle-yaks vs. 5-years-old yaks comparison revealed differentially methylated genes to be associated with GO terms related to protein ubiquitination, ubiquitin ligase complexes, ubiquitin-dependent protein catabolism, and endocytotic activity, as well as with KEGG terms related to apoptosis and the Fanconi anemia pathway. Overall, enrichment analyses for the cattle-yaks vs. 18-months-old yaks comparison were primarily associated with spermatogenesis, whereas those for the cattle-yaks vs. 5-years-old yaks comparison were primarily associated with apoptosis.

Abstract Background and aim Mitochondrial genome has aseries of characteristics such as simple structure, no recombination, maternalinheritance, stable structure, fast evolution rate, and high copy number. Moreover, it is easy to be sequenced,contains high-resolution phylogenetic information, and exists in a wide rangeof taxa. Therefore, it is widely used in the study of biological phylogeny. Atpresent, phylogenetic studies focus mainly on D-loop region, cytochrome b gene,and protein-coding sequence. Phylogenetic studies using the mitochondrialcomplete sequence are rarely reported in yak. Therefore, the present studyaimed to construct phylogenetic tree using yak mitochondrial complete sequenceand compare the subsequent results with previous findings obtained usingpartial sequences. Results Complete mitochondrial sequences of five yakpopulations from Qinghai and Xinjiang were obtained. The mitotype diversity ofthe five populations was Xueduo yak (0.992 ± 0.015), Pamir yak (0.990 ± 0.014),Yushu yak (0.963 ± 0.033), Qilian yak (0.948 ± 0.036), and Huanhu yak (0.905 ±0.048), which showed a higher mitotype diversity compared with other breeds fromthe previous reports, including Jiulong yak, Maiwa yak, Zhongdian yak, andTianzhu yak. A total of 78 mitotypes were obtained from 111 individuals. Amongthese, Yushu yak, Huanhu yak, Xueduo yak, and Qilian yak all shared mitotypes,but the Pamir yak did not share mitotypes with these four populations.Phylogenetic analysis showed that yak populations were separable into threedistinct branches. The analysis identified a new phylogenetic branch containingboth wild and domestic yaks. The 155 mitotypes found in 206 individuals weredivided into 3 haplogroups by mitotype clustering. Thehaplogroup was not associated with the geographical distribution of yaks. Theyaks in the same population or the same ecological environment were distributedin different haplogroups. Among the threehaplogroups, haplogroup A and haplogroup B showed a star-shaped distribution ofmitotypes. The central mitotypes were widely distributed and had a highfrequency. Conclusions Thegenetic diversity of yaks in Qinghai was high. Both domestic and wild yaks clusteredinto three branches.

Background: Transmembrane protein 95 (TMEM95) plays a role in male fertility. Previous studies showed that genes with a significant impact on reproductive traits can also affect the growth traits of livestock. Thus, we speculated that the genetic variation of TMEM95 gene may have effects on growth traits of cattle. Results: Two SNPs were genotyped. The rs136174626 and rs41904693 were in the intron 4 and 3′-untranslated region, respectively. The linkage disequilibrium analysis illustrated that these two loci were not linked. The rs136174626 was associated with six growth traits of Nanyang cattle, four traits of Luxi cattle, and three traits of Ji’an cattle. For rs41904693 locus, the GG individuals had greater body height and abdominal girth in Ji’ an cattle than TT and TG individuals. In Jinnan cattle, GG and TT individuals had greater body height, height at hip cross, body length, and heart girth than TG individuals. The potential splice site prediction results suggest that the rs136174626 may influence the splicing efficiency of TMEM95, and the miRNA binding site prediction results showed that the rs41904693 may influence the expression of TMEM95 by affecting the binding efficiency of Bta-miR-1584 and TMEM95 3′-UTR. Conclusions: The findings of the study suggested that the two SNPs in TMEM95 could be a reliable basis for molecular breeding in cattle.How to cite: Guo X, Zhang S, Yang H, et al. Bovine TMEM95 gene: Polymorphisms detecting in five Chinese indigenous cattle breeds and their association with growth traits. Electron J Biotechnol 2021;51. https://doi.org/10.1016/j.ejbt.2021.03.004

Abstract The current study aimed to investigate the damage of long-term high concentrate diet feeding pattern on Yak jejunal structure, physiological function and protein composition during cold season. Twelve Datong male Yak (Bos grunniens) with the same age from cold season were randomly selected and slaughtered to determine Yak jejunal digestive enzyme activity, morphology and protein composition on different feeding patterns in Tibetan Plateau. The results showed that Yak jejunum digestive enzyme activity and morphology of grazing reared group were better than those in the intensively reared group. A total of 96 differentially expressed proteins were identified by label-free Mass Spectrometry (MS), which could be concluded to two predominant themes: protein structure and inflammatory response. Nine differentially expressed proteins were correlated in Yak jejunum damage in different feeding patterns. According to this research, we found that feeding pattern resulted the differences in Yak jejunum physiological function, morphology and protein composition. This fact was confirmed long-term high dietary concentrate feeding could damage the jejunum epithelial morphology and function.

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.

Epididymis development is the basis of male reproduction and is a crucial site where sperm maturation occurs. In order to further understand the epididymal development of yak and how to regulate sperm maturation, we conducted a multi-omics analysis. We detected 2274 differential genes, 222 differential proteins and 117 co-expression genes in the cauda epididymis of yak before and after sexual maturity by RNA-seq and proteomics techniques, which included TGFBI, COL1A1, COL1A2, COL3A1, COL12A1, SULT2B1, KRT19, and NPC2. These high abundance genes are mainly related to cell growth, differentiation, adhesion and sperm maturation, and are mainly enriched via extracellular matrix receptor interaction, protein differentiation and absorption, and lysosome and estrogen signaling pathways. The abnormal expression of these genes may lead to the retardation of epididymal cauda development and abnormal sperm function in yak. In conclusion, through single and combined analysis, we provided a theoretical basis for the development of the yak epididymal cauda, sperm maturation, and screening of key genes involved in the regulation of male yak reproduction.

The utilization of yak milk is still in a primary stage, and the nutrition composition of yak colostrum is not systematically characterized at present. In this study, the lipids, fatty acids, amino acids and their derivatives, metabolites in yak colostrum, and mature milk were detected by the non−targeted lipidomics based on (ultra high performance liquid chromatography tandem quadrupole mass spectrometer) UHPLC−MS, the targeted metabolome based on gas chromatography−mass spectrometer (GC−MS), the targeted metabolome analysis based on UHPLC−MS, and the non-targeted metabolome based on ultra high performance liquid chromatography tandem quadrupole time of flight mass spectrometer (UHPLC−TOF−MS), respectively. Meanwhile, the nutrition composition of yak colostrum was compared with the data of cow mature milk in the literatures. The results showed that the nutritive value of yak colostrum was higher by contrast with yak and cow mature milk from the perspective of the fatty acid composition and the content of Σpolyunsaturated fatty acids (PUFAs), Σn−3PUFAs; the content of essential amino acid (EAA) and the ratio of EAA/total amino acid (TAA) in yak colostrum were higher than the value in yak mature milk; and the content of functional active lipids including phosphatidylcholines (PC), phosphatidylglycerol (PG), phosphatidylserine (PS), lyso-phosphatidylcholine (LPC), lyso-phosphatidylglycerol (LPG), lyso-phosphatidylinositol (LPI), sphingomyelin (SM), ganglioside M3 (GM3), ganglioside T3 (GT3), and hexaglycosylceramide (Hex1Cer) in yak colostrum, was higher than the value of yak mature milk. Moreover, the differences of nutritive value between yak colostrum and mature milk were generated by the fat, amino acids and carbohydrate metabolism that were regulated by the ovarian hormone and referencesrenin-angiotensin-aldosterone system in yaks. These research results can provide a theoretical basis for the commercial product development of yak colostrum.

N6-methyladenine (m6A) RNA undergoes epigenetic modification, which is the most extensive intermediate chemical modification in mRNA. Although this modification occurs in all living organisms, it is the most widely studied among mammals. However, to date, no study has investigated the m6A transcriptome-wide map of yak and its potential biological functions in muscle development. In this study, the differences of m6A methylation and gene expression in yak muscle development belonging to three age groups, namely 3 years (group A), 6 months (group M), and 90-day-old fetuses (group E), were determined by using methylated RNA immunoprecipitation sequencing (MeRIP-Seq) and RNA sequencing (RNA-Seq). In these three groups, a total of 6,278 (A), 9,298 (E), and 9,584 (M) m6A peaks were identified, with average densities between 1.02 and 2.01. m6A peaks were mostly enriched in the stop codon, 3′ untranslated region (UTR) region, and inner long exon region with consensus motifs of UGACA. In all the three stages, the m6A peak enrichment level was significantly negatively correlated with mRNA abundance (Pearson's correlation coefficient r = −0.22 to −0.32, p < 10−16). The functional enrichment of genes consistently modified by m6A methylation, particularly those genes that regulate cell differentiation as well as muscle growth and development, was observed at all three stages. Moreover, m6A abundance was negatively associated with gene expression levels, indicating that m6A might play a vital role in modulating gene expression during yak muscle development. This comprehensive map thus provides a solid foundation for determining the potential functional role of m6A RNA modification in yak muscle growth.

Tianzhu white yak is a rare local yak breed with a pure white coat in China. In recent years, breeders have discovered long-haired individuals characterized by long hair on the forehead in the Tianzhu white yak, and the length and density of the hair on these two parts of the body are higher than that of the normal Tianzhu white yak. To elucidate the genetic mechanism of hair length in Tianzhu white yak, we re-sequence the whole genome of long-haired Tianzhu White yak (LTWY) (n = 10) and normal Tianzhu White yak (NTWY) (n = 10). Then, fixation index (FST), θπ ratio, cross-population composite likelihood ratio (XP-CLR), integrated haplotype score (iHS), cross-population extended haplotype homozygosity (XP-EHH), and one composite method, the de-correlated composite of multiple signals (DCMS) were performed to discover the loci and genes related to long-haired traits. Based on five single methods, we found two hotspots of 0.2 and 1.1 MB in length on chromosome 6, annotating two (FGF5, CFAP299) and four genes (ATP8A1, SLC30A9, SHISA3, TMEM33), respectively. Function enrichment analysis of genes in two hotspots revealed Ras signaling pathway, MAPK signaling pathway, PI3K-Akt signaling pathway, and Rap1 signaling pathway were involved in the process of hair length differences. Besides, the DCMS method further found that four genes (ACOXL, PDPK1, MAGEL2, CDH1) were associated with hair follicle development. Henceforth, our work provides novel genetic insights into the mechanisms of hair growth in the LTWY.

Abstract Background Long non-coding RNA (lncRNA) as an important regulator has been demonstrated playing an indispensable role in the biological process of hair follicles (HFs) growth. However, their function and expression profile in the HFs cycle of yak are yet unknown. Only a few functional lncRNAs have been identified, partly due to the low sequence conservation and lack of identified conserved properties in lncRNAs. Here, lncRNA-seq was employed to detect the expression profile of lncRNAs during the HFs cycle of yak, and the sequence conservation of two datasets between yak and cashmere goat during the HFs cycle was analyzed. Results A total of 2884 lncRNAs were identified in 5 phases (Jan., Mar., Jun., Aug., and Oct.) during the HFs cycle of yak. Then, differential expression analysis between 3 phases (Jan., Mar., and Oct.) was performed, revealing that 198 differentially expressed lncRNAs (DELs) were obtained in the Oct.-vs-Jan. group, 280 DELs were obtained in the Jan.-vs-Mar. group, and 340 DELs were obtained in the Mar.-vs-Oct. group. Subsequently, the nearest genes of lncRNAs were searched as the potential target genes and used to explore the function of DELs by GO and KEGG enrichment analysis. Several critical pathways involved in HFs development such as Wnt signaling pathway, VEGF signaling pathway, and signaling pathways regulating pluripotency of stem cells, were enriched. To further screen key lncRNAs influencing the HFs cycle, 24 DELs with differ degree of sequence conservation were obtained via a comparative analysis of partial DELs with previously published lncRNA-seq data of cashmere goat in the HFs cycle using NCBI BLAST-2.9.0+, and 3 DELs of them were randomly selected for further detailed analysis of the sequence conservation properties. Conclusions This study revealed the expression pattern and potential function of lncRNAs during HFs cycle of yak, which would expand the knowledge about the role of lncRNAs in the HFs cycle. The findings related to sequence conservation properties of lncRNAs in the HFs cycle between the two species may provide valuable insights into the study of lncRNA functionality and mechanism.

Abstract Background Mammalian hair play an important role in mammals’ ability to adapt to changing climatic environments. The seasonal circulation of yak hair helps them adapt to high altitude but the regulation mechanisms of the proliferation and differentiation of hair follicles (HFs) cells during development are still unknown. Here, using time series data for transcriptome and hormone contents, we systematically analyzed the mechanism regulating the periodic expression of hair development in the yak and reviewed how different combinations of genetic pathways regulate HFs development and cycling. Results This study used high-throughput RNA sequencing to provide a detailed description of global gene expression in 15 samples from five developmental time points during the yak hair cycle. According to clustering analysis, we found that these 15 samples could be significantly grouped into three phases, which represent different developmental periods in the hair cycle. A total of 2316 genes were identified in these three consecutive developmental periods and their expression patterns could be divided into 9 clusters. In the anagen, genes involved in activating hair follicle growth are highly expressed, such as the WNT pathway, FGF pathway, and some genes related to hair follicle differentiation. In the catagen, genes that inhibit differentiation and promote hair follicle cell apoptosis are highly expressed, such as BMP4, and Wise. In the telogen, genes that inhibit hair follicle activity are highly expressed, such as DKK1 and BMP1. Through co-expression analysis, we revealed a number of modular hub genes highly associated with hormones, such as SLF2, BOP1 and DPP8. They may play unique roles in hormonal regulation of events associated with the hair cycle. Conclusions Our results revealed the expression pattern and molecular mechanisms of the seasonal hair cycle in the yak. The findings will be valuable in further understanding the alpine adaptation mechanism in the yak, which is important in order to make full use of yak hair resources and promote the economic development of pastoral plateau areas.

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.

Objective As an iconic symbol of Qinghai-Tibetan Plateau and of high altitude, yak are subjected to hypoxic conditions that challenge aerobic metabolism. Matrix metalloproteinases-3 (MMP3) is assumed to be a key target gene of hypoxia-inducible factor-1α that function as a master regulator of the cellular response to hypoxia. Therefore, the aim of this investigation was to identify the DNA polymorphism of MMP3 gene in domestic yak and to explore its possible association with high-altitude adaptation. Methods The single-nucleotide polymorphisms (SNPs) genotyping and mutations scanning at the MMP3 locus were conducted in total of 344 individuals from four domestic Chinese yak breeds resident at different altitudes on the Qinghai-Tibetan Plateau, using high-resolution melting analysis and DNA sequencing techniques. Results The novel of SNPs rs2381 A→G and rs4331 C→G were identified in intron V and intron VII of MMP3, respectively. Frequencies of the GG genotype and the G allele of SNP rs2381 A→G observed in high-altitude Pali yak were significantly higher than that of the other yak breeds resident at middle or low altitude (p0.05). Haplotype GC was the dominant among the 4 yak breeds, and Pearson correlation analysis showed that the frequencies of GC was significantly lower in Ganan (GN), Datong (DT), and Tianzhu white yaks (TZ) compared with Pali (PL) yak. The two SNPs were in moderate linkage disequilibrium in high-altitude yaks (PL) but not in middle-altitude (GN, DT) and low-altitude (TZ) yaks. Conclusion These results indicate that MMP3 may have been subjected to positive selection in yak, especially that the SNP rs2381 A→G mutation and GC haplotypes might contribute to adaptation for yak in high-altitude environments.

Yak (Bos grunniens) is considered an iconic symbol of Tibet and high altitude, but they suffer from malnutrition during the cold season that challenges the metabolism of energy. Adipocytes perform a crucial role in maintaining the energy balance, and adipocyte differentiation is a complex process involving multiple changes in the expression of genes. N6-methyladenosine (m6A) plays a dynamic role in post-transcription gene expression regulation as the most widespread mRNA modification of the higher eukaryotes. However, currently there is no research existing on the m6A transcriptome-wide map of bovine animals and their potential biological functions in adipocyte differentiation. Therefore, we performed methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) to determine the distinctions in m6A methylation and gene expression during yak adipocyte differentiation. In yak adipocyte and preadipocyte the content of m6A and m6A-associated enzymes was substantially different. In the two groups, a total of 14,710 m6A peaks and 13,388 m6A peaks were identified. For the most part, m6A peaks were enriched in stop codons, 3′-untranslated regions, and coding regions with consensus motifs of GGACU. The functional enrichment exploration displayed that differentially methylated genes participated in some of the pathways associated with adipogenic metabolism, and several candidate genes (KLF9, FOXO1, ZNF395, and UHRF1) were involved in these pathways. In addition to that, there was a positive association between m6A abundance and levels of gene expression, which displayed that m6A may play a vital role in modulating gene expression during yak adipocyte differentiation. Further, in the adipocyte group, several methylation gene protein expression levels were significantly higher than in preadipocytes. In short, it can be concluded that the current study provides a comprehensive explanation of the m6A features in the yak transcriptome, offering in-depth insights into m6A topology and associated molecular mechanisms underlying bovine adipocyte differentiation, which might be helpful for further understanding its mechanisms.

Cattle-yak, as the hybrid offspring of cattle (Bos taurus) and yak (Bos grunniens), demonstrates obvious heterosis in production performance. Male hybrid sterility has been focused on for a long time; however, the mRNAs and non-coding RNAs related to muscle development as well as their regulatory networks remain unclear. The phenotypic data showed that the production performance (i.e., body weight, withers height, body length, and chest girth) of cattle-yak was significantly better than that of the yak, and the economic benefits of the cattle-yak were higher under the same feeding conditions. Then, we detected the expression profiles of the longissimus dorsi muscle of cattle-yak and yak to systematically reveal the molecular basis using the high-throughput sequencing technology. Here, 7,126 mRNAs, 791 lncRNAs, and 1,057 circRNAs were identified to be differentially expressed between cattle-yaks and yaks in the longissimus dorsi muscle. These mRNAs, lncRNA targeted genes, and circRNA host genes were significantly enriched in myoblast differentiation and some signaling pathways related to muscle development (such as HIF-1 signaling pathway and PI3K-Akt signaling pathway). We constructed a competing endogenous RNA (ceRNA) network and found that some non-coding RNAs differentially expressed may be involved in the regulation of muscle traits. Taken together, this study may be used as a reference tool to provide the molecular basis for studying muscle development.

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.

The potential regulatory role of N6-methyladenosine (m6A), the most prominent mRNA modification in eukaryotes, has recently been identified in mammals, plants, and yeast. However, whether and how m6A methylation is involved in sexual maturation in mammals remains largely unexplored. In this study, testicular tissue was obtained from yaks before and after sexual maturation, and m6A maps were generated via preliminary experiments and methylated RNA immunoprecipitation sequencing. Only spermatogonial cells and a few primary spermatocytes were observed in the testicular tissue of yaks before sexual maturation, while spermatogenic cells at different stages of maturity could observed after sexual maturation. Experiments examining the expression of methylation-related enzymes and overall methylation levels showed that the methylation levels in yak testes increased after sexual maturation. Overall, 1,438 methylation peaks were differentially expressed before and after sexual maturation; 1,226 showed significant up-regulation and 212 showed significant down-regulation after sexual maturation. Annotation analysis showed that the differential methylation peaks were most commonly concentrated in the exon region, followed by the 3′UTR and finally the 5′UTR region. KEGG pathway analysis demonstrated that homologous recombination, the Notch signaling pathway, growth hormone synthesis, and other signaling pathways may be involved in testicular development and maturation in yaks. Levels of most m6A modifications were positively correlated with mRNA abundance, suggesting that m6A plays a regulatory role in mammalian sexual maturation. To our knowledge, this is the first report of an m6A transcriptional map of the yak testes, and our study lays the foundation for elucidating the function of m6A in the development of yak testes.

Copy number variation (CNV) is a fundamental type of structural variation of the genome affecting the economic traits of livestock. The SOX6 gene (sex-determining region Y-box 6), as a transcription factor, has multiple functions with regard to sex determination, embryonic growth, the nervous system development, as well as bone, and various organ formation. This study employed quantitative real-time fluorescence quota PCR (qPCR) for detecting the SOX6-CNV of the 311 Ashidan yaks and analyzed the correlation of the SOX6-CNV with four phenotypes (including body weight, withers height, body length, and chest girth) of the yaks aged 6, 12, 18, and 30 months using ANOVA and multiple comparisons. Furthermore, the SOX6 gene expression was identified in seven different tissues of the yaks. The experiment results demonstrated the expression of SOX6 in each tissue, and the kidney and muscle tissue were found to have higher relative expression levels. Based on the processing by IBM SPSS software, SOX6-CNV was significantly correlated with the chest girth of the 6-months old yaks (p < 0.05) and 30-months yaks (p < 0.05), and withers height of 6 months yaks (p < 0.05) and 18-months yaks (p < 0.05), as well as the normal type of CNV, was chosen for yak breeding. In conclusion, SOX6 might be prominently involved in promoting growth and development of yaks, suggesting that the SOX6 gene can be used in breeding yaks by molecular marker-assisted selection (MAS). The study also offered some important insights into the references and clues for the genetic breeding of yaks.

The Jeryak is the hybrid offspring of yaks and Jersey cattle and exhibit improved milk and meat yields. Biomolecules carried within milk exosomes are important for cell growth, development, immune regulation, and various pathophysiological processes. Previous studies showed that miRNAs regulate mammary gland development, lactation, and milk quality. This study explored the relationship between milk exosomes miRNAs and lactation performance. A comparison of the milk content showed that yak milk was of a better quality compared to Jeryak milk (casein, fat, TS, SNF, lactose). Milk collected in December was superior to that collected in June for both yak and Jeryak, except for lactose concentrations. Exosomes were extracted by density gradient centrifugation and miRNA expression profiles in milk exosomes from three yaks and three Jeryaks collected in June and December were detected by small RNA sequencing. In all, 22, 120, 78, and 62 differentially expressed miRNAs (DEMs) were identified in Jun_ JY vs. Jun_ Y (P1: Jeryak in June vs. Yak in June), Jun_ JY vs. Dec_ JY (P2: Jeryak in June vs. Jeryak in December), Dec_ JY vs. Dec_ Y (P3: Jeryak in December vs. Yak in December), and Jun_ Y vs. Dec_ Y (P4: Yak in June vs. Yak in December) groups. These DEMs were enriched in functions and signaling pathways related to lactation performance. In conclusion, these findings are a reference tool to study the molecular basis of lactation performance.

Cattle-yak is a hybrid F1 generation of cattle and yak, which has a history of more than 3000 years and has shown better production performance and higher economic benefits than those of yaks. However, up to now, there has been no study on the transcriptome-wide m6A methylation profile of bovine skeletal muscle and its potential biological function during muscle development. Here, we observed significant changes in the expression levels of muscle-related marker genes and methylation-related enzymes during the development of cattle-yak, and the overall m6A content in the Longissimus dorsi muscle of 18-month-old cattle-yak decreased significantly. A total of 36,602 peaks, 11,223 genes and 8388 lncRNAs were identified in the two groups, including 2989 differential peaks (427 up-regulated peaks and 2562 down-regulated peaks), 1457 differentially expressed genes (833 up-regulated genes and 624 down-regulated genes) and 857 differentially expressed lncRNAs (293 up-regulated lncRNAs and 564 down-regulated lncRNAs). GO and KEGG analysis revealed that they were significantly enriched in some muscle-related pathways (Wnt signaling pathway and MAPK signaling pathway) and high-altitude adaptation-related pathway (HIF-1 signaling pathway). Moreover, m6A abundance was positively correlated with gene expression levels, while it was negatively correlated with lncRNA expression levels. This indicates that m6A modification played an important role in the Longissimus dorsi muscle development of cattle-yak; however, the regulation mechanism of m6A-modified mRNA and lncRNA may be different. This study was the first report of transcriptome-wide m6A-modified mRNAs and lncRNAs atlas in the Longissimus dorsi muscle development of cattle-yak, one which will provide new perspectives for genetic improvement in bovines.

Copy number variations (CNVs) are a result of genomic rearrangement affecting DNA regions over 1 kb in length, and can include inversions, translocations, deletions, and duplications. The molecule interacting with CasL-like protein 2 (MICALL2) gene is primarily associated with mitochondrial protein targeting and exhibits predicted stress fiber colocalization. The monoacylglycerol O-acyltransferase 2 (MOGAT2) gene encodes an enzyme responsible for catalyzing diacylglycerol synthesis from 2-monoacylglycerol and fatty acyl-CoA. For this study, blood samples were obtained from 315 yaks, and the body weight, body length, withers height, and chest girth of these animals were measured at 6, 12, 18, and 30 months of age. Genomic DNA was harvested from the collected blood samples, and CNVs in these samples were detected by qPCR. The resultant data were compared using ANOVAs, revealing significant associations between MICALL2 gene CNVs and body weight at 6 months of age (p < 0.05), body length and chest girth at 30 months of age (p < 0.05), and withers height at 18 months of age (p < 0.01) in Ashidan yaks. Similarly, MOGAT2 CNVs were significantly associated with body weight at 6 and 30 months of age (p < 0.05), and with withers height at 18 months of age (p < 0.01) in these Ashidan yaks. MICALL2 and MOGAT2 gene expression was further analyzed in yak tissue samples, revealing that MICALL2 was most highly expressed in the adipose tissue, whereas MOGAT2 was most highly expressed in the lung. These results thus confirmed the relationship between CNVs in the MICALL2 and MOGAT2 genes and Ashidan yak growth traits, providing a valuable gene locus that can be leveraged for future marker-assisted yak breeding efforts.

Kecai yaks are regarded as an important genetic resource in China owing to their high fecundity and flavorful meat. However, the genetic characteristics of Kecai yaks have not been effectively characterized to date, and the relationship between Kecai yaks and other yak breeds remains to be fully characterized. In this paper, the resequencing of the Kecai yak genome is performed leading to the identification of 11,491,383 high-quality single nucleotide polymorphisms (SNPs). Through principal component, phylogenetic, and population genetic structure analyses based on these SNPs, Kecai yaks were confirmed to represent an independent population of yaks within China. In this study, marker and functional enrichment analysis of genes related to positive selection in Kecai yak was carried out, and the results show that such selection in Kecai yaks is associated with the adaptation to alpine environments and the deposition of muscle fat. Overall, these results offer a theoretical foundation for the future utilization of Kecai yak genetic resources.

Single-cell sequencing technology can fully reflect the heterogeneity of cell populations at the single cell level, making it possible for us to re-recognize various tissues and organs. At present, the sequencing study of hair follicles is transiting from the traditional ordinary transcriptome level to the single cell level, which will provide diverse insights into the function of hair follicle cells. This review focuses on research advances in the hair follicle microenvironment obtained from scRNA-seq studies of major cell types in hair follicle development, with a special emphasis on the discovery of new subpopulations of hair follicles by single-cell techniques. We also discuss the problems and current solutions in scRNA-seq observation and look forward to its prospects.

The plateau adaptability and stress resistance of yaks are widely known based on their capacity to survive under severe habitat conditions. However, a few studies on brain mitochondria have characterized these adaptations at the protein level. We identified and quantified the brain mitochondrial proteins using isobaric tags for relative and absolute quantification (iTRAQ) and Proteomics. Western blotting was used to verify changes in the expression of target proteins. A total of 57 differentially abundant proteins (DAPs) were identified in the yak brain tissue. Gene Ontology (GO) analysis showed molecular functions of these DAPs including downregulated oxidoreductase activity but upregulated coenzyme binding. Significantly enriched biological processes were oxidation–reduction process (downregulated) and small molecule metabolic processes (upregulated). STRING protein interaction analysis indicated a complex interaction between dehydrogenase, transaminase, and ATP synthetase families. Reactome pathway analysis highlighted that the majority of DAPs participated in aerobic metabolic pathways such as metabolism, citric acid cycle, and respiratory electron transport. Immunoblotting confirmed that changes in FKBP4 and ATPAF2 expression were consistent with the results of mass spectrometry. We performed a high-throughput screening to identify DAPs in brain mitochondria between yak and cattle, which could explain the plateau adaptability of yaks.

The change in the feeding system can greatly improve the growth performance of the yak (Bos grunniens), an important livestock species in the plateau region. Here, we comprehensively compared the effects of different feeding systems on the growth performance and ruminal development of yaks, and investigated the effects of ruminal microorganisms and metabolites using the 16S rRNA gene sequencing and liquid chromatograph–mass spectrometer (LC-MS) technologies. We found that compared to traditional grazing feeding, house feeding significantly improved the growth performance (such as average daily gain and net meat weight) and rumen development of the yaks. At the genus level, the abundance of Rikenellaceae RC9 Gut group, Christensenellaceae R-7 group, Lachnospiraceae NK3A20 group, Ruminococcaceae UCG-014, and Prevotellaceae UCG-003 showed significant differences and was closely related to rumen development in the two distinct feeding systems. Also, metabolomics revealed that the change in the feeding system significantly affected the concentration and metabolic pathways of the related rumen metabolites. The metabolites with significant differences were significantly enriched in purine metabolism (xanthine, adenine, inosine, etc.), tyrosine metabolism (L-tyrosine, dopaquinone, etc.), phenylalanine metabolism (dihydro-3-caumaric acid, hippuric acid, etc.), and cAMP signaling pathway [acetylcholine, (-)-epinephrine, etc.]. This study scientifically support the house fattening feeding system for yaks. Also, our results provide new insights into the composition and function of microbial communities that promote ruminal development and in general growth of the yaks.

In this study, high-throughput Illumina sequencing was employed to assemble the complete mitochondrial genome of the Meiren yak (Bos grunniens), a local yak breed from Gansu Province, China. The mitochondrial genome is 16,321 bp long with an A + T-biased nucleotide composition and harbors 13 protein-coding, 22 Trna, and 2 rRNA genes, and a noncoding control region. The mitogenomic organization and codon usage are highly similar to those of previously published congeneric mitochondrial genomes. Bayesian phylogenetic analysis indicates that Meiren yak is most closely related to nine other yak breeds (incl. Datong, Huanhu, Pali, Pamir, Polled, Qilian, Seron, Sunan, and Tianjun yaks).

Abstract Background Copy number variations (CNVs), which are genetic variations present throughout mammalian genomes, are a vital source of phenotypic variation that can lead to the development of unique traits. In this study we used the Illunima BovineHD BeadChip to conduct genome-wide detection of CNVs in 215 polled yaks. Results A total of 1066 CNV regions (CNVRs) were detected with a total length of 181.6 Mb, comprising ~ 7.2% of the bovine autosomal genome. The size of these CNVRs ranged from 5.53 kb to 1148.45 kb, with an average size of 170.31 kb. Eight out of nine randomly chosen CNVRs were successfully validated by qPCR. A functional enrichment analysis of the CNVR-associated genes indicated their relationship to a number of molecular adaptations that enable yaks to thrive at high altitudes. One third of the detected CNVRs were mapped to QTLs associated with six classes of economically important traits, indicating that these CNVRs may play an important role in variations of these traits. Conclusions Our genome-wide yak CNV map may thus provide valuable insights into both the molecular mechanisms of high altitude adaptation and the potential genomic basis of economically important traits in yak.

Copy number variation (CNV) is a structural variant with significant impact on genetic diversity. CNV has been widely used in breeding for growth traits, meat production or quality, and coat color. SRY-like box genes (SOXs) are a class of transcription factors that play a regulatory role in cell fate specification and differentiation. SOX5 and SOX8 belong to subgroups D and E of the SOXs, respectively. Previous studies have shown that SOX5 and SOX8 are essential in the development of bones. In this study, we explored the association between the growth traits and CNVs of SOX5 and SOX8 in 326 Ashidan yaks and detected mRNA expression levels in different tissues. Our results illustrated that CNVs of SOX5 and SOX8 were significantly associated with withers height at 18 months of age and chest girth at 30 months of age (p < 0.05). The CNV combination of SOX5 and SOX8 was significantly associated with withers height at 18 months of age (p < 0.01). SOX5 expression in the lung was significantly higher than in the heart, spleen, kidney, and muscle (p < 0.05). SOX8 expression in the lung was significantly higher than in the liver and muscle (p < 0.05). Our results provide evidence that the CNVs of SOX5 and SOX8 genes could be used as new markers for the selection of yak growth traits.