Living at a high altitude involves many environmental challenges. The combined effects of hypoxia and cold stress impose severe physiological challenges on endothermic animals. The yak is integral to the livelihood of the people occupying the vast, inhospitable Qinghai–Tibetan plateau and the surrounding mountainous region. Due to long-term selection, the yak exhibits stable and unique genetic characteristics which enable physiological, biochemical, and morphological adaptations to a high altitude. Thus, the yak is a representative model for mammalian plateau-adaptability studies. Understanding coping mechanisms provides unique insights into adaptive evolution, thus informing the breeding of domestic yaks. This review provides an overview of genetic adaptations in Bos grunniens to high-altitude environmental stress. Combined genomics and theoretical advances have informed the genetic basis of high-altitude adaptations.

Yak reproductive characteristics have received extensive attention, though the molecular regulation mechanism of its ovarian activity remains to be explored. Therefore, this study initially conducted a comparative analysis of yak ovarian activities in anestrus, estrus, and pregnancy regarding their morphology and histology, followed by implementing RNA sequencing (RNA-seq) technology to detect the overall gene expression and biological mechanism in different reproductive stages. H&E staining showed that there were more growing follicles and mature follicles in ovarian tissue sections during estrus than ovarian tissues during non-estrus. The RNA-seq analysis of yak ovary tissues in three periods showed that DEGs related to follicular development and hormone metabolism were screened in the three comparison groups, such as COL1A2, NR4A1, THBS2, PTGS2, SCARB1, STAR, and WNT2B. Bioinformatics analysis showed that these DEGs are involved in ion binding, cell development, metabolic processes, enriched in ECM–receptor interactions, steroid biosynthesis, together with aldosterone generation/discharge and Wnt/PI3K-Akt signaling pathways. In addition, we speculate alternate splice development events to have important role/s in regulating ovarian functional genomic expression profiles. These results provide essential knowledge aimed at scrutinizing pivotal biomarkers for yak ovarian activity, together with paving the way for enhancing researchers’ focus on improving yak reproductive performance.

The yak (Bos grunniens) is subjected to nutritional deficiency during the whole winter grazing season; deciphering the adipose metabolism and energy homeostasis under cold and nutrients stress conditions could be a novel way to understand the specific mechanism of energy metabolism. Circular RNAs (circRNAs) have elucidated that they play a key role in many biological events, but the regulatory function of adipose development remains mostly unknown. Therefore, the expression pattern of circRNAs were identified for the first time during yak adipocyte differentiation to gain insight into their potential functional involvement in bovine adipogenesis. We detected 7203 circRNA candidates, most of them contained at least two exons, and multiple circRNA isoforms could be generated from one parental gene. Analysis of differential expression circRNAs displayed that 136 circRNAs were differentially expressed at day 12 (Ad) after adipocyte differentiation, compared with the control at day 0 (Pread 0), while 7 circRNAs were detected on day 2. Sanger sequencing validated that six circRNAs had head-to-tail junction, and quantitative real-time PCR (qPCR) results revealed that the expression patterns of ten circRNAs were consistent with their expression levels from RNA-sequencing (RNA-seq) data. We further predicted the networks of circRNA-miRNA-gene based on miRNAs sponging by circRNAs, in which genes were participated in the adipocyte differentiation-related signaling pathways. After that, we constructed several adipocyte differentiation-related ceRNAs and revealed six circRNAs (novel_circ_0009127, novel_circ_0000628, novel_circ_0011513, novel_circ_0010775, novel_circ_0006981 and novel_circ_0001494) were related to adipogenesis. Furthermore, we analyzed the homology among yak, human and mouse circRNAs and found that 3536 yak circRNAs were homologous to human and mouse circRNAs. In conclusion, these findings provide a solid basis for the investigation of yak adipocyte differentiation-related circRNAs and serve as a great reference to study the energy metabolism of high-altitude animals.

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

DNA methylation modifications are implicated in many biological processes. As the most common epigenetic mechanism DNA methylation also affects muscle growth and development. The majority of previous studies have focused on different varieties of yak, but little is known about the epigenetic regulation mechanisms in different age groups of animals. The development of muscles in the different stages of yak growth remains unclear. In this study, we selected the longissimus dorsi muscle tissue at three different growth stages of the yak, namely, 90-day-old fetuses (group E), six months old (group M), and three years old (group A). Using RNA-Seq transcriptome sequencing and methyl-RAD whole-genome methylation sequencing technology, changes in gene expression levels and DNA methylation status throughout the genome were investigated during the stages of yak development. Each group was represented by three biological replicates. The intersections of expression patterns of 7694 differentially expressed genes (DEGs) were identified (padj < 0.01, |log2FC| > 1.2) at each of the three developmental periods. Time-series expression profile clustering analysis indicated that the DEGs were significantly arranged into eight clusters which could be divided into two classes (padj < 0.05), class I profiles that were downregulated and class II profiles that were upregulated. Based on this cluster analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that DEGs from class I profiles were significantly (padj < 0.05) enriched in 21 pathways, the most enriched pathway being the Axon guidance signaling pathway. DEGs from the class II profile were significantly enriched in 58 pathways, the pathway most strongly enriched being Metabolic pathway. After establishing the methylation profiles of the whole genomes, and using two groups of comparisons, the three combinations of groups (M-vs.-E, M-vs.-A, A-vs.-E) were found to have 1344, 822, and 420 genes, respectively, that were differentially methylated at CCGG sites and 2282, 3056, and 537 genes, respectively, at CCWGG sites. The two sets of data were integrated and the negative correlations between DEGs and differentially methylated promoters (DMPs) analyzed, which confirmed that TMEM8C, IGF2, CACNA1S and MUSTN1 were methylated in the promoter region and that expression of the modified genes was negatively correlated. Interestingly, these four genes, from what was mentioned above, perform vital roles in yak muscle growth and represent a reference for future genomic and epigenomic studies in muscle development, in addition to enabling marker-assisted selection of growth traits.

Cattle-yak is a hybrid F₁ 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 m⁶A 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 m⁶A 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, m⁶A abundance was positively correlated with gene expression levels, while it was negatively correlated with lncRNA expression levels. This indicates that m⁶A modification played an important role in the Longissimus dorsi muscle development of cattle-yak; however, the regulation mechanism of m⁶A-modified mRNA and lncRNA may be different. This study was the first report of transcriptome-wide m⁶A-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.