Transcriptome profile and unique genetic evolution of positively selected genes in yak lungs

Effect of the TGF-β/BMP Signaling Pathway on the Proliferation of Yak Pulmonary Artery Smooth Muscle Cells under Hypoxic Conditions

To survive in low-oxygen environments, yaks effectively avoid hypoxia-induced pulmonary arterial hypertension through vascular remodeling. The TGF-β/BMP signaling pathway plays a key role in maintaining the homeostasis of pulmonary artery smooth muscle cells (PASMCs). However, little is known about the molecular regulatory mechanisms by which the TGF-β/BMP signaling pathway contributes to the proliferation of yak PASMCs. In this study, yak PASMCs were cultured in vitro, and a hypoxia model was constructed to investigate the effect of TGFβ/BMP signaling on yak PASMC proliferation. Hypoxia treatment increased the proliferation of yak PASMCs significantly. As the duration of hypoxia increased, the expression levels of TGF-β1 and the phosphorylation levels of Smad2/3 were upregulated significantly. The BMP signaling pathway was transiently activated by hypoxia, with increases in BMPR2 expression and Smad1/5 phosphorylation, and these changes were gradually reversed with prolonged hypoxia exposure. In addition, exogenous TGF-β1 activated the TGF-β signaling pathway, increased the phosphorylation levels of the downstream proteins Smad2 and Smad3, and increased the proliferation and migration rates of yak PASMCs significantly. Finally, treatment with noggin (an inhibitor of BMP signaling) significantly reduced BMPR2 protein expression levels and Smad1/5 phosphorylation levels and increased yak PASMC proliferation and migration rates. In summary, these results revealed that under hypoxic conditions, the dynamic regulation of the TGF-β/BMP signaling pathway promotes the proliferation of yak PASMCs.

Transcriptomics reveals age-related changes in ion transport-related factors in yak lungs

Yaks inhabit high-altitude, low-oxygen regions, where ion transport functions play a crucial role in maintaining intracellular and extracellular ionic balance and regulating pulmonary vascular tension. These functions affect pulmonary ventilation and blood flow rate, aiding tissue development and enhancing oxygen transfer efficiency, thus facilitating better adaptation to hypoxic environments. To investigate the regulatory mechanisms of ion transport-related factors on the growth and development of yak lungs, we employed RNA sequencing (RNA-seq)for sequencing the transcriptome in the lung tissues of neonatal (1-day-old), juvenile (1-year-old), and adult (4-year-old) yaks. We also performed differential gene expression and functional analyses. The results yielded 26 genes associated with ion transport, mainly enriched in the salivary and pancreatic secretion pathways. Finally, we used several methods including quantitative polymerase chain reaction (qRT-PCR), and Western blotting (WB), immunohistochemical (IHC) and immunofluorescence (IF) staining to determine the distribution of the expression of the ion transport genes FOXI1, KCNMA1, and SLC12A2 in yak lung tissues. qRT-PCR and WB results indicated that mRNA and protein relative expression levels of FOXI1 and SLC12A2 were significantly higher in neonatal yaks than in juvenile and adult yaks (all p < 0.05), whereas those of KCNMA1 were significantly higher in adult yaks than in neonatal and juvenile yaks (all p < 0.05). IHC and IF results demonstrated that FOXI1, KCNMA1, and SLC12A2 were distributed among the epithelial mucosal layers (including ciliated, goblet, and Clara cells) of the yaks' bronchi and their branches in the lungs across different age groups of yak. Therefore, our results suggested that FOXI1, KCNMA1, and SLC12A2 may be strongly associated with the development and aging processes in yak lungs. These results provide insights into the molecular mechanisms underlying the yak's adaptation to high-altitude environments and valuable references for further research.

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.

Combined transcriptome and proteome analysis of yak PASMCs under hypoxic and normoxic conditions

Background

Yaks are animals that have lived in plateau environments for generations. Yaks can adapt to the hypoxic plateau environment and also pass this adaptability on to the next generation. The lungs are the most important respiratory organs for mammals to adapt to their environment. Pulmonary artery smooth muscle cells play an important role in vascular remodeling under hypoxia, but the genetic mechanism underpinning the yak's ability to adapt to challenging plateau conditions is still unknown.

Methods

A tandem mass tag (TMT) proteomics study together with an RNA-seq transcriptome analysis were carried out on pulmonary artery smooth muscle cells (PASMCs) that had been grown for 72 hours in both normoxic (20% O2) and hypoxic (1% O2) environments. RNA and TP (total protein) were collected from the hypoxic and normoxic groups for RNA-seq transcriptome sequencing and TMT marker protein quantification, and RT-qPCR validation was performed.

Results

