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.

Background: The Intramuscular fat (IMF) content in meat products, which is positively correlated with meat quality, is an important trait considered by consumers. The regulation of IMF deposition is species specific. However, the IMF-deposition-related mRNA and non-coding RNA and their regulatory network in yak (Bos grunniens) remain unknown. High-throughput sequencing technology provides a powerful approach for analyzing the association between transcriptome-related differences and specific traits in animals. Thus, the whole transcriptomes of yak muscle and adipose tissues were screened and analyzed to elucidate the IMF deposition-related genes. The muscle tissues were used for IMF content measurements. Results: Significant differences were observed between the 0.5- and 2.5-year-old yaks. Several mRNAs, miRNAs, lncRNAs and circRNAs were generally expressed in both muscle and adipose tissues. Between the 0.5- and 2.5-year-old yaks, 149 mRNAs, 62 miRNAs, 4 lncRNAs, and 223 circRNAs were differentially expressed in muscle tissue, and 72 mRNAs, 15 miRNAs, 9 lncRNAs, and 211 circRNAs were differentially expressed in adipose tissue. KEGG annotation revelved that these differentially expressed genes were related to pathways that maintain normal biological functions of muscle and adipose tissues. Moreover, 16 mRNAs, 5 miRNAs, 3 lncRNAs, and 5 circRNAs were co-differentially expressed in both types of tissue. We suspected that these co-differentially expressed genes were involved in IMF-deposition in the yak. Additionally, LPL, ACADL, SCD, and FASN, which were previously shown to be associated with the IMF content, were identified in the competing endogenous RNA (ceRNA) regulatory network that was constructed on the basis of the IMF deposition-related genes. Three ceRNA subnetworks also revealed that TCONS-00016416 and its target SIRT1 “talk” to each other through the same miR-381-y and miR-208 response elements, whereas TCONS-00061798 and its target PRKCA, and TCONS-00084092 and its target LPL “talk” to each other through miR-122-x and miR-499-y response elements, respectively. Conclusion: Taken together, our results reveal the potential mRNA and noncoding RNAs involved in IMF deposition in the yak, providing a useful resource for further research on IMF deposition in this animal species.

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.

Background: On the Qinghai–Tibet Plateau, known as the roof ridge of the world, the yak is a precious cattle species that has been indispensable to the human beings living in this high-altitude area. However, the origin of domestication, dispersal route, and the divergence of domestic yaks from different areas are poorly understood. Results: Here, we resequenced the genome of 91 domestic yak individuals from 31 populations and 1 wild yaks throughout China. Using a population genomics approach, we observed considerable genetic variation. Phylogenetic analysis suggested that the earliest domestications of yak occurred in the south-eastern QTP, followed by dispersal to the west QTP and northeast to SiChuang, Gansu, and Qinghai by two routes. Interestingly, we also found potential associations between the distribution of some breeds and historical trade routes such as the Silk Road and Tang-Tibet Ancient Road. Selective analysis identified 11 genes showing differentiation between domesticated and wild yaks and the potentially positively selected genes in each group were identified and compared among domesticated groups. We also detected an unbalanced pattern of introgression among domestic yak, wild yak, and Tibetan cattle. Conclusions: Our research revealed population genetic evidence for three groups of domestic yaks. In addition to providing genomic evidence for the domestication history of yaks, we identified potential selected genes and introgression, which provide a theoretical basis and resources for the selective breeding of superior characters and high-quality yak.