Yak (Bos grunniens) is an unique ruminant species in the Qinghai-Tibetan Plateau (QTP). The ruminant gastrointestinal tract (GIT) microbiota is not only associated with the nutrients metabolism, but also contributes to the host’s local adaptation. Examining the microbiota between cattle and yak in different geography could improve our understanding about the role of microbiota in metabolism and adaptation. To this end, we compared the microbiota in rumen, reticulum, omasum, and abomasum of dairy cattle, yellow cattle, and three yak herds (WQ yak, SZ yak, and ZB yak) lived in different altitude, based on sequencing the bacterial 16S rRNA gene on Illumina Miseq. The bacterial diversity was significantly different among five breeds, whereas the difference among four stomach regions is limited. The phyla Bacteroidetes and Firmicutes were the dominated bacteria regardless of breeds and regions. The nonmetric multidimensional scaling (NMDS) results showed that the microbiota in dairy cattle, yellow cattle and WQ yak significantly differed from that in SZ yak and ZB yak for all four stomach compartments. Canonical correlation analysis revealed that Prevotella and Succiniclasticum spp. were abundant in dairy cattle, yellow cattle and WQ yak, whereas the Christensenellaceae R7 group and the Lachnospiraceae UCG 008 group were prevalent in SZ yak and ZB yak. Moreover, the microbiota in WQ yak was significantly different from that in SZ yak and ZB yak, which were characterized by the higher relative abundance Romboutsia spp., Eubacterium coprostanoligenes, Acetobacter spp., Mycoplasma spp., and Rikenellaceae RC9 group. Overall, these results improves our knowledge about the GIT microbiota composition of QTP ruminant.

Genomic structural variation represents a source for genetic and phenotypic variation, which may be subject to selection during the environmental adaptation and population differentiation. Here, we described a genome-wide analysis of copy number variations (CNVs) in 16 populations of yak based on genome resequencing data and CNV-based cluster analyses of these populations. In total, we identified 51,461 CNV events and defined 3174 copy number variation regions (CNVRs) that covered 163.8 Mb (6.2%) of yak genome with more “loss” events than both “gain” and “both” events, and we confirmed 31 CNVRs in 36 selected yaks using quantitative PCR. Of the total 163.8 Mb CNVR coverage, a 10.8 Mb region of high-confidence CNVRs directly overlapped with the 52.9 Mb of segmental duplications, and we confirmed their uneven distributions across chromosomes. Furthermore, functional annotation indicated that the CNVR-harbored genes have a considerable variety of molecular functions, including immune response, glucose metabolism, and sensory perception. Notably, some of the identified CNVR-harbored genes associated with adaptation to hypoxia (e.g., DCC, MRPS28, GSTCD, MOGAT2, DEXI, CIITA, and SMYD1). Additionally, cluster analysis, based on either individuals or populations, showed that the CNV clustering was divided into two origins, indicating that some yak CNVs are likely to arisen independently in different populations and contribute to population difference. Collectively, the results of the present study advanced our understanding of CNV as an important type of genomic structural variation in yak, and provide a useful genomic resource to facilitate further research on yak evolution and breeding.

The yak is a valuable species in the Qinghai-Tibet Plateau of China. Nevertheless, the molecular mechanisms underlying its adaptation to high-altitude environments remain largely unknown. In the present study, comparative transcriptome sequencing was performed for lung and gluteus tissues from two species of low-altitude cattle (Sanjiang and Holstein cattle), Tibetan cattle (living at a moderate altitude), and yak (living at a high altitude) and the differentially expressed genes were validated using real-time quantitative PCR. The results showed that CD36 antigen was up-regulated and CD59 antigen was down-regulated in yak in comparison to the other animals, which might promote the development of red blood cells and inhibit the development of lymphocytes in yak. In addition, thrombospondin type 1, coagulation factor 5/8, and fibronectin were all down-regulated, but serpin and alpha 2-macroglobulin (A2M) were up-regulated. These differences would inhibit blood coagulation, thus reducing the risk of pulmonary edema. The expression levels of the calcium-release, potassium, and transient receptor potential channels decreased in yak, minimizing membrane depolarization and the harmful effects of pulmonary edema. Eleven KEGG pathways associated with innate immunity were more activated in yak and Tibetan cattle than in other cattle strains, which should reduce their risk of infection and disease. These changes together might facilitate the adaptation of yak and Tibetan cattle to live in high-altitude habitats.

