Transcriptome profiles revealed the mechanisms underlying the adaptation of yak to high-altitude environments

Yak genome database: a multi-omics analysis platform

BACKGROUND

The yak (Bos grunniens) is a large ruminant species that lives in high-altitude regions and exhibits excellent adaptation to the plateau environments. To further understand the genetic characteristics and adaptive mechanisms of yak, we have developed a multi-omics database of yak including genome, transcriptome, proteome, and DNA methylation data.

DESCRIPTION

The Yak Genome Database ( http://yakgenomics.com/ ) integrates the research results of genome, transcriptome, proteome, and DNA methylation, and provides an integrated platform for researchers to share and exchange omics data. The database contains 26,518 genes, 62 transcriptomes, 144,309 proteome spectra, and 22,478 methylation sites of yak. The genome module provides access to yak genome sequences, gene annotations and variant information. The transcriptome module offers transcriptome data from various tissues of yak and cattle strains at different developmental stages. The proteome module presents protein profiles from diverse yak organs. Additionally, the DNA methylation module shows the DNA methylation information at each base of the whole genome. Functions of data downloading and browsing, functional gene exploration, and experimental practice were available for the database.

CONCLUSION

This comprehensive database provides a valuable resource for further investigations on development, molecular mechanisms underlying high-altitude adaptation, and molecular breeding of yak.

Genomic insights into Yak (Bos grunniens) adaptations for nutrient assimilation in high-altitudes

High-altitude environments present formidable challenges for survival and reproduction, with organisms facing limited oxygen availability and scarce nutrient resources. The yak (Bos grunniens), indigenous to the Tibetan Plateau, has notably adapted to these extreme conditions. This study delves into the genomic basis of the yak's adaptation, focusing on the positive selection acting on genes involved in nutrient assimilation pathways. Employing techniques in comparative genomics and molecular evolutionary analyses, we selected genes in the yak that show signs of positive selection associated with nutrient metabolism, absorption, and transport. Our findings reveal specific genetic adaptations related to nutrient metabolism in harsh climatic conditions. Notably, genes involved in energy metabolism, oxygen transport, and thermoregulation exhibited signs of positive selection, suggesting their crucial role in the yak's successful colonization of high-altitude regions. The study also sheds light on the yak's immune system adaptations, emphasizing genes involved in response to various stresses prevalent at elevated altitudes. Insights into the yak's genomic makeup provide valuable information for understanding the broader implications of high-altitude adaptations in mammalian evolution. They may contribute to efforts in enhancing livestock resilience to environmental challenges.

Single-cell RNA-seq and single-cell bisulfite-sequencing reveal insights into yak preimplantation embryogenesis

Extensive epigenetic reprogramming occurs during preimplantation embryonic development. However, the impact of DNA methylation in plateau yak preimplantation embryos and how epigenetic reprogramming contributes to transcriptional regulatory networks are unclear. In this study, we quantified gene expression and DNA methylation in oocytes and a series of yak embryos at different developmental stages and at single-cell resolution using single-cell bisulfite-sequencing and RNA-seq. We characterized embryonic genome activation and maternal transcript degradation and mapped epigenetic reprogramming events critical for embryonic development. Through cross-species transcriptome analysis, we identified 31 conserved maternal hub genes and 39 conserved zygotic hub genes, including SIN3A, PRC1, HDAC1/2, and HSPD1. Notably, by combining single-cell DNA methylation and transcriptome analysis, we identified 43 candidate methylation driver genes, such as AURKA, NUSAP1, CENPF, and PLK1, that may be associated with embryonic development. Finally, using functional approaches, we further determined that the epigenetic modifications associated with the histone deacetylases HDAC1/2 are essential for embryonic development and that the deubiquitinating enzyme USP7 may affect embryonic development by regulating DNA methylation. Our data represent an extensive resource on the transcriptional dynamics of yak embryonic development and DNA methylation remodeling, and provide new insights into strategies for the conservation of germplasm resources, as well as a better understanding of mammalian early embryonic development that can be applied to investigate the causes of early developmental disorders.

Study of Transcriptomic Analysis of Yak (Bos grunniens) and Cattle (Bos taurus) Pulmonary Artery Smooth Muscle Cells under Oxygen Concentration Gradients and Differences in Their Lung Histology and Expression of Pyruvate Dehydrogenase Kinase 1-Related Factors

The aim of this study was to investigate the molecular mechanisms by which hypoxia affects the biological behavior of yak PASMCs, the changes in the histological structure of yak and cattle lungs, and the relationships and regulatory roles that exist regarding the differences in the distribution and expression of PDK1 and its hypoxia-associated factors screened for their role in the adaptation of yak lungs to the plateau hypoxic environment. The results showed that, at the level of transcriptome sequencing, the molecular regulatory mechanisms of the HIF-1 signaling pathway, glucose metabolism pathway, and related factors (HK2/PGK1/ENO1/ENO3/ALDOC/ALDOA) may be closely related to the adaptation of yaks to the hypoxic environment of the plateau; at the tissue level, the presence of filled alveoli and semi-filled alveoli, thicker alveolar septa and basement membranes, a large number of erythrocytes, capillary distribution, and collagen fibers accounted for all levels of fine bronchioles in the lungs of yaks as compared to cattle. A higher percentage of goblet cells was found in the fine bronchioles of yaks, and PDK1, HIF-1α, and VEGF were predominantly distributed and expressed in the monolayers of ciliated columnar epithelium in the branches of the terminal fine bronchioles of yak and cattle lungs, with a small amount of it distributed in the alveolar septa; at the molecular level, the differences in PDK1 mRNA relative expression in the lungs of adult yaks and cattle were not significant (p > 0.05), the differences in HIF-1α and VEGF mRNA relative expression were significant (p < 0.05), and the expression of PDK1 and HIF-1α proteins in adult yaks was stronger than that in adult cattle. PDK1 and HIF-1α proteins were more strongly expressed in adult yaks than in adult cattle, and the difference was highly significant (p < 0.01); the relative expression of VEGF proteins was not significantly different between adult yaks and cattle (p > 0.05). The possible regulatory relationship between the above results and the adaptation of yak lungs to the plateau hypoxic environment paves the way for the regulatory mechanisms of PDK1, HIF-1α, and VEGF, and provides basic information for studying the mechanism of hypoxic adaptation of yaks in the plateau. At the same time, it provides a reference for human hypoxia adaptation and a target for the prevention and treatment of plateau diseases in humans and plateau animals.

Transcriptomic analysis reveals differentially expressed genes and key immune pathways in the spleen of the yak (Bos grunniens) at different growth stage

Yak is one of the rare and unique cattle species on the Qinghai-Tibetan Plateau, which has strong adaptability to the extreme environment of the plateau. The spleens are important functional organs that enable animals to adapt to their external environment and are vital in the growth and development process. To further investigate changes in immune function during yak development, we compared the transcriptome profiles of spleen tissues among juvenile (1-day old), youth (15-months old), and prime (5-years old) yaks. Immunology of spleen development was evaluated based on histological analyses and global gene expression was examined by using RNA-sequencing (RNA-seq) technology. In this work, we found 6378 genes with significant differences between the spleen of juvenile yak and youth yak, with the largest difference between groups. There were 3144 genes with significant differences between the spleen of young yak and prime yak, with the smallest differences between groups. Further, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted for the functional annotation of these genes. GO and KEGG analysis showed that some of them were related to growth, disease, immune, and metabolism. However, the genetic mechanism underlying the adaptability of yak spleens at different ages to harsh plateau environments remains unknown. These findings are important for studying the mechanisms of spleen development in yaks of different age groups.

Genome-wide comparison of DNA methylation patterns between yak and three cattle strains and their potential association with mRNA transcription

Yak has evolved specific adaptative mechanisms to high-altitude environment. Up to date, only a few studies reported the DNA methylation in yak. In the present study, genome-wide DNA methylome and transcriptome profiles in lung, mammary, and biceps brachii muscle tissues were compared between yak and three cattle breeds (Tibetan cattle, Sanjiang cattle, and Holstein cattle). The association between differentially expressed genes (DEGs) and differentially methylated regions (DMRs) was analyzed, and the biological functions of DEGs potentially driven by DMRs were explored by KEGG enrichment analysis. Finally, we found that yak-specific DMRs-driven DEGs were mainly involved in neuromodulation, respiration, lung development, blood pressure regulation, cardiovascular protection, energy metabolism, DNA repair, and immune functions. The higher levels of the key genes associated with these functions were observed in yak than in cattle, suggesting that DNA methylation might regulate these genes. Overall, the present study contributes basic data at the DNA methylation level to further understand the physiological metabolism in yak.

A de novo assembled genome of the Tibetan Partridge (Perdix hodgsoniae) and its high-altitude adaptation

The Tibetan Partridge (Perdix hodgsoniae) is an endemic species distributed in high-altitude areas of 3600-5600 m on the Qinghai-Tibet Plateau. To explore how the species is adapted to the high elevation environment, we assembled a draft genome based on both the Illumina and PacBio sequencing platforms with its population genetics and genomics analysis. In total, 134.74 Gb short reads and 30.81 Gb long reads raw data were generated. The 1.05-Gb assembled genome had a contig N50 of 4.56 Mb, with 91.94% complete BUSCOs. The 17 457 genes were annotated, and 11.35% of the genome was composed of repeat sequences. The phylogenetic tree showed that P. hodgsoniae was located at the basal position of the clade, including Golden Pheasant (Chrysolophus pictus), Common Pheasant (Phasianus colchicus), and Mikado Pheasant (Syrmaticus mikado). We found that 1014, 2595, and 2732 of the 6641 one-to-one orthologous genes were under positive selection in P. hodgsoniae, detected using PAML, BUSTED, and aBSREL programs, respectively, of which 965 genes were common under positive selection with 3 different programs. Several positively selected genes and immunity pathways relevant to high-altitude adaptation were detected. Gene family evolution showed that 99 gene families experienced significant expansion events, while 6 gene families were under contraction. The total number of olfactory receptor genes was relatively low in P. hodgsoniae. Genomic data provide an important resource for a further study on the evolutionary history of P. hodgsoniae, which provides a new insight into its high-altitude adaptation mechanisms.

Yak rumen microbiome elevates fiber degradation ability and alters rumen fermentation pattern to increase feed efficiency

Rumen microbes play an important role in ruminant energy supply and animal performance. Previous studies showed that yak (Bos grunniens) rumen microbiome and fermentation differ from other ruminants. However, little is understood about the features of the rumen microbiome that make yak adapted to their unique environmental and dietary conditions. This study was to investigate the rumen microbiome and metabolome to understand how yak adapt to the coarse forage and harsh environment in the Qinghai-Tibetan plateau. Nine female Qaidam yellow cattle (Bos taurus), 9 dzomo (hybrids of cattle and yak) and 9 female plateau yak (B. grunniens), about 5 to 6 years old, were used in this study. Rumen fermentation parameters, fibrolytic enzyme activities, and rumen metataxonomic were determined. Then 18 (6 samples per group) were selected for rumen metagenomic and metabolome analysis. Metataxonomic analysis revealed that the rumen microbiota was significantly different among plateau yak, Qaidam yellow cattle, and dzomo (P < 0.05). Metagenomic analysis displayed a larger gene pool encoding a richer repertoire of carbohydrate-active enzymes in the rumen microbiome of plateau yak and dzomo than Qaidam yellow cattle (P < 0.05). Some of the genes encoding glycoside hydrolases that mediate the digestion of cellulose and hemicellulose were significantly enriched in the rumen of plateau yak than Qaidam yellow cattle, but glycoside hydrolase 57 that primarily includes amylases was abundant in Qaidam yellow cattle (P < 0.05). The rumen fermentation profile differed also, Qaidam yellow cattle having a higher molar proportion of acetate but a lower molar proportion of propionate than dzomo and plateau yak (P < 0.05). Based on metabolomic analysis, rumen microbial metabolic pathways and metabolites were different. Differential metabolites are mainly amino acids, carboxylic acids, sugars, and bile acids. Changes in rumen microbial composition could explain the above results. The present study showed that the rumen microbiome of plateau yak helps its host to adapt to the Qinghai-Tibetan plateau. In particular, the plateau yak rumen microbiome has more enzymes genes involved in cellulase and hemicellulase than that of cattle, resulting higher fibrolytic enzyme activities in yak, further providing stronger fiber degradation function.

Genomic Diversity and Selection Signatures for Weining Cattle on the Border of Yunnan-Guizhou

Weining cattle is a Chinese indigenous breed influenced by complex breeding and geographical background. The multi-ethnic breeding culture makes Weining cattle require more attention as livestock resources for its genetic diversity. Here, we used 10 Weining cattle (five newly sequenced and five downloaded) and downloaded another 48 genome data to understand the aspects of Weining cattle: genetic diversity, population structure, and cold-adapted performance. In the current study, a high level of genetic diversity was found in Weining cattle, and its breed comprised two potential ancestries, which were Bos taurus and Bos indicus. The positive selective sweep analysis in Weining cattle was analyzed using composite likelihood ratio (CLR) and nucleotide diversity (θπ), resulting in 203 overlapped genes. In addition, we studied the cold adaptation of Weining cattle by comparing with other Chinese cattle (Wannan and Wenshan cattle) by three methods (FST, θπ-ratio, and XP-EHH). Of the top 1% gene list, UBE3D and ZNF668 were analyzed, and these genes may be associated with fat metabolism and blood pressure regulation in cold adaptation. Our findings have provided invaluable information for the development and conservation of cattle genetic resources, especially in southwest China.

Physiology and Transcriptomics Analysis Reveal the Contribution of Lungs on High-Altitude Hypoxia Adaptation in Tibetan Sheep

The Tibetan sheep is an indigenous species on the Tibetan plateau with excellent adaptability to high-altitude hypoxia and is distributed at altitudes of 2500–5000 m. The high-altitude hypoxia adaptation of Tibetan sheep requires adaptive reshaping of multiple tissues and organs, especially the lungs. To reveal the mechanisms of adaptation at the tissue and molecular levels in the lungs of Tibetan sheep under hypoxic conditions at different altitudes, we performed light and electron microscopic observations, transcriptomic sequencing, and enzyme-linked immunosorbent assay studies on the lungs of Tibetan sheep from three altitudes (2500, 3500, and 4500 m). The results showed that in addition to continuous increase in pulmonary artery volume, thickness, and elastic fiber content with altitude, Tibetan sheep increase the hemoglobin concentration at an altitude of 3500 m, while they decrease the Hb concentration and increase the surface area of gas exchange and capacity of the blood at an altitude of 4500 m. Other than that, some important differentially expressed genes related to angiogenesis (FNDC1, HPSE, and E2F8), vasomotion and fibrogenesis (GJA4, FAP, COL1A1, COL1A2, COL3A1, and COL14A1), and gas transport (HBB, HBA1, APOLD1, and CHL1) were also identified; these discoveries at the molecular level explain to some extent the physiological findings. In conclusion, the lungs of Tibetan sheep adopt different strategies when adapting to different altitudes, and these findings are valuable for understanding the basis of survival of indigenous species on the Tibetan plateau.

Comparative Transcriptome Analyses of Gayal (Bos frontalis), Yak (Bos grunniens), and Cattle (Bos taurus) Reveal the High-Altitude Adaptation

Gayal and yak are well adapted to their local high-altitude environments, yet the transcriptional regulation difference of the plateau environment among them remains obscure. Herein, cross-tissue and cross-species comparative transcriptome analyses were performed for the six hypoxia-sensitive tissues from gayal, yak, and cattle. Gene expression profiles for all single-copy orthologous genes showed tissue-specific expression patterns. By differential expression analysis, we identified 3,020 and 1,995 differentially expressed genes (DEGs) in at least one tissue of gayal vs. cattle and yak vs. cattle, respectively. Notably, we found that the adaptability of the gayal to the alpine canyon environment is highly similar to the yak living in the Qinghai-Tibet Plateau, such as promoting red blood cell development, angiogenesis, reducing blood coagulation, immune system activation, and energy metabolism shifts from fatty acid β-oxidation to glycolysis. By further analyzing the common and unique DEGs in the six tissues, we also found that numerous expressed regulatory genes related to these functions are unique in the gayal and yak, which may play important roles in adapting to the corresponding high-altitude environment. Combined with WGCNA analysis, we found that UQCRC1 and COX5A are the shared differentially expressed hub genes related to the energy supply of myocardial contraction in the heart-related modules of gayal and yak, and CAPS is a shared differential hub gene among the hub genes of the lung-related module, which is related to pulmonary artery smooth muscle contraction. Additionally, EDN3 is the unique differentially expressed hub gene related to the tracheal epithelium and pulmonary vasoconstriction in the lung of gayal. CHRM2 is a unique differentially expressed hub gene that was identified in the heart of yak, which has an important role in the autonomous regulation of the heart. These results provide a basis for further understanding the complex transcriptome expression pattern and the regulatory mechanism of high-altitude domestication of gayal and yak.

Integrated proteomic, phosphoproteomic, and N-glycoproteomic analyses of the longissimus thoracis of yaks

Yaks (Bos mutus) live in the Qinghai–Tibet plateau. The quality of yak meat is unique due to its genetic and physiological characteristics. Identification of the proteome of yak muscle could help to reveal its meat-quality properties. The common proteome, phosphoproteome, and N-glycoproteome of yak longissimus thoracis (YLT) were analyzed by liquid chromatography-tandem mass spectrometry-based shotgun analysis. A total of 1812 common proteins, 1303 phosphoproteins (3918 phosphorylation sites), and 204 N-glycoproteins (285 N-glycosylation sites) were identified in YLT. The common proteins in YLT were involved mainly in myofibril structure and energy metabolism; phosphoproteins were associated primarily with myofibril organization, regulation of energy metabolism, and signaling; N-glycoproteins were engaged mainly in extracellular-matrix organization, cellular immunity, and organismal homeostasis. We reported, for the first time, the “panorama” of the YLT proteome, specifically the N-glycoproteome of YLT. Our results provide essential information for understanding post mortem physiology (rigor mortis and aging) and the quality of yak meat.

•

A total of 2650 proteins were identified in yak longissimus thoracis.

•

Common proteins were involved mainly in myofibril structure and energy metabolism.

•

Phosphoproteins were associated with myofibrils, energy metabolism, and signaling.

•

N-glycoproteins were engaged mainly in ECM organization, immunity, and homeostasis.

Metabolomics reveals potential plateau adaptability by regulating inflammatory response and oxidative stress-related metabolism and energy metabolism pathways in yak

Species are facing strong selection pressures to adapt to inhospitable high-altitude environments. Yaks are a valuable species and an iconic symbol of the Qinghai-Tibet Plateau. Extensive studies of high-altitude adaptation have been conducted, but few have focused on metabolism. In the present study, we determined the differences in the serum metabolomics between yaks and the closely related species of low-altitude yellow cattle and dairy cows. We generated high-quality metabolite profiling data for 36 samples derived from the three species, and a clear separation trend was obtained between yaks and the other animals from principal component analysis. In addition, we identified a total of 63 differentially expressed metabolites among the three species. Functional analysis revealed that differentially expressed metabolites were related to the innate immune activation, oxidative stress-related metabolism, and energy metabolism in yaks, which indicates the important roles of metabolites in high-altitude adaptation in yaks. The results provide new insights into the mechanism of adaptation or acclimatization to high-altitude environments in yaks and hypoxia-related diseases in humans.

Differences in proteomic profiles between yak and three cattle strains provide insights into molecular mechanisms underlying high-altitude adaptation

Yaks display unique properties of the lung and heart, enabling their adaptation to high-altitude environments, but the underlying molecular mechanisms are still largely unknown. In the present study, the proteome differences in lung and heart tissues were compared between yak (Bos grunniens) and three cattle strains (Bos taurus, Holstein, Sanjiang and Tibetan cattle) using the sequential window acquisition of all theoretical mass spectra/data-independent acquisition (SWATH/DIA) proteomic method. In total, 51,755 peptides and 7215 proteins were identified. In the lung tissue, there were 162, 310 and 118 differential abundance proteins (DAPs) in Tibetan, Holstein and Sanjiang cattle compared to yak respectively. In the heart tissue, there were 71, 57 and 78 DAPs in Tibetan, Holstein and Sanjiang cattle compared to yak respectively. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the DAPs were enriched for the retinol metabolism and toll-like receptor categories in lung tissue. The changes in these two pathways may regulate hypoxia-induced factor and immune function in yaks. Moreover, DAPs in heart tissues were enriched for cardiac muscle contraction, Huntington's disease, chemical carcinogenesis and drug metabolism-cytochrome P450. Further exploration indicated that yaks may alter cardiac function through regulation of type 2 ryanodine receptor (RyR2) and Ca²⁺ -release channels. The present results are useful to further develop an understanding of the mechanisms underlying adaptation of animals to high-altitude conditions.

Adaptation Mechanisms of Yak (Bos grunniens) to High-Altitude Environmental Stress

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.

Changes in the expression levels of elastic fibres in yak lungs at different growth stages

BACKGROUND

Yaks have a strong adaptability to the plateau environment, which can be attributed to the effective oxygen utilization rate of their lung tissue. Elastic fibre confers an important adaptive structure to the alveolar tissues in yaks. However, little research has been focused on the structural development of lung tissues and the expression levels of elastic fibres in yaks after birth. Therefore, this study aimed to investigate the morphological changes of elastic fibers and expression profiles of fibre-formation genes in yak lungs at different growth stages and the relationship between these changes and plateau adaptation.

RESULTS

Histological staining was employed to observe the morphological changes in the lung tissue structure of yaks at four different ages: 1 day old, 30 days old, 180 days old and adult. There was no significant difference in the area of a single alveolus between the 1-day-old and 30-day-old groups (P-value > 0.05). However, the single alveolar area was gradually increased with an increase in age (P-value < 0.05). Elastic fibre staining revealed that the amount of elastic fibres in alveolar tissue was increased significantly from the ages of 30 days to 180 days (P-value < 0.05) and stabilized during the adult stage. Transcriptome analysis indicated that the highest levels of differentially expressed genes were found between 30 days of age and 180 days of age. KEGG analysis showed that PI3K-Akt signalling pathway and MAPK pathway, which are involved in fibre formation, accounted for the largest proportion of differentially expressed genes between 30 days of age and 180 days of age. The expression levels of 36 genes related to elastic fibre formation and collagen fibre formation were also analysed, and most of these genes were highly expressed in 30-day-old and 180-day-old yaks.

CONCLUSIONS

The content of elastic fibres in the alveolar tissue of yaks increases significantly after birth, but this change occurs only from 30 days of age to 180 days of age. Our study indicates that elastic fibres can improve the efficiency of oxygen utilization in yaks under harsh environmental conditions.

Environmental gradients of selection for an alpine-obligate bird, the white-tailed ptarmigan (Lagopus leucura)

The warming climate will expose alpine species adapted to a highly seasonal, harsh environment to novel environmental conditions. A species can shift their distribution, acclimate, or adapt in response to a new climate. Alpine species have little suitable habitat to shift their distribution, and the limits of acclimation will likely be tested by climate change in the long-term. Adaptive genetic variation may provide the raw material for species to adapt to this changing environment. Here, we use a genomic approach to describe adaptive divergence in an alpine-obligate species, the white-tailed ptarmigan (Lagopus leucura), a species distributed from Alaska to New Mexico, across an environmentally variable geographic range. Previous work has identified genetic structure and morphological, behavioral, and physiological differences across the species' range; however, those studies were unable to determine the degree to which adaptive divergence is correlated with local variation in environmental conditions. We used a genome-wide dataset generated from 95 white-tailed ptarmigan distributed throughout the species' range and genotype-environment association analyses to identify the genetic signature and environmental drivers of local adaptation. We detected associations between multiple environmental gradients and candidate adaptive loci, suggesting ptarmigan populations may be locally adapted to the plant community composition, elevation, local climate, and to the seasonality of the environment. Overall, our results suggest there may be groups within the species' range with genetic variation that could be essential for adapting to a changing climate and helpful in guiding conservation action.

Evolution and domestication of the Bovini species

Domestication of the Bovini species (taurine cattle, zebu, yak, river buffalo and swamp buffalo) since the early Holocene (ca. 10 000 BCE) has contributed significantly to the development of human civilization. In this study, we review recent literature on the origin and phylogeny, domestication and dispersal of the three major Bos species - taurine cattle, zebu and yak - and their genetic interactions. The global dispersion of taurine and zebu cattle was accompanied by population bottlenecks, which resulted in a marked phylogeographic differentiation of the mitochondrial and Y-chromosomal DNA. The high diversity of European breeds has been shaped through isolation-by-distance, different production objectives, breed formation and the expansion of popular breeds. The overlapping and broad ranges of taurine and zebu cattle led to hybridization with each other and with other bovine species. For instance, Chinese gayal carries zebu mitochondrial DNA; several Indonesian zebu descend from zebu bull × banteng cow crossings; Tibetan cattle and yak have exchanged gene variants; and about 5% of the American bison contain taurine mtDNA. Analysis at the genomic level indicates that introgression may have played a role in environmental adaptation.

Transcriptome analysis identified long non-coding RNAs involved in the adaption of yak to high-altitude environments

The mechanisms underlying yak adaptation to high-altitude environments have been investigated using various methods, but no report has focused on long non-coding RNA (lncRNA). In the present study, lncRNAs were screened from the gluteus transcriptomes of yak and their transcriptional levels were compared with those in Sanjiang cattle, Holstein cattle and Tibetan cattle. The potential target genes of the differentially expressed lncRNAs between species/strains were predicted using cis and trans models. Based on cis-regulated target genes, no KEGG pathway was significantly enriched. Based on trans-regulated target genes, 11 KEGG pathways in relation to energy metabolism and three KEGG pathways associated with muscle contraction were significantly enriched. Compared with cattle strains, transcriptional levels of acyl-CoA dehydrogenase, acyl-CoA-binding protein, 3-hydroxyacyl-CoA dehydrogenase were relatively higher and those of glyceraldehyde 3-phosphate dehydrogenase, phosphoglycerate mutase 1, pyruvate kinase and lactate/malate dehydrogenase were relatively lower in yak, suggesting that yaks activated fatty acid oxidation but inhibited glucose oxidation and glycolysis. Besides, NADH dehydrogenase and ATP synthase showed lower transcriptional levels in yak than in cattle, which might protect muscle tissues from deterioration caused by reactive oxygen species (ROS). Compared with cattle strains, the higher transcriptional level of glyoxalase in yak might contribute to dicarbonyl stress resistance. Voltage-dependent calcium channel/calcium release channel showed a lower level in yak than in cattle strains, which could reduce the Ca²⁺ influx and subsequently decrease the risk of hypertension. However, levels of EF-hand and myosin were higher in yak than in cattle strains, which might enhance the negative effects of reduced Ca²⁺ on muscle contraction. Overall, the present study identified lncRNAs and proposed their potential regulatory functions in yak.

Signature of high altitude adaptation in the gluteus proteome of the yak

Yak is the unique Bovidae species in the Qinghai-Tibetan Plateau. A previous proteomic study has compared the yak muscle tissue to one cattle strain using the isobaric tags for relative and absolute quantification approach. In this study, to further investigate the molecular mechanisms underlying yak adaptation, the proteomic profiles of gluteus were compared between yak and one moderate-altitude cattle strain (Tibetan cattle) and two low-altitude cattle strains (Holstein and Sanjiang cattle) using a label-free quantitative method. The comparisons identified 20, 364, 143 upregulated proteins and 4, 6, 37 downregulated proteins in yak, compared with Tibetan, Holstein, and Sanjiang cattle, respectively. Protein-protein interaction analysis indicated that these differentially expressed proteins were mainly related to "oxidative phosphorylation" and "electron transport chain." Further analysis revealed that NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 11, NADH dehydrogenase [ubiquinone] 1 beta subcomplex subunit 4, cytochrome C oxidase subunit 6A2, mitochondrial and cytochrome c oxidase subunit NDUFA4 were all increased in the yak, suggesting that yak might increase mitochondrial capacity to sustain metabolic rates under high altitude conditions, which might be a long-term adaptive mechanism underlying adaptation to high altitude environments. Yak increased the level of thioredoxin reductase 2 to protect themselves from oxidative damages. Moreover, the increased expression levels of phosphatidylinositol 4,5-bisphosphate 3-kinase catalytic subunit alpha isoform and caveolin-1 in yak suggested that yaks promoted glucose uptake for adaptation to high altitude. These results provided more information to better understand the molecular mechanisms underlying yak adaption.

Insights into the genetic characteristics and population structures of Chinese two Tibetan groups using 35 insertion/deletion polymorphic loci

Studying the genetic structure of each ethnic group is helpful to clarify the genetic background and trace back to the ethnic origin. Tibetan people have lived in the Qinghai-Tibet Plateau (mean elevation over 4500 m) for generations, and have well adapted to the high-altitude environment. Due to the relatively closed geographical environment, Tibetans have preserved their representative physical characteristics and genetic information, thereby become an important research group in human genetics. In this study, genetic characteristics and population structures of two Tibetan groups (Qinghai Tibetans and Tibet Tibetans) were revealed by 35 insertion/deletion polymorphism (DIP) loci, aiming to provide valuable genetic information for population genetic differentiation analyses and forensic identifications. The combined discrimination power, cumulative exclusion probability and combined match probability of the 35 DIP loci in Qinghai Tibetan and Tibet Tibetan groups were 0.9999999999999945, 0.9988, 5.56623 × 10^-15; and 0.9999999999999904, 0.9990, 9.69071 × 10^-15, respectively, indicating that the panel possessed a strong capability for Tibetan personal identifications. Population differentiations and genetic relationship analyses among the two studied Tibetan groups and other 27 comparison populations were carried out using the Nei's D_A genetic distances, population pairwise genetic distances F-statistics (F_ST), analysis of molecular variance (AMOVA), phylogenetic tree reconstruction, principal component analysis and STRUCTURE methods. Results demonstrated that the most intimate genetic relationships existed in these two Tibetan groups; and genetic similarities between two Tibetan groups and the populations from East Asia were much stronger than that between the Tibetan groups and other geographical populations. Furthermore, forensic ancestral informativeness assessments suggested that several loci could be regarded as ancestry informative markers inferring individual biogeographic origins as well as contributing to forensic anthropology and population genetic researches.

Reference gene selection and myosin heavy chain (MyHC) isoform expression in muscle tissues of domestic yak (Bos grunniens)

Domestic yak (Bos grunniens) is the most crucial livestock in the Qinghai-Tibetan Plateau, providing meat and other necessities for local people. The skeletal muscle of adult livestock is composed of muscle fibers, and fiber composition in muscle has influence on meat qualities, such as tenderness, pH, and color. Real-time quantitative polymerase chain reaction (RT-qPCR) is a powerful tool to evaluate the gene expression of muscle fiber, but the normalization of the data depends on the stability of expressed reference genes. Unfortunately, there is no consensus for an ideal reference gene for data normalization in muscle tissues of yak. In this study, we aimed to assess the stability of 14 commonly used candidate reference genes by using five algorithms (GeNorm, NormFinder, BestKeeper, Delat Ct and Refinder). Our results suggested UXT and PRL13A were the most stable reference genes, while the most commonly used reference gene, GAPDH, was most variably expressed across different muscle tissues. We also found that the extensor digitorum lateralis (EDL), trapezius pars thoracica (TPT), and psoas major (PM) muscle had the higher content of type I muscle fibers and the lowest content of type IIB muscle fibers, while gluteobiceps (GB) muscle had the highest content of type IIB muscle fibers. Our study provides the suitable reference genes for accurate analysis of yak muscle fiber composition.

Selection signatures for high-altitude adaptation in ruminants

High-altitude areas are important socio-economical habitats with ruminants serving as a major source of food and commodities for humans. Living at high altitude, however, is extremely challenging, predominantly due to the exposure to hypoxic conditions, but also because of cold temperatures and limited feed for livestock. To survive in high-altitude environments over the long term, ruminants have evolved adaptation strategies, e.g. physiological and morphological modifications, which allow them to cope with these harsh conditions. Identification of such selection signatures in ruminants may contribute to more informed breeding decisions, and thus improved productivity. Moreover, studying the genetic background of altitude adaptation in ruminants provides insights into a common molecular basis across species and thus a better understanding of the physiological basis of this adaptation. In this paper, we review the major effects of high altitude on the mammalian body and highlight some of the most important short-term (coping) and genetically evolved (adaptation) physiological modifications. We then discuss the genetic architecture of altitude adaptation and target genes that show evidence of being under selection based on recent studies in various species, with a focus on ruminants. The yak is presented as an interesting native species that has adapted to the high-altitude regions of Tibet. Finally, we conclude with implications and challenges of selection signature studies on altitude adaptation in general. We found that the number of studies on genetic mechanisms that enable altitude adaptation in ruminants is growing, with a strong focus on identifying selection signatures, and hypothesise that the investigation of genetic data from multiple species and regions will contribute greatly to the understanding of the genetic basis of altitude adaptation.

The Selection of Reference Genes for Quantitative Real-Time PCR in the Ashidan Yak Mammary Gland During Lactation and Dry Period

Simple Summary

The Ashidan yak is a new cultivated breed which has polled characteristics and a mild temperament. Improving milk yield is an important aspect of a breeding program for this breed. The mammary gland undergoes dramatic physiological and metabolic changes during the transition from lactation to dry periods, which involves the expression and regulation of a great number of genes. Quantification of gene expression levels by real-time quantitative polymerase chain reaction (RT-qPCR) is important to reveal the molecular mechanisms of mammary gland development and lactation. The accuracy of RT-qPCR is strongly influenced by the expression stability of reference genes, however, a systematic approach for selecting reference genes used for analyzing gene expression of the Ashidan yak has not been developed. In this study, we selected reference genes and analyzed their expression stability at different physiological stages (lactation and dry period). We found the hydroxymethylbilane synthase gene (HMBS) and the tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide gene (YWHAZ) were the most stable genes of the mammary gland of the Ashidan yak. These results help to improve the accuracy of gene expression analysis and provide a basis for future functional studies of target gene expression in the mammary gland of the Ashidan yak.

Abstract

Investigating the critical genes related to milk synthesis is essential for the improvement of the milk yield of the yak. Real-time quantitative polymerase chain reaction (RT-qPCR) is a reliable and widely used method to measure and evaluate gene expression levels. Selection of suitable reference genes is mandatory to acquire accurate normalization of gene expression results from RT-qPCR. To select the most stable reference genes for reliable normalization of mRNA expression by RT-qPCR in the mammary gland of the Ashidan yak, we selected 16 candidate reference genes and analyzed their expression stability at different physiological stages (lactation and dry period). The expression stability of the candidate reference genes was assessed using geNorm, NormFinder, BestKeeper, Delta Ct, and RefFinder methods. The results showed that the hydroxymethylbilane synthase gene (HMBS) and the tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide gene (YWHAZ) were the most stable genes across all treatment samples. The reliability of selected reference genes was validated by normalizing relative expression of the lactation-related 60S ribosomal protein L35 gene (RPL35). The relative expression of RPL35 varied considerably according to the different reference genes. This work provides valuable information to further promote research in the molecular mechanisms involved in lactation and mammary gland development and provides a foundation for the improvement of the milk yield and quality of the Ashidan yak.

Comparisons of lung and gluteus transcriptome profiles between yaks at different ages

The yak, Bos grunniens, is the only large mammal in the Qinghai-Tibet Plateau and has been bred to provide meat, milk, and transportation. Previous studies indicate that the immune system contributes to the yak's adaptation to high-altitude environments. In order to further investigate changes in immune function during yak development, we compared the transcriptome profiles of gluteus and lung tissues among yaks at 6, 30, 60, and 90 months of age. Analyses of significantly differentially expressed genes (DEGs) in lung tissues revealed that immune function was more activated at 6-months and less activated at 90-months than in the 30 and 60-month-old animals. DEG exploration in gluteal tissues revealed that immune functions were more highly activated at both 6 and 90-months, compared with 30 and 60-months. Immune system activation in the muscle and lung tissues of 30-month-old yaks may increase their resistance to infections, while decreased may be due to aging. Furthermore, the higher immune activation status in the gluteal tissues in 90-month-old yaks could be due to muscle injury and subsequent regeneration, which is supported by the fact that 5 unigenes related with muscle injury and 3 related to muscle regeneration displayed greater expression levels at 90-months than at 30 and 60-months. Overall, the present study highlights the important role of the immune system in yak development, which will facilitate future investigations.