A total of 17,711 genes and 6,859 proteins were identified. Further, 5,969 differentially expressed genes (DEGs) and 531 differentially expressed proteins (DEPs) were identified in the comparison group, including 2,924 and 186 upregulated genes and proteins and 3,045 and 345 down-regulated genes and proteins, respectively. The transcriptomic and proteomic analyses revealed that 109 DEGs and DEPs were highly positively correlated, with 77 genes showing the same expression trend. Nine overlapping genes were identified in the HIF-1 signaling pathway, glycolysis / gluconeogenesis, central carbon metabolism in cancer, PPAR signaling pathway, AMPK signaling pathway, and cholesterol metabolism (PGAM1, PGK1, TPI1, HMOX1, IGF1R, OLR1, SCD, FABP4 and LDLR), suggesting that these differentially expressed genes and protein functional classifications are related to the hypoxia-adaptive pathways. Overall, our study offers abundant data for further analysis of the molecular mechanisms in yak PASMCs and their adaptability to different oxygen concentrations.

Transcriptome analysis of pika heart tissue reveals mechanisms underlying the adaptation of a keystone species on the roof of the world

High-altitude environments impose intense stresses on living organisms and drive striking phenotypic and genetic adaptations, such as hypoxia resistance, cold tolerance, and increases in metabolic capacity and body mass. As one of the most successful and dominant mammals on the Qinghai-Tibetan Plateau (QHTP), the plateau pika (Ochotona curzoniae) has adapted to the extreme environments of the highest altitudes of this region and exhibits tolerance to cold and hypoxia, in contrast to closely related species that inhabit the peripheral alpine bush or forests. To explore the potential genetic mechanisms underlying the adaptation of O. curzoniae to a high-altitude environment, we sequenced the heart tissue transcriptomes of adult plateau pikas (comparing specimens from sites at two different altitudes) and Gansu pikas (O. cansus). Differential expression analysis and weighted gene co-expression network analysis (WGCNA) were used to identify differentially expressed genes (DEGs) and their primary functions. Key genes and pathways related to high-altitude adaptation were identified. In addition to the biological processes of signal transduction, energy metabolism and material transport, the identified plateau pika genes were mainly enriched in biological pathways such as the negative regulation of smooth muscle cell proliferation, the apoptosis signalling pathway, the cellular response to DNA damage stimulus, and ossification involved in bone maturation and heart development. Our results showed that the plateau pika has adapted to the extreme environments of the QHTP via protection against cardiomyopathy, tissue structure alterations and improvements in the blood circulation system and energy metabolism. These adaptations shed light on how pikas thrive on the roof of the world.

Plateau Adaptation Gene Analyses Reveal Transcriptomic, Proteomic, and Dual Omics Expression in the Lung Tissues of Tibetan and Yorkshire Pigs

Elevated environments such as plateaus are often classified as low oxygen environments. The hypoxic adaptation mechanisms utilized by organisms in these conditions are not well understood. To address this, the differentially expressed genes (DEGs) involved in hypoxia adaptation were assessed using two pig breeds (Tibetan pig [TP] and Yorkshire sow [YY]). Genes related to lung tissue responses to hypoxia were assessed using transcriptomic (using RNA-seq) and proteomic (using iTRAQ) analysis. A total of 1021 DEGs were screened out. In the iTRAQ omics data, a total of 22,100 peptides were obtained and 4518 proteins were found after filtering. A total of 271 differentially expressed proteins [DEPs] were screened using the conditions of p < 0.05; FC ≤ 0.833; and FC ≥ 1.2. A total of 14 DEGs at the mRNA and protein levels were identified and found to be associated with regulation of the inflammatory response; blood particles; and MAPK cascade response regulation. Among the DEGs, six were associated with hypoxia adaptation function (mitochondria and glycolysis) in pigs. The results of this study identify novel candidate genes involved in porcine hypoxia adaptation mechanisms.

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.

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.

Characteristics of pulmonary microvascular structure in postnatal yaks

Yaks are typical plateau-adapted animals, however the microvascular changes and characteristics in their lungs after birth are still unclear. Pulmonary microvasculature characteristics and changes across age groups were analysed using morphological observation and molecular biology detection in yaks aged 1, 30 and 180 days old in addition to adults. Results: Our experiments demonstrated that yaks have fully developed pulmonary alveolar at birth but that interalveolar thickness increased with age. Immunofluorescence observations showed that microvessel density within the interalveolar septum in the yak gradually increased with age. In addition, transmission electron microscopy (TEM) results showed that the blood–air barrier of 1-day old and 30-days old yaks was significantly thicker than that observed at 180-days old and in adults (P < 0.05), which was caused by the thinning of the membrane of alveolar epithelial cells. Furthermore, Vegfa and Epas1 expression levels in 30-day old yaks were the highest in comparison to the other age groups (P < 0.05), whilst levels in adult yaks were the lowest (P < 0.05). The gradual increase in lung microvessel density can effectively satisfy the oxygen requirements of ageing yaks. In addition, these results suggest that the key period of yak lung development is from 30 to 180 days.

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.

Comparative Analysis of Skeleton Muscle Proteome Profile between Yak and Cattle Provides Insight into High-Altitude Adaptation

Background::