Cattleyak are the hybrid offspring between cattle and yak and combine yak hardiness with cattle productivity. Much attempt has been made to examine the mechanisms of male sterility caused by spermatogenic arrest, but yet there is no research systematically and precisely elucidated testis gene expression profiling between cattleyak and yak. To explore the higher resolution comparative transcriptome map between the testes of yak and cattleyak, and further analyze the mRNA expression dynamics of spermatogenic arrest in cattleyak. We characterized the comparative transcriptome profile from the testes of yak and cattleyak using high-throughput sequencing. Then we used quantitative analysis to validate several differentially expressed genes (DEGs) in testicular tissue and spermatogenic cells. Testis transcriptome profiling identified 6477 DEGs (2919 upregulated and 3558 downregulated) between cattleyak and yak. Further analysis revealed that the marker genes and apoptosis regulatory genes for undifferentiated spermatogonia were upregulated, while the genes for differentiation maintenance were downregulated in cattleyak. A majority of DEGs associated with mitotic checkpoint, and cell cycle progression were downregulated in cattleyak during spermatogonial mitosis. Furthermore, almost all DEGs related to synaptonemal complex assembly, and meiotic progression presented no sign of expression in cattleyak. Even worse, dozens of genes involved in acrosome formation, and flagellar development were dominantly downregulated in cattleyak. DEGs indicated that spermatogenic arrest of cattleyak may originate from the differentiation stage of spermatogonial stem cells and be aggravated during spermatogonial mitosis and spermatocyte meiosis, which contributes to the scarcely presented sperms in cattleyak.

Yak is an important livestock animal for the people indigenous to the harsh, oxygen‐limited Qinghai‐Tibetan Plateau and Hindu‐Kush ranges of the Himalayas. The yak genome was sequenced in 2012 but its assembly was fragmented because of the inherent limitations of the Illumina sequencing technology used to analyze it. An accurate and complete reference genome is essential for the study of genetic variations in this species. Long‐read sequences are more complete than their short‐read counterparts and have been successfully applied toward high‐quality genome assembly for various species. In this study, we present a high‐quality chromosome‐scale yak genome assembly (BosGru_PB_v1.0) constructed with long‐read sequencing and chromatin interaction technologies. Compared to an existing yak genome assembly (BosGru_v2.0), BosGru_PB_v1.0 shows substantially improved chromosome sequence continuity, reduced repetitive structure ambiguity, and gene model completeness. To characterize genetic variation in yak, we generated de novo genome assemblies based on Illumina short reads for seven recognized domestic yak breeds in Tibet and Sichuan and one wild yak from Hoh Xil. We compared these eight assemblies to the BosGru_PB_v1.0 genome, obtained a comprehensive map of yak genetic diversity at the whole‐genome level, and identified several protein coding genes absent from the BosGru_PB_v1.0 assembly. Despite the genetic bottleneck experienced by wild yak, their diversity was nonetheless higher than that of domestic yak. Here we identified breed‐specific sequences and genes by whole genome alignment, which may facilitate yak breed identification.

Skeletal muscle is a complex heterogeneous tissue and characterizing its cellular heterogeneity and transcriptional and epigenetic signatures are important for understanding the details of its ontogeny. In our study, we applied scRNA‐seq and scATAC‐seq to investigate the cell types, molecular features, transcriptional and epigenetic regulation, and patterns of developing bovine skeletal muscle from gestational, lactational and adult stages. Detailed molecular analyses were used to dissect cellular heterogeneity, and we deduced the differentiation trajectory of myogenic cells and uncovered their dynamic gene expression profiles. SCENIC analysis was performed to demonstrate key regulons during cell fate decisions. We explored the future expression states of these heterogeneous cells by RNA velocity analysis and found extensive networks of intercellular communication using the toolkit CellChat. Moreover, the transcriptomic and chromatin accessibility modalities were confirmed to be highly concordant, and integrative analysis of chromatin accessibility and gene expression revealed key transcriptional regulators acting during myogenesis. In bovine skeletal muscle, by scRNA‐seq and scATAC‐seq analysis, different cell types such as adipocytes, endothelial cells, fibroblasts, lymphocytes, monocytes, pericyte cells and eight skeletal myogenic subpopulations were identified at the three developmental stages. The pseudotime trajectory exhibited a distinct sequential ordering for these myogenic subpopulations and eight distinct gene clusters were observed according to their expression pattern. Moreover, specifically expressed TFs (such as MSC, MYF5, MYOD1, FOXP3, ESRRA, BACH1, SIX2 and ATF4) associated with muscle development were predicted, and likely future transcriptional states of individual cells and the developmental dynamics of differentiation among neighbouring cells were predicted. CellChat analysis on the scRNA‐seq data set then classified many ligand–receptor pairs among these cell clusters, which were further categorized into significant signalling pathways, including BMP, IGF, WNT, MSTN, ANGPTL, TGFB, TNF, VEGF and FGF. Finally, scRNA‐seq and scATAC‐seq results were successfully integrated to reveal a series of specifically expressed TFs that are likely to be candidates for the promotion of cell fate transition during bovine skeletal muscle development. Overall, our results outline a single‐cell dynamic chromatin/transcriptional landscape for normal bovine skeletal muscle development; these provide an important resource for understanding the structure and function of mammalian skeletal muscle, which will promote research into its biology.