Mechanisms underlying yak adaptation to high-altitude environments have been investigated at the levels of morphology, anatomy, physiology, genome and transcriptome, but have not been explored at the proteome level.

Objective:

The protein profiles were compared between yak and cattle to explore molecular mechanisms underlying yak adaptation to high altitude conditions.

Methods:

In the present study, an antibody microarray chip was developed, which included 6,500 mouse monoclonal antibodies. Immunoprecipitation and mass spectrometry were performed on 12 selected antibodies which showed that the chip was highly specific. Using this chip, muscle tissue proteome was compared between yak and cattle, and 12 significantly Differentially Expressed Proteins (DEPs) between yak and cattle were identified. Their expression levels were validated using Western blot.

Results:

ompared with cattle, higher levels of Rieske Iron-Sulfur Protein (RISP), Cytochrome C oxidase subunit 4 isoform 1, mitochondrial (COX4I1), ATP synthase F1 subunit beta (ATP5F1B), Sarcoplasmic/ Endoplasmic Reticulum Calcium ATPase1 (SERCA1) and Adenosine Monophosphate Deaminase1 (AMPD1) in yak might improve oxygen utilization and energy metabolism. Pyruvate Dehydrogenase protein X component (PDHX) and Acetyltransferase component of pyruvate dehydrogenase complex (DLAT) showed higher expression levels and L-lactate dehydrogenase A chain (LDHA) showed lower expression level in yak, which might help yak reduce the accumulation of lactic acid. In addition, higher expression levels of Filamin C (FLNC) and low levels of AHNAK and Four and a half LIM domains 1 (FHL1) in yak might reduce the risks of pulmonary arteries vasoconstriction, remodeling and hypertension.

Conclusion:

Overall, the present study reported the differences in protein profile between yak and cattle, which might be helpful to further understand molecular mechanisms underlying yak adaptation to high altitude environments.

NOD-like receptors mediate inflammatory lung injury during plateau hypoxia exposure

Background

The lung is an important target organ for hypoxia treatment, and hypoxia can induce several diseases in the body.

Methods

We performed transcriptome sequencing for the lungs of rats exposed to plateau hypoxia at 0 day and 28 days. Sequencing libraries were constructed, and enrichment analysis of the differentially expressed genes (DEGs) was implemented using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Subsequently, experimental validation was executed by quantitative real-time PCR (qRT-PCR) and western blot.

Results

The results showed that the nucleotide-binding oligomerization domain (NOD)-like receptor (NLR) signaling pathway that was involved in immunity may play a crucial function in lung injury caused by plateau hypoxia. And the expressions of NOD1, NOD2, IL-1β, TNF-α, IL-6, and IL-18 were higher at 28 days of exposure to plateau hypoxia than that at 0 day. Similarly, CARD9, MYD88, p38 MAPK, and NF-κB p65, which are related to the NF-κB and MAPK signaling pathways, also demonstrated increased expression at 28 days exposure to plateau hypoxia than at 0 day.

Conclusions

Our study suggested that the NF­κBp65 and p38 MAPK signaling pathways may be activated in the lungs of rats during plateau hypoxia. Upregulated expression of NF­κBp65 and p38 MAPK can promote the transcription of downstream inflammatory factors, thereby aggravating the occurrence and development of lung tissue remodeling.

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.

Genome-wide selective sweep analysis of the high-altitude adaptability of yaks by using the copy number variant

The domestic yak (Bos grunniens) from the Qinghai-Tibet Plateau is an important animal model in high-altitude adaptation studies. Here, we performed the genome-wide selective sweep analysis to identify the candidate copy number variation (CNV) for the high-altitude adaptation of yaks. A total of 531 autosomal CNVs were determined from 29 yak genome-wide resequencing data (15 high- and 14 low-altitude distributions) by using a CNV caller with a CNV identification interval > 5 kb, CNV silhouette score > 0.7, and minimum allele frequency > 0.05. Most high-frequency CNVs were located at the exonic (44.63%) and intergenic (46.52%) regions. In accordance with the results of the selective sweep analysis, 7 candidate CNVs were identified from the interaction of the top 20 CNVs with highest divergence from the F _ST and V _ST between the low (LA) and high (HA) altitudes. Five genes (i.e., GRIK4, IFNLR1, LOC102275985, GRHL3, and LOC102275713) were also annotated from the seven candidate CNVs and their upstream and downstream ranges at 300 kb. GRIK4, IFNLR1, and LOC102275985 were enriched in five known signal pathways, namely, glutamatergic synapse, JAK-STAT signaling pathway, cytokine-cytokine receptor interaction, neuroactive ligand-receptor interaction, and olfactory transduction. These pathways are involved in the environmental adaptability and various physiological functions of animals, especially the physiological regulation under a hypoxic environment. The results of this study advanced the understanding of CNV as an important genomic structure variant type that contributes to HA adaptation and helped further explain the molecular mechanisms underlying the altitude adaptability of yaks.

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.

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

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.